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FOREWORD 


The  research  work  in  this  report  was  performed  by  Western  Electro-Acoustic 
Laboratory,  Inc.,  Los  Angeles,  California,  for  the  Flight  Dynamics  Laboratory,  Directorate 
of  Aeromechanics,  Deputy  for  Technology,  Aeronautical  Systems  Division,  Wright-Patterson 
Air  Force  Base,  Ohio,  under  AF  Contract  Nr  AF33(616)-7772.  This  research  is  part  of  a 
continuing  effort  to  obtain  more  accurate  prediction  methods  for  noise  suppression  in  Jet 
f'i-'W'  for  flight  vehicles  which  is  part  of  the  Air  Force  Systems  Command's  Applied  Research 
Program  750A,  the  Mechanics  of  Flight.  The  Project  Nr  is  1370  "Dynamic  Problems  in 
Fl'  jht  Vehicles"  and  the  Task  Nr  is  137005  "Methods  of  Noise  Prediction,  Control  and 
Measurement",  David  Smith  and  later  Phillip  Hermes  of  the  Flight  Dynamics  Laboratory 
wvre  the  Project  Engineers,  and  their  helpful  comments  have  been  appreciated  by  the  authors. 
The  research  was  conducted  from  December  1960  to  January  1962. 

The  individual  participants  in  this  report  include:  R.  Kosin  of  Norair,  who  acted  as 
a  consultant  to  the  program  and  contributed  the  major  part  of  Section  V  entitled  "Performance 
Aspects  of  Noise  Suppression  Nozzles  and  Turbofan  Engines";  R.  White,  who  wrote  Appendix 
B  entitled  "Derivation  of  Modified  Acoustic  Equation  for  Turbulent  Jet  Exhausts";  M.  Cottis, 
who  wrote  Appendix  C  entitled  "Sound  Propagation  in  a  Cylindrically  Symmetrical  Tempera¬ 
ture  Distribution";  M.  Mann,  who  derived  the  near  field  prediction  model  in  Section  III;  and 
K.  Eldred,  who  is  primarily  responsible  for  the  remainder,  including  Appendix  A. 
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ABSTRACT 


This  report  gives  methods  for  computing  the  acoustic  power  spectra  generated  by  jets 
from  coaxial  turbofans,  generalized  mixing  nozzles  and  ejectors,  and  compares  the  aero¬ 
dynamic  performance  gains  or  penalties  associated  with  these  various  devices.  The  results 
give  the  limits  of  acoustic  power  reduction  which  can  be  achieved  by  mixing  nozzles. 
Further,  the  results  show  that  major  new  reduction  in  jet  noise  can  be  obtained  only  with 
the  development  of  new  turbofan  engines  which  have  the  maximum  practical  bypass  ratios. 

These  prediction  methods  for  noise  suppression  in  jet  flows  are  based  in  part  on  deriva¬ 
tions  of  the  flow  from  axisymmetric  and  slot  jets  contained  in  the  report.  These  flow 
derivations  are  based  on  the  Prandtl  mixing  length  concept  and  include  an  analysis  of  the 
mixing  length  determined  from  measured  jet  turbulence,  and  the  interrelationship  between 
flow  Mach  number,  mixing  length,  and  the  resultant  jet  core  length. 

The  report  also  derives  a  new  model  of  the  angular  distribution  of  the  radiation  of 
noise  from  a  1000°F  Mach  1  jet  as  a  function  of  both  frequency  and  axial  position  from 
empirical  data.  This  derivation  sheds  considerable  new  light  on  the  radiation  properties  and 
generation  of  noise  in  the  jet.  These  results  are  found  consistent  with  a  new  basic  differen¬ 
tial  equation  governing  the  radiation  and  generation  of  sound  within  the  jet,  which  is 
derived  in  this  report.  In  addition,  theoretical  consideration  is  given  to  the  radiation  of 
sound  in  a  cylindrical  temperature  gradient  and  to  the  acoustic  power  generation  of  the  jet 
as  a  function  of  axial  position. 
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SECTION  I 

INTRODUCTION 


The  work  described  in  this  report  represents  a  continuation  of  efforts  by  this  and 
other  organizations  toward  practical  understanding  and  reduction  of  jet  noise. 

The  primary  purpose  of  the  present  study  is  the  evaluation  of  methods  to  reduce  jet 
noise  in  the  near  field  where  it  is  responsible  for  the  sonic  fatigue  of  aircraft  structure.  The 
validity  and  practicality  of  these  evaluations  presuppose  an  adequate  understanding  of  the 
noise  generating  and  radiating  characteristics  of  arbitrary  jet  flows.  In  turn,  this  under¬ 
standing  depends  upon  a  satisfactory  definition  of  the  jet  flow  itself,  and  upon  the  relation¬ 
ships  between  flow  parameters  and  noise  generation. 

Unfortunately,  these  interrelationships  are  not  well  documented.  Therefore,  consid¬ 
erable  emphasis  in  this  study  was  directed  toward  the  development  of  similarity  parameters  for 
jet  flow,  near  field  noise  generation,  and  refraction  of  noise  within  the  jet.  These  efforts 
assisted  in  the  semi-empirical  determination  of  the  spacial  and  directional  distribution  of 
noise  along  the  outer  boundary  of  a  conventional  jet  at  military  power.  In  addition,  they 
provided  a  basis  for  an  evaluation  of  the  effect  of  flow  alteration  on  jet  noise  generation, 
supplemented  with  a  discussion  of  the  aerodynamic  considerations  implicit  in  multiple  nozzles 
and  bypass  engines. 

These  topics  are  developed  throughout  the  report  as  follows:  Section  II  relates  jet 
flow  and  turbulence  similarity  parameters.  Section  ill  documents  the  near  noise  field  of  the 
1000°F,  Mach  1 ,  axisymmetric  jet  and  derives  an  empirical  model  for  the  angular  distribution 
of  noise  as  a  function  of  frequency  and  axial  position.  The  performance  of  various  noise 
control  devices  is  treated  in  Section  IV  comparing  measured  results  with  those  expected  on 
the  basis  of  prediction  techniques  developed  in  Sections  II,  III  and  IV.  Finally,  the  aero¬ 
dynamic  performance  penalties  or  gains  of  noise  suppression  nozzles  and  turbofan  engines  are 
discussed  in  Section  V.  Three  appendixes  follow  the  text.  Appendix  A  derives  solutions  for 
several  representative  jet  flows.  A  theory  of  generation  and  radiation  of  sound  is  discussed 
in  Appendix  B  and  a  new  equation  is  derived  which  removes  some  of  the  restrictions  of  earlier 
theory.  The  properties  of  sound  radiating  from  the  axis  of  a  cylindrically  symmetric  tempera¬ 
ture  distribution  are  developed  in  Appendix  C. 


Manuscript  released  by  the  authors  September  1962  for  publication 
os  an  ASD  Technical  Documentary  Report. 
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SECTION  II 

TURBULENCE  AND  FLOW  SIMILARITY 


The  noise  generated  by  a  jet  flow  results  fiw-i  the  turbulent  fluctuations  within 
the  flow.  The  magnitude  of  this  noise  is  dependent  upon  the  intensity  and  frequency  of 
the  turbulent  fluctuations  and  the  size  of  fluctuating  eddies.  Further,  the  magnitude  of 
the  radiated  noise  is  dependent  upon  the  mean  flow  parameters  and  their  gradients.  Thus, 
alteration  of  one  or  more  of  the  basic  flow  characteristics  may  be  expected  to  alter  the 
characteristics  of  the  noise  radiated  from  the  flow.  This  is  the  basic  premise  of  the  in¬ 
flight  noise  suppressor. 


The  first  step  toward  the  evaluation  of  the  noise  radiated  from  a  modified  jet  flow  is 
the  estimation  of  the  characteristics  of  the  flow  itself.  Since  no  purely  analytical  theory 
or  method  is  known  which  can  predict  the  noise  radiated  from  a  jet  flow,  the  evaluation  of 
the  noise  radiation  of  a  modified  flow  is  heavily  dependent  upon  analogy  to  the  noise 
radiation  of  we  1 1 -documented  flows.  In  this  latter  category,  the  majority  of  available 
data  applies  to  the  axisymmetric  jet,  and  therefore  it  provides  a  primary  basis  for  evalua¬ 
tion  of  the  effects  of  flow  modification. 


Figures  1  to  3  illustrate  the  gross  flow  parameters  for  an  axisymmetric  jet.  As  shown 
in  figure  1,  the  jet  core  extends  for  several  nozzle  diameters  from  the  nozzle.  This 
essentially  laminar  core  is  surrounded  by  a  mixing  region  which  grows  linearly  with  axial 
distance.  The  maximum  turbulent  intensity  occurs  in  the  portion  of  the  mixing  region 
where  the  velocity  gradient  is  steepest,  as  shown  in  figure  2.  Since  the  noise  generated 
in  the  flow  is  directly  dependent  upon  a  power  of  the  turbulent  intensity,  the  region  wherein 
the  mean  square  turbulence  (u1)2  exceeds  one-half  of  its  maximum  value  affords  a  concep¬ 
tual  definition  of  the  actual  noise  source  region.  Thus,  the  region  of  primary  noise 
generation  in  an  axisymmetric  jet  is  an  annular  volume  which  grows  in  thickness  with 
increased  axial  distance,  as  illustrated  in  figure  3. 


Downstream  of  the  core  tip  the  axial  velocity  decreases  and  the  jet  continues  to 
spread.  This  decrease  in  velocity  reduces  the  velocity  gradients  and  the  intensity  of 
turbulence  generated  in  the  downstream  region.  Here,  as  pointed  out  by  Ribner  (Ref.  1), 
the  intensity  of  the  noise  generated  downstream  of  the  core  should  decrease  with  a  rela¬ 
tively  high  power  of  the  axial  distance.  Therefore,  as  will  be  shown  later  in  a  more 
conclusive  fashion,  the  mixing  region  external  to  the  core  is  the  major  source  of  the 
noise  generated  by  the  jet  flow. 


This  conceptual  model  of  the  relationships  between  turbulence  and  gross  flow 
parameters  is  continued  in  the  following  subsections  to  provide  a  practical  definition  of 
the  similarity  parameters  required  for  later  consideration  of  the  noise  generation  of  modified 
flows. 


Note: 


The  references  are  numbered  consecutively  in  each  section  and  are  listed  at  the  end 
of  each  section  together  with  a  list  of  symbols. 
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Figure  1 .  Mean  flow  for  constant  density  axisymmetric  |et  at  Mach  .7  from  Laurence 
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Jet  Turbulence 


Although  turbulent  phenomena  have  been  observed  in  both  wind  and  sea  throughout 
history,  a  great  deal  remains  to  be  understood  about  the  basic  mechanisms  of  turbulence, 
particularly  for  shear  flows.  Therefore,  it  is  necessary  to  rely  primarily  on  experimental 
results  with  only  general  theoretical  guidance.  Fortunately,  considerable  turbulence  data 
have  been  obtained  in  jet  flows  by  Corrsin  (Refs.  2,  3,  41,  Liepman  (Ref.  5),  Laurence 
(Refs.  6,  7,  8)  and  others.  Since  Laurence's  data  were  token  with  primary  regard  to  noise 
generation,  the  major  emphasis  in  his  experiments  was  placed  on  the  mixing  region  surround¬ 
ing  the  core.  The  resulting  data  form  the  most  complete  set  of  flow  and  turbulence  data 
available  for  the  study  of  jet  noise. 


The  longitudinal  velocity  spectra  for  an  axisymmetric  constant  density  jet  of  1.75 
in.  radius  and  exit  Mach  number  of  .4  measured  by  Laurence  (Ref.  6)  are  given  in 
figure  4, 


where 


J  (j  t  Af-*o  Af 


Af 


u2  is  the  mean  square  turbulence  velocity  measured  in  a 


bandwidth  of  Af(cps) 


and 


u0Qis  the  total  mean  square  turbulence  over  the  total 
^bandwidth 

/  if  Cf Hf  =  I 


The  spectra  for  the  non-dimensional  axial  distances  (x/re)  between  2.29  and  16 
have  been  normalized  by  Ue/x  and  are  given  as  a  function  of  a  Strouhal  number  fx/Ue  in 


figure  5.  Despite  the  unique  low  frequency  characteristic  of  the  spectra  at  x/re  of  2.29 


and  4.58  with  respect  to  the  spectra  farther  downstream,  there  is  an  apparent  gross  simi¬ 
larity  of  the  normalized  spectra  throughout  the  entire  mixing  region. 


The  lateral  turbulence  spectra  (Ref.  6),  measured  at  two  stations,  have  been  normal¬ 
ized  by  U e/x  and  are  compared  to  the  longitudinal  spectra  in  figure  6.  These  data 
indicate  that  the  lateral  velocity  spectra  are  probably  also  similar;  however,  they  appear 
shifted  in  frequency  by  a  factor  of  two  above  the  longitudinal  spectra. 


Although  the  data  shown  in  figures  4  to  6  offer  evidence  of  gross  similarity  for  the 
turbulence  spectra  as  a  function  of  axial  distance,  they  do  not  test  similarity  as  a  function 
of  velocity  because  Ue  was  not  varied  in  the  experiment.  Furthermore,  these  turbulence 
spectra  do  not  give  the  actual  frequency  spectra  of  the  turbulent  fluctuations  themselves, 
when  viewed  by  an  observer  moving  at  the  mean  flow  velocity.  For,  as  Taylor  pointed  out, 
(Ref.  9)  and  Laurence  (Ref.  6)  proved  for  the  majority  of  these  data,  these  turbulent  power 
spectra  density  curves  are  simply  the  Fourier  transforms  of  the  longitudinal  spacial  velocity 
correlation  coefficient  (H.x)«  Therefore,  the  turbulent  spectra  measured  with  the  hot  wire 
are  exactly  a  measure  of  the  mean  square  fluctuating  velocities  associated  with  eddies  of 
varying  size  ( ^ ),  which  are  convected  past  the  measuring  device  with  velocity  Uc.  Thus, 


J 

■J 


the  frequency  scale  In  figure  4  Is  equivalent  to  Uc/X/  and  the  similarity  found  In 
figures  5  and  6  really  applies  to  the  distribution  of  turbulent  energy  with  respect  to  eddy 
size/  rather  than  with  respect  to  actual  turbulent  frequencies. 

The  rate  of  the  spread  of  the  mixing  region  in  the  flow  Is  governed  by  the  turbulent 
shear  stress  which,  for  an  axisymmetric  flow.  Is  given  primarily  bypuv  where  u  and 
v  are  the  instantaneous  velocity  fuctuations  in  the  longitudinal  and  lateral  directions, 
respectively.  There  is  considerable  evidence  in  Anderson  and  Johns'  work  (Ref.  10) 
that  the  length  of  the  core  of  an  axially  symmetric  jet  increases  significantly  at  higher 
Mach  numbers,  implying  a  reduction T  with  an  increase  of  Mach  number.  The  relation¬ 
ship  between and  M  is  also  significant  for  noise  generation  since  rhe  primary  noise 
source  term  Is  the  second  time  derivative  of't'.  Unfortunately,  no  measurements  of"T 
itself  could  be  found  for  flows  with  nozzle  Mach  numbers  greater  than  approximately  .2. 
However,  Laurence  (Ref.  6)  gives  profiles  of  the  longitudinal  velocity  fluctuations  over 
a  Mach  range  from  .2  to  .7  which  show  a  decrease  in  the  ratio  of  rms  turbulent  velocity 
to  exit  velocity  (u'/Ue)  as  Mach  number  is  increased.  These  data  afford  an  opportunity 
to  determine  the  approximation  variation  q'  with  Mach  number. 


Several  investigators,  including  Prandtl,  Taylor  and  vonKnrmen,  have  proposed 
expressions  which  relate  the  turbulent  shear  stress  to  the  gross  flow  parameters  (Ref.  11). 
0(  these  expressions,  the  concept  of  a  constant  exchange  coefficient  E  across  the  mix- 


dU 


ing  regior,  such  that  *X  /  and  the  original  Prandtl  mixing  length  hypotheses. 


are  most  readily  applied  to  jet  flows. 


Figure  7.  Illustration  of  Prandtl  mixing  length  hypothesis. 


The  mixture  length  (i  )  in  Prandtl's  hypothesis  is  analogous  to  the  mean  free  path 
in  the  kinetic  theory  of  gasses.  If  a  small  mass  of  fluid  initially  at  radius  r],  and 
velocity  U],  as  illustrated  in  figure  7,  is  displaced  radially  outward  to  r,  it  gives 
rise  to  an  increase  of  velocity  at  T,  the  magnitude  of  the  increase  Au  being  equal  to 
U]  -  U. .  Similarly/  a  small  mass  of  fluid  at  r 2  displaced  inward  to  r  results  in  a 
decrease  in  the  velocity  at  r  of  magnitude  0  -  U2*  With  this  assumption/  the  rms 
fluctuating  longitudinal  velocity  component  u'  at  r  is: 

uWl/2(|AU]|  + 


where  AU]  =  U]  -  U  =  “U  (  r  -  £  )  -  VL(  r  )  =  JdU  +  higher  order 

dr 

terms. 

Taking  the  lateral  velocity  component  proportional  to  the  longitudinal  component/ 
the  shear  stress  becomes: 


dU 

ar 


dUl  dU 


Tr\Tr 


Jar  "W 

where  the  absolute  value  of  ^  is  utilized  to  insure  that  “X  has  the  same  sign 

dr 

as  dU 
Tr 

Defining  a  non-dimensional  mixing  length  parameter  £  ,  where 

-l- 

re 


*  f  (r/r  e,  x/r  e,  M,  Re) 


<J7 


the  Laurence  data  were  examined  to  determine  the  functional  dependence  of  ^  on  radial 
and  axial  position,  Mach  number  and  Reynolds  number. 

Figure  8  gives  £  /§  0v  a*  a  function  of  radial  position  and  at  various  axial  distances 
where  5  av  I*  the  average  value  of  §  across  the  central  mixing  zone  at  each  station  and 
for  each  Mach  number.  Superimposed  on  the  scatter  of  approximately  ±10%  appears  to  be 
a  systematic  variation  in  the  mean  line  with  the  value  of  §  at  the  center  of  the  mixing 
zone  approximately  4%  above  its  average  value. 

The  variation  of  5  av  with  QX^a*  position  is  given  in  figure  9.  The  data  clearly  show 
that  §  av  is  a  linear  function  of  axial  distance  and  o  function  of  Moch  and/or  Reynolds 
number.  Although  it  is  not  possible  to  separate  the  latter  effects  in  the  present  data 
because  only  one  si*^  iet  was  utilized,  the  Reynolds  number,  which  wos  above  192,000  for 
the  various  experiments,  is  sufficiently  high  to  expect  that  fully  turbulent  flow  develops  in 
the  mixing  region  in  a  relatively  short  distance  downstream  from  the  nozzle.  Therefore, 
it  Is  probable  that  the  major  factor  is  Mach  number.  The  magnitude  of  the  dependence  of 
^aV  on  Moch  number  is  shown  in  figure  10,  where  ^  ^  ov  ^  is  given  os  a  function 

*Ae 

of  Mach  number.  The  relationship  between  the  lateral  scale  of  turbulence  found  by 
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Laurence  in  his  spacial  correlations  and  the  defined  by  the  mixing  length  hypothesis 
Is  shown  in  figure  11.  At  Mach  .3,  H  is  approximately  2/3  of  the  conventional  lateral 


scale  of  turbulence  ^  where 


^=Jjh  ^  andT^is  the  correlation  of  longitudinal  velocity  fluctuations 
as  a  function  of  lateral  distance  (tj)  between  the  measurement  points. 


To  the  extent  that  the  mixing  length  hypothesis  represents  the  gross  turbulent  mix¬ 
ing  process  in  the  region  of  maximum  jet  shear,  the  frequency  of  the  turbulence  (f  )  in 
die  coordinate  system  moving  at  the  mean  eddy  convection  velocity  is  of  the  order  of 


f  aX  =  1  dU 

fl  23T 


Figure  11.  Comparison  of  lateral  scale  of  turbuience  at  Mach  .3  with  mixing  length  as 
a  function  of  axial  distance,  data  from  Laurence. 
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A*  will  be  shown  later,  the  maximum  velocity  gradient  at  the  tip  of  the  core  Is  equal  to 
1 .2  Ue  where  d  is  the  jet  exit  diameter.  Hence,  with  these  assumptions,  the  predicted 

T  U 

center  frequency  of  the  turbulence  at  the  core  tip  would  be  of  the  order  of  .6  ®  ,  or 

a 

equivalently,  the  Strouhal  number  based  on  exit  diameter  and  velocity  would  be  approxi¬ 
mately  .6  at  the  core  tip. 


Jet  Flow  Similarity 

The  flow  of  a  fully  turbulent  axisymmetric  and  a  two-dimensional  slot  Jet  with  zero 
external  velocity  (Uc  =  0),  and  constant  density  has  been  investigated  theoretically  by 
several  individuals  (Refs,  li,  12).  More  recently,  Szablewski  (Ref.  13)  undertook  solu¬ 
tion  of  the  downstream  fully  developed  flow  for  variable  jet  density  and  UQ>  0. 


The  flow  in  the  mixing  region  adjacent  to  the  core  was  first  solved  by  Kuethe  (Ref.  14) 
for  a  constant  density  jet  with  Uo  =  0.  His  solution  contained  an  asymptotic  velocity 
profile  obtained  through  an  iterative  procedure.  Later,  Squire  and  Trouncer  (Ref.  15) 
obtained  solutions  for  the  complete  constant  density  jet  with  variable  external  velocity 
utilizing  a  (1  -  cosine)  velocity  profile.  For  the  purposes  of  this  study,  it  was  desirable  to 
have  a  uniform  set  of  flow  solutions  which  included  both  axisymmetric  and  slot  jets,  and 
which  allows  extension  to  high  Mach  number  and/or  variable  density  flows.  It  was  also 
desirable  that  the  solutions  utilize  an  asymptotic  profile  which  could  be  readily  adapted  to 
concentric  and  more  complicated  jet  flows.  Further,  the  profile  should  match  the  mean 
flow  data  and,  where  possible,  agree  with  the  profile  of  turbulent  fluctuations  so  that  it 
would  have  general  utility  in  describing  the  turbulence  in  the  mixing  region  in  terms  of 
noise  generation. 
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An  examination  of  the  mean  flow  profiles  in  the  region  adjacent  to  the  core  shows 
that  they  can  be  represented  by  an  exponential  of  the  form  whereat  is  a  non- 

dimensional  shape  parameter.  The  maximum  error  occurs  at  the  outermost  boundary  of  the 
flow  and  Is  of  the  order  of  4%.  Since  the  low  value  of  shear  at  the  outermost  boundary  is 
not  important  to  the  overall  mixing  process,  this  error  is  acceptable.  This  profile  function 
has  been  previously  utilized  In  the  fully  developed  downstream  region  by  Reichardt  (Ref.  16) 
and  others,  so  that  Its  use  enables  a  simple  approximate  extension  from  the  core  region  to 
the  entire  flow. 

The  flow  for  a  constant  density  axisymmetric  jet  with  variable  external  velocity, 
and  a  constant  density  slot  jet  has  been  derived  in  Appendix  A.  The  derivations  utilize 
the  integral  form  of  the  momentum  equations  described  by  Goldstein  (Ref.  17)  and  follow 
the  general  approach  of  Squire  and  Trouncer  (Ref.  1^5).  In  the  derivations  the  velocity 
profile  was  assumed  to  be  U  =  UQ  +  (Un,  -  U0)e“  ^  /* 

where  U0  is  the  external  velocity, 

U(n  is  the  centerline  velocity,  and  equal  to  Ue  adjacent  to  the 
core  region 

and  =  f  ~  o  where  a  is  the  distance  from  the  centerline  to  the  outer 
"T>  boundary  of  the  core  and  b  is  a  width  parameter  which 
controls  the  magnitude  of  the  velocity  gradient. 
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Figure  12.  Comparison  of  assumed  velocity  profile  for  axisymmetric  constant 
density  jet  with  Mach  .7  data  from  Laurence. 


This  profile  is  compared  with  Laurence's  (Ref.  6)  Mach  .7  data  at  various  axial 
stations  in  figure  12.  Note  that  when  r  -  a  =  b,  1,  and  U/(Um-U0)  =  .606. 


The  Appendix  A  solution  for  relative  length  of  the  core  (Xi )  of  the  axisymmetric  jet 


is 


xt  _  ^ 

—  =  k  c2 

where  c  =  and  k  is  a  constant  =  .88. 


It  was  noted  previously  that  the  mixing  length  X  varied  as  a  function  of  Mach  number 
in  the  range  of  M  =  .2  to  .7.  Further,  the  data  of  Anderson  and  Johns  (Ref.  10)  have 
demonstrated  that  the  relative  core  length  becomes  much  longer  as  the  exit  Mach  number 
is  increased.  These  results  have  been  combined  in  figure  13  to  give  the  variation  of  c 
with  Mach  number  (from  Refs.  3,  10,  18,  19),  Figure  14  shows  the  equivalent  variation  in 
relative  core  length  with  exit  Mach  number.  Perhaps  the  most  interesting  feature  of  these 
two  curves  is  the  small  effect  which  large  changes  in  temperature  have  on  the  core  length, 
relative  to  the  significant  role  of  Mach  number.  In  general,  the  hot,  low  density,  jet  core 
appears  to  be  no  more  than  20%  shorter  than  the  core  of  a  jet  of  the  same  exit  Mach  number 
which  has  a  total  temperature  equal  to  ambient. 
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Several  types  of  theoretical  derivations  of  the  variable  density  jet  were  conducted 
under  this  program.  However,  the  results  were  inconclusive  in  that  they  were  excessively 
dependent  upon  initial  assumptions,  particularly  in  regard  to  the  ratio  of  heat  to  momentum 
transfer  which  has  been  variously  reported  (Ref.  1 1)  through  a  range  of  1.3  -  2.0.  These 
values  result  from  previous  theoretical  and  experimental  efforts  which  were  restricted  to 
temperature  ratios  (jet  to  ambient)  of  approximately  2  and  in  general  are  limited  to  low 
Mach  number.  From  the  present  investigation  it  appears  that  the  theoretical  solution  of 
the  effects  of  density  on  jet  flow  in  the  primary  range  of  interest  requires  additional 
experimental  effort  to  gain  a  better  definition  of  the  relationship  between  heat  and 
momentum  transfer  at  high  Mach  number.  However,  for  the  present,  the  indication  that  a 
hot  jet  core  length  is  reduced  no  more  than  20%  compared  with  the  ambient  jet,  is  suffi¬ 
cient. 


The  flow  contours  calculated  from  the  Appendix  A  solution  for  an  axisymmetric  jet 
with  an  exit  Mach  number  of  1  are  given  in  figure  15.  The  solid  contours  give  the  ratio 


U 

of  tt —  and  the  dashed  contours  downstream  of  the  core  give  the  ratio  of  -j-j — ,  where 
ue  um 

Um  Is  the  centerline  velocity.  These  contours  can  be  extended  to  any  other  Mach  number 
by  revising  the  axial  scale  by  the  ratio  of  xt  where  xf  is  obtained  for  the  new  Mach 
number  from  figure  14.  13 


U 


Figure  16  shows  that  the  calculated  axial  decay  of  the  mean  velocity  In  the  axi¬ 
symmetric  jet  downstream  of  the  nozzle  agrees  with  the  data  from  Laurence  (Ref.  6).  In 
addition,  figure  16  shows  the  relative  turbulent  velocities  u'm  from  Laurence  as  a 


u  fm 


function  of  axial  position.  Here  u^  is  the  maximum  value  of  rms  turbulent  velocity  at 
each  axial  station  and  u'fm  Is  the  maximum  value  of  u'm  in  the  entire  flow.  From  the 
simple  flow  theory,  the  axial  variation  of  u'm  should  be  similar  to  the  variation  of  center- 


u'fm 


line  velocity.  This  relationship  appears  to  hold  to  x  of  24.  However,  farther  down¬ 


stream  the  turbulence  is  somewhat  greater  than  would  be  expected  from  the  local  value  of 
the  flow  gradients.  These  higher  values  of  u'm  probably  result  from  the  convection  of 


u'fm 


the  more  intense  turbulence  generated  upstream  past  the  downstream  positions. 


Figure  17  gives  the  calculated  flow  contours  for  an  axisymmetric  jet  at  Mach  1, 
with  an  external  velocity  U0  equal  to  one-half  of  the  nozzle  exit  velocity.  The  value  of 
c  utilized  in  the  equation  was  chosen  at  Mach  .5,  the  relative  Mach  number  for  the  two 
flows.  It  Is  interesting  to  note  that  the  external  flow  has  lengthened  the  core  to  twice  the 
axial  distance  for  zero  external  velocity. 


This  result  Is  generalized  in  figure  18  where  the  ratio  of  Xf/(xf  for  U0  =  0)  Is  given 
as  a  function  of  U0.  Since  this  curve  depends  upon  Mach  number,  through  the  dependence 
of  c  on  Mach  number,  actual  computations  of  ewe  length  should  be  based  directly  upon  the 


17 


VV1  VLN  V*  V".  \%  V 


-  * 


•.'  s'  V 


31 


■jj— =  0  for  U0  -  0 


1  for  Uo=0 


10  12  14  16  18  20  22  24  26 

x/re 


Figure  17.  Calculated  flow  contours  for  constant  density  axisymmetric  jet  with 
external  velocity  UQ  =  .5  Ue/  Mg  =  1,  and  using  c  based  on  relative  Machinumber. 


Velocity  ratio  U^/t^ 

Figure  18.  Calculated  variation  in  the  core  length  as  a  function  of  the 
ratio  of  external  velocity  (UQ)  to  jet  exit  velocity  {Ue)(not  including 
variation  of  c  with  Mach  number.) 
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Figure  20.  Axial  decay  of  velocity  computed  for  Mach  .3  slot  jet  using 
(k  c2)  from  axisymmetric  jet.  Data  are  for  3. 14  x  1  in.  nozzles  with 
two  spacings  from  Laurence,  NACA  TN  4029. 
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LIST  OF  SYMBOLS  USED  IN  SECTION 
Roman 


a 

b 


d 

E 

h 


k 

£ 


Xi 

M 


■R. 

Re 


r 

U 

u 

v 

X 

V 
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A 
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c 

e 

f 

m 

o 

oa 

t 

x 
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Distance  from  centerline  to  boundary  of  core 
Velocity  gradient  width  parameter 

XA 

Diameter 

Momentum  exchange  coefficient 

Normalized  longitudinal  velocity  spectral  density 

Frequency 

Thickness  of  slot  jet  flow 
Proportionality  coefficient 
Prandtl  mixing  length 
Lateral  scale  of  turbulence 
Mach  number 

Correlation  of  longitudinal  velocity  fluctuations 

Reynolds  number 

Radius 

Flow  velocity 

Instantaneous  longitudinal  velocity  fluctuations 
Instantaneous  lateral  velocity  fluctuations 
Axial  distance  from  exit  plane 
Lateral  distance  from  centerline 


Greek 


Non-dimensional  mixing  length 
G;aracteri$tic  eddy  size 
Non-dimensional  shape  parameter 
Density 

Turbulent  shear  stress 


Subscripts 


Characteristic 

Exit 

Flow 

Maximum  at  a  given  axial  station 

Ambient 

Overall  level 

Core  tip 

As  a  function  of  longitudinal  distance  between  measurement  points 
As  a  function  of  lateral  distance  between  measurement  points 


Superscripts 


Mean 

Root  mean  square 
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SECTION  III 


THE  NEAR  NOISE  FIELD 

The  near  noise  field  of  a  jet  is  very  complex  in  comparison  to  the  far  noise  field, 
since  the  spacial  distribution  of  the  noise  sources  in  the  jet  and  other  possible  near  field 
effects  must  be  considered.  Furthermore,  since  the  measurements  required  for  definition 
of  the  near  noise  field  are  extensive,  only  a  few  relatively  complete  investigations  of  a 
full  scale  jet  engine  have  been  undertaken  (Refs.  1,  2,  3).  Similar  detailed  measure¬ 
ments  were  undertaken  for  small  rockets  (Ref.  4)  and  model  cold  jets  (Refs.  5,  6),  but  the 
major  volume  of  other  data  comes  from  more  limited  measurements  on  a  variety  of  types  of 
jets. 

Several  efforts  (Refs.  2,  6,  7,  8)  have  been  made  to  generalize  some  of  the  existing 
near  field  data  for  use  as  a  prediction  model.  References  7  and  8  presented  simplified 
empirical  models  for  the  noise  radiated  forward  from  rocket  flows.  In  addition.  Ref.  7 
presented  a  method  for  rotating  the  generalized  noise  contours  for  a  sonic  1000*F  jet,  to 
account  for  higher  jet  temperatures  and  to  predict  aircraft  noise  environments.  Wolfe 
(Ref.  2)  presented  the  spacial  variation  of  the  exponent  (n),  from  assuming  sound  pressure 
at  a  point  in  the  near  field  proportional  to  U'  where  U  is  the  velocity,  for  a  series  of  hot 
jet  velocities  ranging  from  790  -  1800  ft/sec.  Howes  (Raf,  6)  applied  and  extended  similar 
concepts,  which  he  ascribes  to  Greatrex  (Ref.  9),  in  the  analysis  of  detailed  subsonic  cold 
jet  flows. 

All  of  the  above  efforts  to  determine  useful  prediction  models  are  predicated  on  the 
similarity  of  the  near  noise  field  from  an  axisymmetric  jet.  In  this  manner.  It  Is  unnecessary 
to  postulate  the  mechanism  of  noise  generation,  radiation,  and  source  distribution.  All 
that  is  required  Is  a  demonstration  of  similarity  of  the  noise  field  itself.  Since  various 
investigators,  including  Morgan,  et  al,  (Ref.  10),  have  shown  that  the  near  noise  field  of 
a  jet  scales  directly  with  jet  size,  as  long  os  the  flow  parameters  are  held  constant,  it  is 
only  necessary  to  determine  the  change  in  the  noise  field  resulting  from  density,  tempera¬ 
ture,  and  velocity  changes  in  the  jet.  For  the  majority  of  practical  applications  to  aircraft, 
the  jet  exit  velocity  (Ue)  is  sonic;  hence,  the  temperature,  density  and  velocity  are  specifi¬ 
cally  related,  and  the  variation  of  a  single  velocity  exponent  is  sufficient  to  give  o  practi¬ 
cal  result.  For  the  more  general  case,  when  the  exit  velocity  is  subsonic,  the  variation  of 
the  noise  field  with  temperature  and  velocity  must  be  considered  separately. 

In  the  present  study,  the  objective  Is  the  prediction  of  the  noise  generation  and  the 
resulting  near  field  for  modified  jet  flows.  Here,  the  flow  Is  a  variable,  violating  the  basic 
postulate  of  the  foregoing  near  noise  field  models.  Therefore,  it  is  necessary  to  determine 
the  relationships  between  the  local  jet  flow  parameters  and  the  characteristics  of  the  local 
generation  and  radiation  of  noise.  In  the  following  subsection  the  noise  radio. Ion  from  a 
axisymmetric  sonic  jet  will  be  reviewed,  and  a  generalized  mode*  of  source  generation  will 
be  developed. 
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CharOvterisHcs  of  the  Near  Noise  Field  of  a  IQOO^F  Axisymmetrlc  Sonic  Jet 

The  most  Intensive  measurements  close  to  a  full  scale  jet  were  made  by  NACA 
(Ref.  1)  In  a  region  extending  from  the  jet  boundary  to  approximately  10  feet  from  the  axis 
at  the  nozzle  end  to  approximately  30  feet  from  the  axis  at  a  distance  of  60  feet  down¬ 
stream  from  the  nozzle.  These  measurements  were  made  on  an  engine  with  9600  lb. 
thrust,  1000°F  exit  temperature,  1850  ft/sec.  exit  velocity  and  1.85  ft.  diameter. 

Because  these  data  are  the  most  complete  documentation  of  the  near  noise  field  of  any  jet 
known  to  the  author,  they  have  been  utilized  extensively  in  this  discussion,  supplemented, 
when  necessary,  by  additional  data. 

Smoothed  contours  of  equal  overall  noise  levei  for  the  1000°F  sonic  jet  are  given  tn 
figure  21,  based  on  data  from  Refs.  1*  and  3.  The  figure  shows  the  steep  noise  gradient 
forward  and  to  the  side  of  the  nozzle,  an  extended  source  region  along  the  flow,  and  a 
low  noise  gradient  in  the  direction  of  the  maximum  noise  generation.  Figure  22  gives  the 
spacial  distribution  (Ref.  1)  for  three  frequencies  illustrating  the  gross  shift  in  the  angle  of 
maximum  noise  generation  as  function  of  frequency.  As  can  be  seen,  noise  at  100  cps  is 
directed  at  a  small  angle  relative  to  the  jet  boundary  and  appears  to  originate  over  an 
extended  region  of  the  jet.  In  contrast,  th*  noise  at  2000  cps  appears  to  ortglnufe  close 
to  the  nozzle  and  Is  primarily  radiated  at  a  large  angle  from  the  jet  boundary. 

The  overall  noise  levels  measured  (Ref.  1)  along  a  line  extending  from  the  nozzle  at 
an  angle  of  10°  to  the  jet  axis  are  given  in  figure  23.  The  third  octave  spectra  for  these 
data  are  given  In  figure  24  and  the  locations  of  maximum  sound  pressure  levels  for  these 
spectra  are  presented  in  figure  25.  Fr*m  figure  23  it  is  clear  that  the  overall  noise  level 
along  the  10°  boundary  is  constant  as  a  function  of  axial  distance  until  a  short  distance 
beyond  the  tip  of  the  core,  and  from  there  on  downstream  the  noise  level  decreases  r-uidly. 
|r,  figure  24,  the  one-third  octave  band  spectra  all  tend  to  exhibit  sharp  cutoffs  beyond  the 
core  tip,  and  a  gross  overall  similarity.  This  gross  similarity  of  the  noise  spectra  in  the 
core  region  would  be  expected  from  the  similarities  previously  seen  for  the  turbulence  and 
flow  in  the  same  region.  From  these  considerations  it  would  be  expected  that  the  axial 
location  of  the  maximum  level  for  each  frequency  band  would  vary  inversely  with  axial 
distance  as  shown  in  figure  25.  The  departure  from  the  l/x  relationship  exhibited  at  the 
low  frequencies  in  figure  25  results  from  the  general  decrease  in  level  and  change  in  simi¬ 
larity  downstream  of  the  core. 

It  would  be  convenient  to  utilize  the  concept  of  figure  25  to  describe  the  most 
probable  axial  locations  for  the  sources  at  the  various  frequencies.  However,  the  extended 
axial  distribution  of  noise  in  each  frequency  band  shown  In  figure  24  indicates  that  this 
procedure  would  fail  at  positions  close  to  the  jet,  where  the  length  of  the  source  distribu¬ 
tion  is  within  an  order  of  magnitude  of  the  distance  from  the  position  to  the  jet.  Therefore, 

*AII  data  from  Ref.  1  have  been  corrected  to  rms  from  original  quasi  peak  Bruel  &  Kjaer 
readings  by  subtracting  2.5  db. 

26 


Figure  21.  Contours  of  overall  sound  pressure  level  (rms)  in  db  re  .0002  dyne/cm2  fc  J57  engine  at  military 
power,  from  unpublished  data  by  Bioacoustics  Branch  ASD,  and  from  Howes  et  ai. 
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Figure  22.  Contours  of  sound  pressure  'evel  in  db  re  .0002  dyne/cmJ  in  one-third 

octave  bandwidths  for  J57  engine  at  military  power,  data  from  Howes,  ef  al 
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Figure  25.  location  of  maximum  sound  pressure  level  in  one-third  octave 
frequency  bands  as  a  function  of  a  modified  Strouhal  number 
and  axial  distance.  Data  from  Howes,  ei  ai. 


It  is  necessary  to  include  this  axial  distribution  in  a  prediction  method  and  to  determine 
the  directional  properties  of  the  radiation  for  each  frequency  as  a  function  of  axial  posi¬ 
tion. 

It  was  anticipated  that  the  axial  distribution  of  source  power  would  be  exceedingly 
difficult  to  obtain  from  boundary  measurements  because  of  acoustic  near  field  effects  super¬ 
imposed  on  the  geometrical  near  field.  From  consideration  of  probable  acoustic  near 
field  effects  it  was  expected  that  the  sound  pressures  along  the  jet  boundary  would  include 
not  only  the  pressures  associated  with  the  radiated  acoustic  power,  but  additional  pressures 
associated  with  the  non-radiating  inductive  near  field,  discussed  in  detoil  by  Franz 
(Ref.  12).  Furthermore,  it  was  expected  that  the  low  frequency  pressures  very  close  to 
the  actual  jet  boundary  would  contain  components  of  convected  large  scale  eddies  as 
found  by  Franklin  and  Foxwell  (Ref.  13).  These  latter  pressures  from  convected  turbulence 
have  been  identified  approximately  on  figure  24. 

In  order  to  obtain  a  first  estimate  of  the  magnitude  of  the  near  inductive  field  effect 
on  the  boundary  pressure  measurements,  the  approximate  acoustic  power  was  computed  in 
each  octave  frequency  band,  based  on  the  maximum  pressures  along  the  10°  boundary  and 
the  area  of  the  truncated  10°  boundary  cone  bounded  by  the  axial  stations  where  the 
measured  pressure  was  3  db  below  maximum.  The  result  of  this  crude  approximation  is 
compared  in  figure  26  to  the  generalized  J57  far  field  power  spectrum  (Ref.  14). 

The  most  striking  feature  of  this  comparison  is  the  agreement  at  low  frequencies, 
since  it  was  expected  (Ref.  12)  that  the  low  frequency  power  computed  from  the  near  field 
would  be  10  to  15  db  above  the  far  field  result,  with  agreement  obtained  only  for  the  high 
frequencies.  This  result  indicates  that  a  measure  of  the  true  acoustic  power  can  be  ob¬ 
tained  from  the  10°  boundary  measurements  of  this  jet  if  proper  attention  is  given  to  the 
actual  boundary  area  normal  to  the  direction  of  radiation  as  determined  by  the  correlation 
measurements  for  the  various  frequencies  and  axial  locations. 

A  qualitative  estimate  of  the  exit  angle  for  noise  radiated  from  a  source  located  in 
the  mixing  region  at  the  point  of  maximum  shear  (figures  1,  2,  3)  is  seen  In  figure  27. 

This  ray  diagram  is  one  of  several  constructed  to  observe  the  qualitative  effects  of  the 
refraction  of  sound  within  the  flow  by  temperature  and  velocity  gradients.  The  rays  were 
constructed  for  radiation  from  a  simple  point  source  In  a  plane  containing  the  jet  axis, 
and  are  based  on  the  ray  velocity  (Ref.  15)  given  by  the  vector  sum  of  the  local  flow 
velocity  and  speed  of  sound  along  the  normal  to  the  phase  front.  As  is  clear  from  the 
figure,  a  very  pronounced  refraction  occurs,  with  the  result  that  a  significant  portion  of 
the  sound  Is  directed  at  an  angle  of  110°  to  115°  from  the  forward  axis.  It  Is  noted  that 
this  refraction  is  sufficiently  severe  so  as  to  dominate  the  directional  characteristic  of  the 
radiated  sound  for  a  relatively  wide  variety  of  assumed  source  directivities. 

Of  course,  refraction  estimated  from  a  ray  diagram  of  this  type  is  accurate  only  at 
high  frequencies,  or  more  generally  when  the  value  of  kb  is  large  in  comparison  with  unity, 
where  k  =  2TT f/a  and  b  is  the  gradient  width  parameter  utilized  in  Section  II.  Further,  the 


diagram  assumes  a  simple  point  source  located  in  the  center  of  the  mixing  zone,  rather 
than  a  finite  quadrupole  source.  Clearly,  the  exit  angle  would  shift  further  forward  for 
sources  located  toward  the  center  since  the  overall  temperature  change  is  greater  and, 
conversely,  the  exit  angle  would  shift  aft  for  sources  located  nearer  the  outer  boundary 
where  the  overall  temperature  change  is  less.  Thus,  in  the  region  of  high  kb  when  the 
gross  effects  from  ray  considerations  might  actually  apply  to  the  associated  small  sources 
(^/b^<l),  the  angular  distribution  of  radiated  noise  should  be  relatively  broad. 

With  this  qualitative  ray  diagram  as  background,  it  is  pertinent  to  review  the 
spaclal  correlation  measurements  (Ref.  1)  made  along  the  10°  boundary.  The  spacial 
correlation  coefficient  (gr*)  was  determined  from  the  time  average  product  of  the  pressure 
at  a  fixed  point  (a)  and  the  pressure  at  a  moving  point  (b)  divided  by  the  product  of  the 
rms  pressures  at  the  two  points.  Thus, 


For  a  sound  wave  in  free  space,  the  angle  of  propagation,  relative  to  the  base  line  for  the 
spacial  correlations,  can  be  determined  by  the  distance  between  the  origin  and  the  first 
zero  crossing  as  shown  in  figure  28. 

This  technique  was  applied  to  the  correlation  presented  in  Ref.  1.  Typical  results 
are  shown  in  figure  29  for  the  one-third  octave  frequency  band  centered  at  400  cps.  The 
primary  direction  of  propagation  determined  from  the  correlations  are  shown  as  arrows  from 
each  measurement  position  along  the  10°  boundary.  As  can  be  seen,  the  apparent  primary 
direction  of  propagation  at  the  location  near  the  nozzle  makes  only  a  small  acute  angle 
with  the  jet  boundary,  whereas  the  primary  direction  for  downstream  locations  makes  a 
considerably  greater  angle  with  the  jet.  It  should  be  noted  that  these  interpretations  of 
the  near  field  pressures  differ  from  those  of  Franklin  and  Foxw*?!!  (Ref.  13),  where  it  was 
shown  that  the  correlations  close  to  the  jet  boundary  indicated  the  convection  of  eddies 
past  the  measurement  microphone,  rather  than  the  radiation  of  sound  energy.  However, 
their  measurement  locations  were  close  to  the  core  boundary  of  a  model  cold  jet  and  a  full 
scale  engine,  both  of  which  had  exit  velocities  of  710  ft/sec.  This  exit  velocity  is  only 
38%  of  the  velocity  of  the  present  jet.  If  the  mean  square  radiated  sound  pressure  is 
assumed  to  vary  as  the  eighth  power  of  velocity,  whereas  the  direct  fluctuating  turbulent 
mean  square  pressures  (pseudo  sound)  vary  as  only  the  fourth  power  of  velocity,  this 
Increase  of  a  factor  of  2.6  in  exit  velocity  would  result  in  an  increase  in  the  ratio  of  mean 
square  sound  pressure  to  pseudo  sound  pressure  of  approximately  46.  This  large  potential 
increase  in  the  ratio  could  easily  account  for  measurement  of  true  radiated  sound  along  the 
10°  boundary  of  the  high  speed  jet,  whereas  measurements  on  the  lower  speed  jets  repre¬ 
sented  convected  turbulence.  It  should  also  be  noted  that  these  correlations  do  not  include 
the  data  below  x/^values  of  .25  where  pseudo  sound  is  suspected,  as  indicated  in  figure 
24. 
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Example  of  relation  between  zero  crossing  of 
wavelength  and  direction  of  propagation  . 


correlation  coefficient 


Igur®  29.  Contour*  of  equal  sound  pressure  level  in  5  db  increments  for  one-third  octave  band  centered  at 
400  cps,  together  with  axial  variation  in  level  along  10°  measurement  boundary  and  primary 
direction  from  special  pressure  correlations.  Data  from  Howes/  et  al. 
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Figure  30.  Variation  of  the  predominant  propagation  angle  for  various 
frequencies  from  measurements  along  the  10°  boundary  as 
a  function  of  a  non-dimensional  axial  frequency  parameter 
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distribution  of  the  power  spectrum.  Furthermore,  the  angle  of  maximum  radiation  in  the 
far  field  and  breadth  of  the  associated  directivity  pattern  are  dependent  on  frequency 
(Ref.  ,16),  with  the  angle  of  maximum  radiation  varying,  as  shown  in  figures  31  to  34 
from  Ref.  16,  between  approximately  150°  and  115°  for  center  frequencies  between  53 
and  6800  cps,  respectively.  Note  that  these  far  field  directivities  were  obtained  at 
distances  of  100  and  250  feet  from  the  |et  nozzle,  so  that  the  jet  is  assumed  to  approximate 
a  point  source.  The  far  field  variation  in,  maximum  angle  is  consistent  with  the  angular 
variation  of  the  primary  directions  along  the  10°  boundary,  since  the  major  portion  of  the 
low  frequency  power  is  radiated  from  regions  of  low  x/^  6and  the  majority  of  the  high 
frequency  power  is  radiated  from  regions  extending  to  high  values  of  x/a*  where  the 
primary  angle  becomes  110°  to  115°.  These  generalized  directivities,  originally  derived 
from  a  composite  of  hot  jet  measurements,  were  utilized  directly  through  an  iterative  pro¬ 
cedure  in  the  derivation  of  angular  distribution  of  sound  radiation  along  the  jet. 

For  convenience  in  the  iteration,  the  jet  was  divided  into  several  segments  by 
transverse  cuts,  and  the  power  from  each  segment  was  calculated  from  a  flux  formula,  as 
follows:  for  the  i^  segment  the  power  level  in  db  re  1 0~1 3  watts  is  given  by 


PWL|  =  SPLj  +  10  log  (A|  cos^i) 


l-D 


where  SPLj  Is  the  average  sound  pressure  level, 

Aj  the  radiating  surface  area,  and 

the  direction  of  the  sound  pressure  maxi¬ 
mum  relative  to  the  normal  of  Aj  at  its 
mean  lengthwise  coordinate. 

To  obtain  a  correspondence  between  the  for  field  directivity  and  the  angular  distri¬ 
bution  of  power  radiated  from  the  jet,  it  is  necessary  that  the  directivity  obtained  by 
summing  the  relative  directivities  of  each  segment  in  each  frequency  band  (k)  equal  the 
far  field  directivity  in  figures  31  to  34.  Hence: 


(HI-2) 
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where  Fj  (  0)  is  the  angular  distribution  function  for 
power  radiated  from  the  1^  segment  in 
frequency  band  (k).  F'^)  (9  )  is  the  far  field 
directivity  function  in  the  frequency  band  (k), 
and  the  summation, is  logarithmic. 

Note  that  this  equation  assumes  that  the  far  field  directivity  functions  Were  deter¬ 
mined  at  a  sufficient  distance  so  that  the  jet  approximated  a  point  source.  Hence,  no 
correction  was  made  for  the  effect  of  axial  displacement  of  the  segment  on  the  definition 
of  the  angle  3  . 
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Figure  34.  Generalized  for  field  directivity  in  the  2400  -  4800  and 
4800  -  9600  cps  octave  bands  for  various  turbojet  engines 
operating  at  military  power. 


The  cne-third  octave  band  sound  pressure  level  measurements  (Ref.  1)  given  in 
figure  24  were  converted  to  full  octave  band  sound  pressure  levels  as  shown  in  figure  35. 

The  three  highest  frequency  bands  have  been  extrapolated  toward  smaller  x//(0to  account 
for  the  power  radiated  in  the  1  ft.  distance  between  the  first  measurement  position  and  the 
nozzle.  Also,  some  smoothing  has  been  done  in  the  300  to  1200  cps  range  to  avoid  what 
is  believed  to  be  ground-reflected  noise.  The  onset  of  significant  reflection  appears  to 
occur  between  15  and  20  ft.  downstream  and  continues  from  there  on  to  approximately  40  to 
45  ft. ,  as  can  be  seen  In  figure  29.  These  reflections  can  be  qualitatively  evaluated  from 
inspection  of  figure  36  where  the  effect  of  the  ground  on  the  sound  pressure  maxima  is 
indicated.  The  directions  associated  with  the  sound  pressures  in  each  segment  are  derived 
from  the  longitudinal  correlation  studies  (Ref.  1)  which  are  summarized  in  figure  30. 


In  order  to  solve  equation  (1 1 1—2)  it  was  assumed  that  the  near  field  directivity  f|(0) 
would  vary  smoothly  with  x/xo.  It  was  further  assumed  that  for  a  sufficiently  low  frequency 
(below  37.5  cps  in  these  computations),,  the  whole  jet  could  be  characterized  by  one  va'ue 
of  X/A0'  **e* /  fhat  the  whole  jet  would  radiate  similarly.  With  this  latter  assumpt'on,  It 
was  possible  to  obtain  an  initial  series  of  points  on  the  fj(Q  )  versus  x/)^ curves,  utilizing 
unpublished  far  field  data  for  a  1/8  scale  J57  jet  engine  (Ref.  17)  for  an  octave  frequency 
band  centered  on  a  modified  Strouhal  number  (fd/ao)  of  .041 .  For  the  first  (37.5  -  75  cps) 
octave  band  of  the  full  scale  data,  the  jet  was  divided  into  two  segments  and  the  power 
calculated  for  each.  Applying  the  above  condition: 

z  z 

(.)  v->.  0) 


iolo<j  F  (6)  =  f  ;(e)j  -  L^pwq 


Vi) 


where  fj(0)  is  the  initial  series  of  points  at  10°  intervals 
and  f2(0  )  is  unknown. 


Solving  for  f2(0)  yields  a  second  series  of  points  on  the  f;(&)  versus  x /xocurves.  The 
procedure  for  the  second  octave  band  (75  -  150  cps)  was  similar,  except  that  the  jet  was 
divided  into  three  segments.  Now 

P-%)  =  lVJpw  L?  +  f,(e)]  -  |V  Pw 

where  f]  (  0)  and  fn($)  are  known  and  f^(0  )  is  unknown. 

This  iteration  continued  throughout  the  eighF  octave  bands  with  the  subscripts  maintaining 
their  identity  with  x/^,  as  the  distance  of  the  ■  enter  of  the  segment  from  the  nozzle  in 
wavelengths. 

During  the  c  ourse  of  the  iteration  it  often  became  apparent  that  a  value  of  fj(0  ) 
would  not  apply  to  other  octave  bands  or  would  not  produce  a  smoothly  varying  function. 
In  such  cases  other  values  were  estimated  and  tested  so  that  the  condition,  Eq.  (111-2) 
would  be  satisfied  to  within  one  decibel.  Figure  3.  shows  the  initial  stages  of  this  "trial 
and  error"  iterative  approach,  together  with  the  initial  ooints  for  two  values  of  O.  When 
the  first  iteration  was  completed  for  the  whole  range  of  0  ,  os  well  as  */\t ,  f;(B  )  was 
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Figure  36.  Illustration  of  effect  of  ground  plane  for  the  600  -  1200  cps 
octave  band. 


plotted  versus  0  with  x//4as  a  parameter.  These  curves  were  then  checked  to  insure  that 
the  integral  of.5f;($)  stnQ  between  0  and  IT was  unity,  as  it  should  be.  The  final  fj(Q) 
is  plotted  in  figure  38  as  a  function  of  Q  and  cross-plotted  in  figure  39  as  a  function  of 

x/Ao- 

It  would  be  desirable  to  confirm  fj(0)  through  experiment,  but  the  correlations  which 
are  required  are  one  order  of  magnitude  more  laborious  than  any  attempted  to  date.  How¬ 
ever,  lacking  direct  experimental  verification,  the  best  estimate  of  the  uniqueness  of  the 
fj(0)  given  In  figures  38  and  39  is  the  strong  tendency  toward  convergence  which  was 
noted  during  the  iteration  procedure.  Therefore,  with  the  exception  of  the  function  for 
x/\9  =  .375,  it  is  believed  that  these  functions  represent  a  relatively  unique  solution  to 
Eq.  (111-2),  and  are  thus  primarily  limit  d  in  accuracy  by  the  accuracy  of  the  generalized 
far  field  directivities  from  Ref.  16  In  figures  31  to  34  and  the  near  field  data  from  Ref.  1 
in  figure  35.  The  ability  of  these  functions  to  assist  in  the  prediction  of  acoustic  power 
will  be  discussed  in  the  next  subsection  and  the  prediction  of  near  field  noise  will  be 
considered  in  the  following  subsection. 
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Power  Distribution  Throughout  the  Jet 


In  the  derivation  of  the  angular  distribution  functions,  fj(0),  for  the  power  radiated 
at  various  values  of  x/^#  ,  there  was  no  requirement  that  the  total  power  computed  from 
the  near  field  equal  the  power  computed  from  far  field  data.  Hence,  a  comparison  of  these 
powers  affords  a  gross  test  of  the  model. 

For  the  first  iteration  for  f;(8),  the  powers  were  computed  using  Eq.  (lil-1)  for  each 
segment  based  on  the  primary  angle  of  radiation  along  the  boundary  (figure  30).  However, 
for  the  final  iteration,  the  angular  distribution  of  power  for  each  segment  was  considered, 
together  with  the  requirement  that  the  mean  square  pressure  at  any  point  on  the  10°  bound¬ 
ary  must  equal  the  sum  of  the  mean  square  pressures  computed  from  the  directed  powers 
flowing  through  the  boundary. 

The  results  of  the  final  power  computations  from  the  10°  boundary  measurements  are 
presented  in  figure  40,  The  agreement  between  the  power  computed  from  the  boundary  and 
that  computed  from  many  far  field  measurements  (Ref.  14)  is  almost  exact  below  300  cps. 
However*  above  300  cps  the  far  field  generalized  power  is  1.5  to  2.5  db  less  than  that 
computed  from  the  near  field.  There  are  many  possible  reasons  for  this  difference,  in  addi¬ 
tion  to  the  assumptions  inherent  in  the  basic  jet  radiation  concept.  These  reasons  include; 
the  ground  attenuation  which  lowers  the  far  field  sound  pressure  at  frequencies  above  300 
cps,  and  hence  the  computed  power;  atmospheric  attenuation  which  lowers  many  of  the  far 
field  sound  pressure  levels  at  frequencies  above  1200  cps;  and  the  assumption  of  uniform 
radiation  over  the  entire  hemisphere  which  is  common  to  all  far  field  powers  computed 
from  ground  measurement  of  horizontal  jets. 

With  these  additional  considerations,  it  is  clear  that  no  further  fundamental  conclusion 
can  be  reached  in  regard  to  the  correctness  of  either  measure  of  power.  To  obtain  a  more 
precise  conclusion  would  require  an  extremely  careful  experiment  in  free  spherical  space, 
under  Ideal  atmospheric  conditions,  correlating  both  near  and  far  fie  Id  measurements . 
Therefore,  It  is  felt  that  the  agreement  shown  in  figure  40  indicates  that  both  measures  of 
power  are  equally  valid,  and  that  the  gross  distribution  of  power  radiated  from  the  jet, 
predicted  by  the  previous  concepts,  is  valid. 

For  the  purpose  of  predicting  the  power  radiated  from  an  arbitrary  flow.  It  is  necessary 
to  examine  the  acoustic  power  generation  as  a  function  of  axial  distance  and  frequency  in 
relation  to  the  actual  flow  parameters.  The  power  per  unit  axial  length  along  the  jet 
(Wx)ob  where  OB  signifies  summation  over  each  octave  frequency  band  is  given  in 
figure  41 .  The  maximum  values  of  ( W  x  )qB  occur  at  values  of  x/re  less  than  1  for 
frequencies  above  2400  cps,  between  1  and  2  for  1200  and  2400  cps,  and  consistently 
further  downstream  for  each  successively  lower  frequency  band.  It  should  be  noted  that 
the  slope  of  (W  x)oB  *'fh  respect  to  axial  distance  for  the  highest  frequency  band  is 
approximately  (  x/r  )“^ ,  indicating  that  significant  contributions  to  the  total  radiated 
power  in  these  frequency  bands  results  from  ail  axial  locations  in  the  mixing  region  adjacent 
to  the  core.  In  contrast,  the  major  contribution  to  the  power  radiated  in  the  lowest  fre¬ 
quency  bands  crosses  the  10°  boundary  downstream  of  the  core. 


Throughout  the  mixing  region  adjacent  to  the  core,  the  typical  frequency  of  the 
turbulence  should  be  proportional  to  U/x.  Hence,  it  would  be  expected  that  the  spectra 
would  be  proportional  to  a  modified  Strouhal  number  with  the  usual  diameters  replaced  by 
x.  Figure  42  gives  the  relative  power  spectra!  density  ( Wfo/Wx  formalized  for  the  band¬ 
width,  as  a  function  of  a  modified  axial  Strouhal  number  for  the  data  of  figure  41  in  the 
region  of  (x/rg4.  13). 

In  the  region  far  downstream  of  the  core,  the  centerline  velocity  is  proportional  to 
(x/ rg  )“1  and  the  jet  width  parameter  b  is  proportional  to  x.  Therefore,  the  characteristic 
frequency  which  is  proportional  to  U/b  should  vary  as  (x/  re)-^  for  the  cold  jet,  rather 
than  relationship  (x/  rg)-^  found  in  the  core  region.  For  hot  jets,  the  variation  of  the 
speed  of  sound  in  this  region  would  be  approximately  proportional  to  (x/ re)“  ^  There¬ 
fore,  the  modified  Strouhal  number  should  contain  (x/re)^\  However,  the  data  given  in 
figure  43  do  not  agree  with  this  reasoning.  Rather,  it  is  seen  that  the  best  correlation  gives 
a  spectrum  which  is  independent  of  axial  distance. 


Apparent  angle  of  primary 
radiation  given  by  ray  acoustics 
in  this  region. 


Modified  axial  Strouhal  number  c~* 

U+oZ o 

Figure  42.  Normalized  power  spectral  density  along  10°  boundary  for  J57 
in  the  core  -egion. 
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The  relationships  between  power  and  axial  position  given  tn  the  previous  three  figures 
apply  directly  to  the  measurements  along  the  10°  boundary.  However,  if  these  relationships 
are  to  be  useful  in  the  prediction  of  power  generation  from  more  arbitrary  flows,  the  axial 
location  should  be  the  actual  tjrbulent  source  position  which  is  upstream  of  the  measurement 
position.  For  this  purpose  the  high  frequency  curves  of  figure  41  have  been  shifted  upstream 
in  accordance  with  the  average  axiol  distance  between  the  source  and  the  10°  boundary 
measurement  position,  as  illustrated  in  figure  27.  The  lower  frequency  curves  have  been 
shifted  by  a  greater  amount  in  accordance  with  their  higher  primary  radiation  angle  such 
that  the  maximum  value  of  Wx  occurs  at  the  presumed  core  tip.  The  results  in  figure  44 
shoula  be  a  better  approximation  of  the  true  source  locations. 


Figure  43.  Normalized  power  spectra  along  10°  boundary  for  J57 

in  downsrream  region,  x/r  >13. 

© 
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Figure  44.  Acoustic  power  per  unit  axial  length  in  octave  bands  shifted  to  probable 
source  locations 


The  overall  power  per  unit  oxial  length  computed  from  figures  41  and  44  is  presented 
in  figure  45.  It  is  interesting  to  note  that  Wx  is  essentially  constant  along  the  core  as  pre¬ 
dated  by  Ribner  (Ref.  18)  and  Dyer  (Ref.  19)  from  dimensional  considerations.  In  the 
downstream  region  Wx  is  falling  with  x/re  at  an  increasing  rate  approaching  (x/re)"3  at 
x/re  =  40.  When  the  calculatea  values  of  jet  breadth  and  centerline  are  utilized  in  this 
transition  region,  together  with  thepoU8b2  relationship,  the  curve  predicts  a  much  faster 
decrease  of  Wx  with  (x/re)  than  found  from  the  data.  Ribner's  analysis  (Ref.  18)  leads  to 
a  prediction  of  (x/re)~7  |n  the  final  downstream  region  of  the  jet. 

Figure  43  gives  a  good  indication  of  the  reason  why  the  noise  generated  in  the  transi¬ 
tion  region  does  not  approach  the  (*/te)-7  or  the  equivalent  DU2b2  relationships.  These 
relationships  assume  that  the  power  spectrum  at  any  station  throughout  the  jet  is  a  function 
of  U/b.  In  the  mixing  region  U  is  constant  and  b  is  directly  proportional  with  x  and  hence 
the  center  frequency  should  vary  inversely  with  (x/re).  Since  this  relationship  is  found 
true  in  figure  42,  Wx  is  constant  as  piadicted. 

In  the  downstream  region  the  same  consideration  for  the  frequency  spectra  proportional 
to  U/b  implies  that  the  frequency  spectre  should  vary  inversely  with  (x/re)2  since  U  is 
eventually  proportional  to  l/x  and  b  is  proportional  to  x.  However,  as  figure  43  shows, 
the  spectrum  is  apparently  independent  of  axial  position  for  13<x/re<50.  This  implies 


Figure  45.  Overall  acoustic  power  per  unit  axial  distonce  and  per  unit  nozzle 
area,  shifted  to  probable  source  location. 
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that  the  more  Intense  turbulence  generated  in  the  mixing  region  adjacent  to  the  core  does 
not  decay  as  fast  as  the  mean  jet  centerline  velocity,  which  is  implied  somewhat  in  figure 
16.  Therefore,  the  primary  generation  of  noise  in  this  transition  region  probably  results 
from  turbulence  developed  upstream  and  convected  downstream,  with  only  small  alteration 
through  the  transition  region,  rather  than  locally  generated  turbulence  which  follows  the 
local  flow  parameters. 

Figures  46  and  47  give  the  normalized  power  spectra  for  the  core  and  downstream 
regions,  respectively,  associated  with  the  "probable"  source  locations.  These  curves  are 
essentially  similar  to  figures  42  and  43,  except  for  the  slight  axial  shift;  however,  they  are 
believed  to  be  more  accurate  for  prediction  of  noise  generated  in  non-axisymmetric  flows. 

The  relative  proportions  of  total  power  radiated  from  the  regions  x/re  greater  and  less 
than  13  are  given  in  figure  48.  As  would  be  expected,  the  high  frequencies  are  primarily 
radiated  from  the  core  region  and  the  low  frequencies  from  further  downstream.  However,  it 
is  noted  that  significant  quantities  of  low  frequency  power  are  generated  in  the  core  region, 
and  high  frequencies  in  the  downstream  region.  Hence,  elimination  of  either  mixing  results 
in  a  decrease  of  approximately  6  db  in  the  power  radiated  in  one  of  Ithe  t^o  frequency  ranges. 
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Figure  48.  Relative  contribution  to  total  power  for  core  region  ^  6.5) 


and  downrtream  region  £  )6.5)as  a  function  of  dimensionless 

d 


frequency. 


Prediction  of  Near  Field  Noise  from  the  Radiation  Model 


The  preceding  analysis  of  the  radiation  of  noise  from  the  jet  has  shown  that  the  power 
computed  from  the  10°  boundary  is  approximately  equal  to  the  power  computed  from  far 
field  measurements.  Furthermore,  this  identity  was  utilized  to  determine  the  angular 
distribution  of  power  as  function  of  x/x  along  the  jet. 


A  more  severe  test  of  the  prediction  technique  is  its  ability  to  predict  near  field 
noise.  However,  before  discussing  the  result  of  this  comparison.  It  Is  necessary  to  consider 
the  influence  of  the  ground  plane  on  noise  measured  in  the  near  field  of  conventional  hori¬ 
zontal  jets. 


The  effects  of  reflections  of  sound  horn  a  ground  plane  of  variable  Impedance  were 
discussed  by  Ingard  (Refs.  27,  28)  and  portions  of  these  results  were  specialized  for  the  far 
field  for  both  random  and  sinusoidal  noise  by  Franken  (Ref.  29).  These  latter  results  were 
discussed  in  detail  by  Howes  (Ref.  30),  and  extended  by  Morgan  ef  al  (Ref.  10)  to  the  near 
field  by  Including  the  effect  of  the  difference  in  path  length  between  direct  ond  reflected 
rays  on  the  pressure  amplitude.  In  order  to  generalize  the  final  results  in  form  suitable  for 
application  to  arbitrary  locations  of  observation  and  source  points,  Franken's  approach  can 
be  extended  directly  to  the  near  field  os  follows. 
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Consider  a  source  height  (h)  and  receiver  height  (H)  with  a  separation  distance  (s)  as 
shown  in  figure  49.  The  path  length  difference  (A)  between  direct  and  reflected  sound  is 
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A  =  V  d2  +  4  Hh'  -  d  »»  2Hh 


(in  the  far  field)  (111-3) 


The  relative  path  length  (  ^  )  is  given  by 

Y-l/2 

i  i  “♦ini 

S  - 


d _ 

d  +  A 


,  4  Hh' 

1 


(in  the  far  field)  (1 1 1—4) 


where  the  time  delays  between  reflected  and  direct  wave  is  simplyA/ao, 
where  aQ  is  the  velocity  of  sound. 

The  sound  pressure  at  the  observation  point  resulting  from  a  single  phase  coherent 
source  is  the  sum  of  the  pressure  of  the  direct  wave  (pj)  and  the  pressure  of  the  reflected 
wave  (pr)/  and  the  resulting  mean  square  pressure  (p2)  at  the  observation  point  is: 


P2  =(Pd  +  Pr)2  =  Pd2  +  Pr2  +  2  PdPr 


(lli-5) 


Noting  that  pr  equals  §pd  when  the  reflecting  plane  has  infinite  impedance,  the  mean 
square  pressure  at  the  observation  point  is: 
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Figure  49.  Sketch  of  geometry  between  source 
and  observation  points. 
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p2  =pd2(i+82  +8&). 
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where  B  equals  twice  the  correlation  coefficient  between  pj  and  pr. 


As  shown  by  Franken  (Ref.  22),  B  is  given  by 

sin  2lUb'Y  -  sin  2TT fQT 

ITT  (fb  -  y 

where  fa  and  f^  are  the  upper  and  lower  frequency  limits  of  the  filter  bandwidth. 
For  white  noise  in  octave  bands,  this  expression  reduces  to 


B  = 


where 


pVT) 

+  f  b 


cos  2  iTf  X 

m  v 


(HI-7) 


f  =  a 
m 


When'fapproaches  zero,  B  approaches  two,  and  if  at  the  same  time  8  approaches 
unity,  the  mean  square  pressure  at  the  observation  point  is  four  times,  or  6  db  greater  than 
the  direct  mean  square  pressure,  as  would  be  expected.  On  the  other  hand,  when  f^ 
becomes  large,  B  approaches  zero,  so  that  the  observed  mean  square  pressure  Is  twice,  or 
3  db,  greater  than  the  mean  square  pressure  for  the  direct  path.  The  relationships  for  the 
near  field  case  between and  s  are  given  in  figure  50,  and  the  relationship  between  B 
and  fmX'is  given  In  figure  51. 

An  experimental  verification  of  these  relationships  is  summarized  in  figure  52.  The 
data  were  obtained  at  25  feet  distance  from  a  3  inch  diameter  hot  jet  whose  centerline  was 
9  inches  above  a  concrete  surfuce.  The  parameter  Hh  was  varied  by  varying  the  microphone 
height  between  9  and  48  inches,  and  fm  was  varied  in  accordance  with  the  octave  bands 
measured. 

It  should  be  noted  that  these  relationships  only  apply  when  the  ground  surface  is  hard 
and  non-absorbing.  When  this  restriction  is  not  satisfied,  the  reflected  wave  is  partially 
absorbed  by  the  ground,  and  at  great  distances  the  direct  wave  is  also  partially  absorbed. 
(Refs.  27,  28)  This  ground  absorption  is  clenriy  illustrated  in  figure  53,  which  gives  the 
average  sound  pressure  level  measured  in  octave  bands  at  250  ft.  from  several  J57  engines 
operating  at  military  power  at  different  locations.  The  dip  in  response  predicted  by  figure 
51  occurs  in  the  1200  -  2400  cps  frequency  band,  and  is  present  in  the  measurements  over 
concrete.  However,  the  large  variation  in  overage  sound  pressure  level  evident  at  fre¬ 
quencies  above  150  cps  for  the  measurements  over  absorbing  ground  surfaces  are  not 
predictable  by  simple  reflection  theory. 
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20  50 

s  feet 


Figure  50.  The  relations  between  relative  path  length  &  t  and  time  delay  ^ 
and  the  distances  between  observation  and  source  points  as  a 
function  of  the  product  of  source  and  receiver  heights  (Hh)  for 
computation  of  the  sum  of  direct  and  reflected  noise. 
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'p‘=?/a-s‘*s?) 

where  ^2‘  is  the  mean  square  pressure  for  the  direct  path 


is  the  mean  square  pressure  resulting  from  both  direct 
and  reflected  paths 


Center  freqo*  <*«'  <imes  delay  <wr> 


Figure  51.  The  parameter  B  (twice  hie  correlation  coefficient  between 
direct  and  reflected  wave)  as  a  function  of^/^for  octove 
band  white  noise  and  infinite  impedance  reflecting  surface 
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Figure  52.  Comparison  of  data  obfained  with  1/8  scale  jet  at  25  ft.  radius,  9  in. 

above  ground  and  for  various  microphone  heights  between  9  and  48  in. 
with  computed  reflection  effect. 


Extreme  caution  must  be  exercised  In  applying  these  results  to  near  field  jet  noise, 
even  when  the  reflecting  surfaces  are  hard,  because  of  the  assumption  of  a  phase  coherers? 
source.  In  the  near  field  the  extended  jet  source  becomes  large  relative  to  the  distance 
between  observation  point  and  receiver.  Hence,  it  is  necessary  to  divide  the  jet  into 
several  segments  and  compute  the  noise  at  the  observation  point  resulting  from  each  segment. 
Since  the  correlation  of  the  noise  generated  by  adjacent  segments  can  be  assumed  to  be'zero, 
the  contributions  from  all  segments  can  be  summed  on  an  energy  basis,  where  the  effect  of 
ground  reflections  is  accounted  for  separately  in  the  computation  of  the  noise  from  each 
segment.  However,  when  the  observation  point  approaches  close  to  the  jet,  the  angle 
between  the  direct  and  reflected  path  (  oc)  becomes  large,  and  there  is  no  longer  any 
guarantee  of  source  coherence,  since  the  noise  generated  in  the  mixing  region  around  the 
jet  axis  at  any  station  is  only  in  phase  over  a  limited  angular  region.  Only  very  limited 
data  exist  in  Ref.  1  relative  to  this  point,  and  it  must  remain  an  uncertain  factor  at  this  time. 

Predictions  of  near  field  noise  were  first  made  for  six  positions  extending  from  40  feet 
forward  to  60  feet  aft  of  the  nozzle  along  a  line  parallel  to,  and  20  feet  from,  the  jet  axis. 
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Legend:  \\  " 

■  North  Island  NAS,  Mar.  61 
O  Wright-Patterson  AF6,  Nov.  I9604  ~© 

©  Miramar  NAS,  February  1959** 

A  George  AFB,  September  1961**  V'v 

A  Miramar  NAS,  May  I960** 

□  Truax  AFB,  January  1959**  \ 

*Data  measured  over  concrete  surface 

**Dato  measured  over  surface  partly  covered  with  weeds  and  grass 


pressure 


This  choice  enabled  a  comparison  with  data  from  independent  noise  r.  }efs.  1,  3, 

20  to  26)  which  were  made  on  J57  engines  in  conjunction  with  the  ev  .  ‘’~.n  of  jet  engine 
ground  runup  noise  suppressors. 

For  the  purpose  of  prediction,  the  |et  was  divided  into  five  segments  with  average 
distances  from  the  nozzle  of  2.1,  4.2,  8.4,  16.8  and  33.6  feet,  respectively,  and  the 
power  per  segment  was  calculated  for  each  octave  band. 

Each  point  in  the  near  field  was  characterized  by  its  distance  from  the  center  of  the 
segments  and  angle  relative  to  the  axis  of  the  jet.  Thus,  for  the  k^1  octave  the  sound 
pressure  level  is  5 

spl;(k>  =  L X  [PWL|(k)  +  f;  (8)  -  10  log  (4Tf  dj2)+  10  log  (1  +§2  +  $  B)1 
i=1  (111-8) 

where  PWLj^)  is  the  power  of  the  i^b  slice  for  the  k*b  octave,  f.(8)  Is  the 
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Axial  distance  relative  to  the  nozzle,  in  feet 

Figure  54.  Predicted  values  of  sound  pressure  level  in  four  octave  bonds  along  a  line 
20  feet  from  and  parallel  to  the  axis  of  the  jet,  considering  both  simple 
hemispheric  rodiofion  and  correlated  ground  reflections,  compared  to  the 
range  and  median  values  of  nine  measurements  on  o  J57  engine  at  military 
power . 


Sound  pressure  level  in  db  re  .0002  dyne/cm 


directivity  for  the  appropriate  *4, of  the  sl  ice,  and  dj  is  the  distance  from  the  center  of 
the  i^  slice. 

The  results  are  shown  in  figures  54  and  55,  for  both  hemispheric  radiation,  where  the 
reflected  ray  is  assumed  uncorrelated  with  the  direct  ray  (B  =  0),  and  for  correlated  reflected 
and  direct  waves.  Only  one  value  is  given  in  figure  55  because  the  value  of  fm'C  is  suffi¬ 
ciently  high  at  these  frequencies  to  insure  B  =  0  and,  therefore,  simple  hemispherical  radia¬ 
tion  exists.  On  each  graph  the  average  error  between  prediction  and  median  of  the  nine 
measurements  is  given.  In  general,  the  predictions  are  within  the  range  of  the  measured 
data,  and  the  fit  appears  most  encourcging,  particularly  between  600  and  4800  cps.  It  is 
interesting  to  note  that  in  the  frequency  region  above  600  cps  where  the  model  predicted 
far  field  power  2  db  greater  than  computed  from  the  far  field  measurements,  the  20  foot 
line  predictions  average  less  than  1  db  higher  than  the  .  -dian.  Inclusion  of  the  effect  of 
correlation  between  direct  and  reflected  waves  (B  -  0)  in  the  low  frequency  range  appears 
particularly  appropriate  in  the  150  -  300  cps  band  where  the  correlation  gives  cancellation. 


Axial  distance  relative  to  the  nozzle,  in  feet 

Figure  55.  Predicted  values  of  sound  pressure  level  in  four  octave  bands  along  a  line 
20  feet  from  and  parallel  to  the  axis  of  the  jet,  compared  to  the  range 
and  median  values  of  nine  measurements  on  a  J57  engine  at  military  power. 
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However,  this  correction  increases  the  average  error  in  the  37  -  75  cps  band  from  1.5  to 
4  db. 


Contours  of  equal  sound  pressure  level  were  computed  for  three  octave  bands,  75  - 
150,  300  -  600,  and  1200  -  2400  cps,  as  shown  in  figure  56.  For  these  computations, 
ground  reflections  were  assumed  to  be  incoherent  (B=  0),  giving  a  3  db  increase  over  the 
direct  wave.  It  is  clear  that  this  example  of  idealized  contours  behaves  in  the  same  manner 
as  the  measured  contours  in  figures  21  and  22,  with  the  low  frequencies  radiating  farther  aft 
than  the  high  frequencies  and  with  the  steep  gradient  of  sound  pressure  level  forward  from 
the  nozzle. 

The  method  also  enables  computation  of  the  apparent  jet  source  location  of  the  sound 
reaching  any  point  in  the  near  field.  This 

splj.  +  10  log 

where  spl|  is  the  sound  pressure  contributed  by  each  segment  of  the  jet, 
spl).  is  the  total  sound  pressure  level  at  the  point,  and  x;  is  the 
distance  downstream  for  the  center  of  each  slice. 

These  computations  were  performed  for  the  two  octaves  between  300  and  1200  cps  for 
the  observation  points  located  at  12  and  32  feet  downstream  and  22.5  feet  from  the  jet  axis. 
These  are  the  points  for  which  correlation  data  are  available  (Ref.  1),  as  illustrated  in 
figure  29.  For  a  center  frequency  of  1000,  the  angles  derived  from  data  are  102°  and  124° 
for  the  two  positions,  respectively,  as  compared  to  102°  and  137°  calculated  from  the  model 
using  the  octave  band  centered  at  850  cps.  The  angles  from  data  for  a  center  frequency  of 
400  cps  are  103°  and  130°.  as  compared  to  95°  and  132°  calculated  for  the  octave  band 
centered  at  425  cps.  Thus,  it  can  be  seen  from  the  model  that  the  axial  location  of  the 
apparent  source  depends  upon  the  angle  between  the  observer  and  the  jet  axis.  This  explains 
why  the  correlations  performed  in  the  field  often  have  indicated  different  apparent  source 
locations  which  were  often  ascribed  to  inaccuracy. 

Effect  of  Jet  Exit  Temperature  and  Mach  Number  on  Directivity 

The  relationship  between  the  directional  properties  of  jet  noise  and  the  jet  flow  para¬ 
meters  have  been  a  subject  of  considerable  speculation  for  many  years.  The  gross  effect  of 
temperature  on  the  far  field  directivity  is  shown  in  figure  57  for  Mach  1  axisymmetric  jets. 

It  is  clear  from  the  figure  that  the  increase  of  temperature  shifts  the  angle  of  maximum  radia¬ 
tion  forward. 

This  shift  in  the  maximum  angle  of  noise  from  jets  at  constant  Mach  number  results 
from  refraction  in  the  jet,  as  first  suggested  by  Hubbard  (Ref.  11).  A  step  toward  improved 
understanding  of  this  refraction  is  the  theoretical  analysis  in  Appendix  C  for  propagation 
from  a  source  on  the  axis  of  a  zero  velocity  jet  which  has  a  cylindrically  symmetrical  temper¬ 
ature  gradient.  These  results  are  not  directly  applicable  to  the  mixing  region  adjacent  to 
the  core  of  a  real  jet  which  has  important  velocity  gradients  and  has  sources  located  approxi¬ 
mately  at  the  maximum  portion  of  these  gradients,  as  shown  in  figures  1  to  3.  However,  the 


requires  solving  for  x  in 

•  5 


r 


X  =  log  spl;  +  1C  log  x; 

i  =1 


(HI-9) 
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assumptions  In  the  theory  are  approximated  In  the  downstream  region  of  a  jet,  say  beyond 
20  diameters.  Further,  these  results  are  useful  in  demonstrating  the  potential  magnitude  of 
temperature  refraction  in  affecting  the  far  field  directivity. 

Because  of  the  time-consuming  difficulties  involved  in  removing  the  restrictions 
inherent  in  the  Appendix  C  theory,  a  series  of  ray  diagrams  were  constructed  for  sources 
located  at  approximately  the  maximum  velocity  gradient  and  for  radiation  in  the  plane 
containing  the  jet  axis.  As  noted  in  conjunction  with  the  discussion  of  figure  27,  the 
temperature  and  velocity  gradients  were  assumed  coincident  and  were  taken  directly  from 
the  Laurence  data  of  Ref.  31 . 

Figures  58  and  59  give  two  additional  examples  of  ray  diagrams  for  Mach  1  jets  at 
exit  static  temperatures  of  60°F  and  2700°F  respectively.  The  results  of  these  two  figures, 
figure  27,  and  several  other  similar  constructions,  are  summarized  in  figure  60.  It  is  inter¬ 
esting  to  note  that  ray  refraction  shifts  the  angle  of  maximum  radiation  forward,  away  from 
the  jet  boundary,  both  for  an  increase  in  flow  Mach  number  at  constant  temperature,  and  for 
an  increase  in  temperature  at  constant  flow  velocity.  Recognizing,  in  accordance  with 


Figure  57.  Comparison  of  the  directivity  of  axisymmetric  jets  with  exit 

Mach  numbers  of  unity  for  three  values  of  exit  static  temperature. 
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source  located  at  the  center  of  the  jet  velocity  and  temperature 
'sec.  and  exit  temperature  of  60°F. 


Ray  diagram  for  radiation  from  simple  source  located  at  the  center  of  the  {et  velocity  a 
temperature  gradient,  for  exit  velocity  of  2750  ft/sec.  and  exit  temperature  of  2700°  F 


from  upsfream  fef  axis,  in  degrees 


These  results  offer  an  opportunity  to  compare  the  various  terms  in  the  new  modified 
wave  equation  for  jet  noise  propagation  derived  in  Appendix  B,  as  given  on  the  left  hand 
side  of  equation  (B-92).  For  propagation  of  sound  in  a  jet  with  a  two  dimensional  gradient 
with  no  sources  present,  this  equation  becomes 


py 

Dt1 


lie 
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(111-10) 


where  'P  is  the  acoustic  pressure 

and  -JpLis  the  total  differential  operator 

D* 

If  the  right  hand  gradient  term  in  equation  (111-10)  is  much  smaller  than  either  of  the  other 
two  bracketed  terms,  then  it  may  be  neglected  and  equation  (Ili-iO)  reduces  to  the  standard 
second  order  differential  equation  for  propagation  of  sound  in  moving  medium: 


1  _ 

Dt" 


CL1 


=  O 


(111-11) 


As  shown  in  Appendix  B,  the  gradient  term  in  equation  (111-10)  is  of  the  same  order  as 
the  first  term  inside  the  brackets  for  a  typical  Fourier  component  of  sound  when 
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For  the  velocity  gradient  given  in  figure  12,  the  maximum  value  of 

OC  U  ___  •  4?  Ue 

CL  ^  ^  and  occur£  at  y  =  b. 

Assuming  a  typical  value  of  M  for  the  jet  as  .5  Mg,  equation  (ill-12)  becomes 


(111-12) 


J  _  .(o  Ue  __  ,0>  Q.&  Me  _  1.2.  KA*- 

'  1tfk'(l+..5Mey  -rrfb(l'V-.'SMe)'5  «"-'3> 

where  ke  is  the  wave  number^~^~t  phased  on  the  speed  of  sound  ae 
at  the  jet  exit  temperature. 

When  Mq  =  1 ,  the  solution  of  equation  (111-13)  gives  keb  =  .36,  or  in  terms  of  a  wave 
number  (k0)defined  for  the  ambient  speed  of  sound; 
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(111-14) 


k0b  -  .36  °e 


For  the  axisymmetric  Mach  1  Jet,  b  =  re  at  the  tip  of  the  core,  b/x  Is  essentially  constant 
between  the  nozzle  and  the  core,  and  the  length  of  the  core  (x*)  is,  from  figure  14,  approx¬ 
imately  13  re.  Combining  these  relationships,  equation  (111-14)  may  be  restated  in  terms  of 
*A>/  giving 


—  =  13  (k0b)  =  .75  fi 
Ao  2  TT  co 


(111-15) 


When  ae  =  1850  ft/sec. ,  the  value  of  x/^  is  1,25  for  approximate  equality  between 
the  gradient  term  and  either  one  of  the  other  terms  in  equation  (111-10).  This  is  consistent 
with  figure  30,  for  at  small  values  of  x/^Q  (of  the  order  of  1 .25  and  lower),  where  the 
gradient  term  has  importance,  the  angle  of  maximum  radiation  from  the  jet  lies  close  to  the 
jet  boundary.  However,  when  the  gradient  term  becomes  much  smaller  than  either  of  the 
other  terms  in  equation  (I li— 10),  (x/^j  ^  5  implies  the  gradient  term  is  less  than  25%  of 
either  of  the  other  two  forms)  the  angle  of  primary  radiation  becomes  a  constant  which 
approximates  the  angle  given  by  simple  ray  considerations.  This  result  is  consistent  with 
equation  (ill— 1 1 ),  since  the  ray  solution  is  an  exact  solution  for  equation  (111-11)  when  the 
bending  of  the  ray  is  small  in  the  distance  of  one  wavelength. 


This  simple  evaluation  of  equation  (111-10)  can  also  demonstrate  the  relationship  of 
the  gradient  term  of  equation  (111-10)  with  the  generalized  power  spectrum  given  in  figure 
46.  For  a  Mach  1  jet,  Ue  =  ae,  and  the  modified  Strouhal  number  absissa  in  figure  46 
becomes  equal  to  x/^,.  The  spectrum  in  figure  46  has  a  low  frequency  asymptote  where  the 
acoustic  power/cps  increases  with  increasing  x^0,  which  is  valid  below  xA0  of  approxi¬ 
mately  1.0.  Above  the  low  frequencies  is  a  transition  region,  approximately  1 .0Cx^o^2.5, 
and  a  high  frequency  asymptote,  approximately  x^D  >2.5.  It  is  pleasantly  surprising  that 
these  three  regions  agree  so  well  with  the  simple  order  of  magnitude  analysis  of  the  im¬ 
portance  of  the  gradient  term  in  equation  (111-10).  The  clear  implication  is  that  the  full 
equation  (111-10)  controls  radiation  of  sound  in  the  jet  mixing  region  adjacent  to  the  core 
for  nondimensional  frequencies  beiow  and  including  the  spectral  maximum,  and  that 
equation  (111-11)  is  applicable  only  at  nondimensional  frequencies  above  the  spectral 
maximum. 


Therefore,  the  theoretical  solution  of  noise  radiation  and  generation  in  jet  flows  must 
rest  on  the  solution  of  equation  (B-92),  of  which  equation  (111-10)  was  a  simplified  version. 
Until  these  solutions  are  available,  it  is  not  possible  to  extend  the  near  field  prediction 
model  developed  in  this  section  for  the  1000°F  Mach  1  jet  to  other  temperatures  and  Mach 
numbers,  except  through  semi-empirical  derivations  similar  to  that  utilized  for  the  near 
field  model  described  for  the  1000°F  jet.  These  derivations  should  be  based  on  a  careful 
set  of  new  model  measurements  made  in  free  space  without  a  ground  plane,  where  exit 
Mach  number  and  temperature  are  varied  independently.  In  the  absence  of  data,  the 
first  approximation  to  an  extension  of  the  model  in  figure  38  to  include  temperature  effects 
would  be  to  shift  the  maximum  angles  and  curves  by  the  ratio  of  ray  angle  for  the  new 
temperature  from  figure  60  to  the  ray  angle  for  1000°F. 
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In  order  to  make  preliminary  prediction  of  the  near  field  noise  from  jets  with  exit 
velocities  lower  than  1850  ft/sec.  and  1000°F,  it  is  best  to  extrapolate  from  empirical 
data.  For  this  purpose  the  most  comprehensive  set  of  near  field  measurements  for  heated 
jets  was  made  by  Wolfe  (Ref.  2)  over  a  velocity  range  of  1000  -  1800  ft/sec.  Figure  61, 
taken  directly  from  Wolfe's  figure  31  (Ref.  2),  gives  the  value  of  n  in  equation  (111-16) 
below: 

spl‘  =  splref +  10  n  log  -ggQ  ■  (111-16) 

where  splref  refers  to  the  values  obtained  for  the  1000°F  jet 

at  a  point,  either  by  the  prediction  model  or  directly 
from  the  data  of  Ref.  1,  and  spl'  is  the  spl  estimated 
for  the  jet  of  differing  velocity. 

Unfortunately,  this  method  does  not  separate  the  effect  of  flow  Mach  number  from 
flow  temperature  and  for  accurate  prediction  of  the  spacial  distribution  of  noise  from  a  jet 
flow,  these  two  parameters  must  be  considered  separately.  However,  until  the  series  of 
experiments  mentioned  above  are  accomplished,  these  data  are  probably  the  most  reliable 
and  useful  for  this  purpose. 


LIST  OF  SYMBOLS  USED  IN  SECTION  III 


Roman 


A 

a 

B 

b 

d 

F(9) 

m 

f 

H 

h 

k 

z 

n 

P 

PWL 

s 

spl 

r 

Sx 

U 

(Wx)OB 

x 

y 

M 


a 

b 

d 


Radiating  area 
Speed  of  sound 

Twice  correlation  coefficient 
Gradient  width  parameter 
Exit  diameter 

Far  field  directivity  function 

Near  field  angular  power  distribution  function 

Frequency 

Receiver  height 

Source  height 

=  2  TT  f/a 

Prandtl  mixing  length 
Power  of  U 
Pressure 

Power  level  re  10"^  watts 

Source  separation  distance 

Sound  pressure  level  re  .0002  dyne/cm2 

Radius 

Strouhal  number,  modified 
Velocity 

Power  per  unit  axial  length  summed  over  all  octave  bands 
Axial  distance  from  exit  plane 
Latercl  distance  from  centerline 
Mach  number 


Greek 


Relative  path  length 

Angle  measured  from  upstream 

Acoustic  wavelength 

Logarithmic  sum 

Spacial  correlation  coefficient 

Delay  time  and  reflected  signal  re  direct  signal 

Angle  of  sound  pressure  maximum  re  normal  to  radiating  area,  A 

Subscripts 

Lower  octave  band  cutoff 
Upper  octave  band  cutoff 
Direct 

(continued) 
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List  of  Symbols  Used  in  Section  ill  (continued) 


Subscripts 


e  Exit 

k  Octave  band  number 

m  Center  of  frequency  band 

o  Characteristic  of  ambient  air 

r  Reflected 
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SECTION  IV 


CONTROL  OF  JET  NOISE  THROUGH  DESIGN 


O  re  of  the  primary  motivations  for  jet  noise  research  is  the  history  of  sonic  fatigue 
w 'd  equipment  malfunctions  of  current  flight  vehicles.  These  noise-induced  failures  are  not 
only  costi/  in  an  economic  sense,  but  threaten  a  vehicle's  overall  reliability  and  effective- 
nes.  as  a  weapors  system.  Although  experience  has  shown  that  satisfactory  reliability  can 
be  achieved  after  lengthy  testing,  redesign,  and  upgrading  of  structure  and  components, 
this  process  is  rcstly  in  time,  money,  and  often  performance.  If  this  history  is  not  to  be 
repeated  in  the  design  of  future  vehicles,  the  noise  exposure  of  flight  vehicles  must  be 
minimized  in  the  design  process. 

The  various  methods  or  devices  which  have  been  or  can  be  proposed  for  altering 
the  jet  noise  field  involve  variation  of  one  or  more  of  four  major  factors: 

Engine  location 

Total  acoustic  power  generated 

Frequency  spectra 

Radiation  characteristics  (directional  properties, 
reflecting  or  absorbing  surfaces,  and  spacial 
distribution  of  noise  sources) 

The  first  of  these  factors  can  be  controlled  by  the  vehicle  designer,  and  the 
remaining  are  controlled  primarily  by  the  engine  designer,  often  in  conjunction  with  the 
vehicle  designer. 

This  section  discusses  these  basic  factors  and  examines  the  modifications  of  these 
factors  which  can  be  achieved  through  design.  Much  of  this  discussion  will  be  given  directly 
in  terms  of  a  1 0, 000— lb.  thrust  engine  to  facilitate  comparison  with  current  experience. 
Furthermore,  the  discussion  is  directed  primarily  toward  jets  from  convergent  nozzles,  and  does 
not  include  supersonic  rockets. 

Engine  Location 

It  U  obundently  clear  from  the  vast  literature  on  jet  noise  and  the  discussions  of  the 
previous  section  that  the  noise  radiated  by  a  jet  to  loca.'ions  forward  of  the  nozzle  plane  is 
an  order  of  magnitude  lower  than  that  radiated  to  locations  aft  of  the  nozzle  plane.  Therefore, 
it  behooves  the  designer  to  accept  this  bonus,  wherever  possible  in  the  initiation  of  new  design 
concepts. 

To  illustrate  this  principle,  the  noise  exposure  of  a  twin  engine  jet  flight  vehicle 
was  estimated  for  wing  mounted  and  tail  mounted  jet  engines.  The  contours  of  overall  spl, 
given  Ir,  figure  21  for  a  1 0, 000- 'b. thrust  engine,  were  superimposed  over  the  vehicle  plan 
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and  fhe  approximate  percentage  of  total  vehicle  surface  exposed  to  6  db  increments  of 
noise  was  computed.  The  results,  given  in  figure  62,  show  that  fhe  overall  noise  exposure 
of  the  vehicle  is  reduced  over  12  db  by  moving  the  engines  aft  and  that  areas  of  extreme 
noise  exposure  are  eliminated  almost  entirely.  Although  these  results  are  only  approximate 
since  ground  reflections  were  not  considered,  they  demonstrate  conclusively  that  a  large 
reduction  in  vehicle  noise  exposure  and  resulting  vibration  can  be  achieved  by  placing  the 
engines  at  the  tail  of  the  vehicle.  Hence,  configurations  incorporating  tail  mounted  engines 
will  have,  in  general,  minimal  sonic  fatigue,  and  equipment  malfunction  from  noise  induced 
vibration. 


Fundamental  Jet  Noise  Parameters 

The  three  remaining  basic  factors  which  can  be  utilized  by  the  designer  to  control 
vehicle  exposure  to  jet  noise  are  the  total  acoustic  power,  frequency  spectra,  and  directivity 
of  the  jet  noise  source.  These  factors  are  determined  by  the  jet  flow  parameters,  and  are, 
except  for  variation  in  nozzle  geometry,  almost  entirely  established  by  the  engine  designer 
at  the  initiation  of  each  new  engine  design.  Clearly,  if  minimum  noise  is  to  be  achieved, 
low  noise  must  be  a  requirement  in  advanced  engine  design  to  insure  that  it  is  considered 
together  with  all  other  engine  objectives. 

The  total  acoustic  power  for  jet  flows  from  convergent  circular  nozzles  has  been 
shown  by  many  investigators,  including  (Ref.  1)  to  be  proportional  to  the  eighth  power  of 
the  nozzle  velocity  and  to  the  nozzle  area. 

An  empirical  equation  for  predicting  the  total  acoustic  power  level  (PWL)  is:* 

PWl  =  146  +  20  log  d1  +80  log  ,  in  db  re  10"13  watts  (IV-1) 

1000 

where  Ue  is  the  actual  nozzle  exit  velocity  in  ft/sec. 

and  d]  is  the  nozzle  diameter  (de)  in  fee*  when  the  nozzle 
pressure  ratio  (P.R.)  is  less  thon  1.89,  and  is  otherwise 
given  by 

d,  =de  |l  +  1.71  (.53  P.R.  -  1)j  l/2;  1 .89  ^  PR  <  3 

Note  that  d|,  as  defined  for  nozzles  operating  at  pressure  ratios  greater  than  critical,  is 
simply  the  diameter  of  a  nozzle  operating  at  1.89  P.R.  with  the  same  total  temperature  and 
thrust  as  the  actual  nozzle . 

The  overriding  effect  of  exit  velocity  on  the  acoustic  power  generated  by  the  jet  is 

•Note  that  the  first  approximation  of  the  total  power  radiated  in  flight  where  the  external 
velocity  U<j  &  0,  is  obtained  by  replacing  Ue  by  (Ue  -  UQ)  and  adding  ten  times  the  logarithm 
of  the  ratio  of  the  core  length  for  UQ  *  0  to  the  core  length  for  U0  =  0  from  figure  18  to 
account  for  the  increased  length  of  the  jet  noise  source. 
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summarized  for  a  10,000  lb.  thrust  engine  in  figure  63  for  various  values  of  P.  R.  Ail  exit 
velocities  along  the  1.89  P.  R.  line  are  sonic  and,  thus,  along  this  line,  lower  velocities 
imply  lower  exit  total  temperature  (T^e)  and  higher  weight  flows  (we)  through  the  engine. 
For  a  non-bypass  engine  with  a  fixed  engine  P.R.  and  exit  area,  the  lower  temperature 
generally  results  in  lower  thermal  efficiency.  Further,  the  increase  in  weight  flow,  either 
from  lower  temperatures  or  reduced  nozzle  P.R.,  results  in  increased  engine  size  and  high 
speed  drag  and  weight.  Therefore,  for  the  conventional  non-bypass  jet  it  is  necessary  to 
maintain  maximum  practical  temperatures,  notwithstanding  the  large  noise  reductions  attain¬ 
able  with  lower  exit  velocities.  As  will  be  seen  later,  the  practical  methods  by  which  the 
designer  can  obtain  the  acoustic  benefits  of  reduced  velocity  include  the  bypass  or  turbofan 
engine  and  to  a  lesser  extent,  the  various  types  of  mixing  nozzles  with  or  without  ejectors. 
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Figure  63.  Relative  total  acoustic  power  level  for  a  10,000-lb.  thrust 
engine  as  a  function  of  exit  velocity.  Data  from  Lee,  et  al 


The  normalized  frequency  spectrum  of  the  total  acoustic  power  radiated  from  axi- 
symmetric  jets  is  given  in  figure  64.  For  jets  from  convergent  nozzles  at  a  P.R.  of  1.89,  the 
exit  velocity  ( Ue)  equals  the  critical  sonic  velocity  (aj  at  the  exit  plane  or  throat  of  the 
nozzle.  Hence,  for  standard  ambient  conditions  the  frequency  spectrum  is  dependent  only 
on  the  diameter  of  the  jet,  with  smaller  jets  characterized  by  higher  frequencies.  Thus,  the 
entire  frequency  spectrum  of  a  jet  engine  could  be  shifted  to  much  higher  frequencies  if  the 
single  nozzles  were  replaced  by  a  number  of  small  nozzles  with  equal  total  area.  In  fact,  as 
first  noted  by  Tyler  (Refs.  3,  4)  and  proved  in  his  ground  runup  muffler  (Ref.  5),  if  the  nozzles 
were  sufficiently  small,  the  characteristic  center  frequencies  of  the  noise  would  be  in  the  ultra¬ 
sonic  frequency  range,  which  is  both  inaudible  and  rapidly  attenuated  by  atmospheric  absorp¬ 
tion. 

The  power  spectra  for  a  10,000  lb.  thrust  engine  at  a  P.R.  of  2.2,  overall  acoustic 
power  level  of  174  db,  and  a  variety  of  nozzle  diameters,  are  given  in  figure  65.  Since  the 
number  of  nozzles  required  for  constant  thrust  is  inversely  proportional  to  the  square  of  the 
nozzle  diameter,  it  is  clear  that  rather  large  numbers  of  nozzles  are  required  to  obtain  an 
appreciable  shift  In  center  frequency.  For  example,  32  four  inch  nozzles  are  required  to 
shift  the  frequency  by  a  factor  of  approximately  5.5  and  126  two  inch  nozzles  are  required 
to  give  a  frequency  shift  by  a  factor  of  approximately  11.  The  potential  reductions  in  the 
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=  1850  ft/sec. 
as  noted 


a  of  non-interfering  axisymmetric  jets  having  a  thrust  of 
nozzle  P.  R.  2.2,  exit  total  temperature  1720°R,  together 
s  of  lower  temperature  and  exit  velocity. 


lower  frequency  noise  ore  impressive,  approximately  21  and  30  db  for  the  four  inch  and  two 
inch  nozzles,  respectively.  However,  these  low  frequency  reductions  are  predicated  on  the 
assumption  that  the  individual  jets  do  not  interfere  with  each  other  and  re-combine  in  a  large 
jet  of  lower  velocity.  This  assumption  requires  large  center-to-center  distance,  and  hence, 
a  rather  large  envelope  area  for  the  entire  configuration.  If  the  jets  are  spaced  closer  than 
approximately  three  diameters,  interference  results  which  increases  the  low  frequency  radia¬ 
tion  above  that  given  in  figure  65,  and  decreases  the  high  frequency  noise  rar'!ation  relative 
to  figure  65. 

In  contrast  to  the  large  potential  improvements  obtainable  at  low  frequencies  from 
multiple  non-interfering  jets,  there  is  an  increase  in  noise  radiated  at  high  frequencies  over 
that  radiated  by  the  single  jet  of  equal  thrust.  This  increase  is  of  the  order  of  7  db  and  10  db 
for  the  four  inch  and  two  inch  jets  respectively.  However,  in  this  frequency  region  interfer¬ 
ence  between  the  jet  flows  is  beneficial,  reducing  the  magnitude  of  the  increase. 

Figure  65  also  shows  the  estimated  spectra  for  two  single  10,000  lb.  thrust  jets  oper¬ 
ating  at  the  same  P.  R.,  but  at  lower  temperature,  and  thus,  at  lower  velocity.  Note  that  the 
center  frequency  of  these  lower  velocity  jets  is  unaltered  from  the  original,  since  Ue  =  a*  in 
all  three  cases.  Therefore,  the  reduction  in  total  acoustic  power  estimated  in  figure  63  applies 
to  all  frequencies. 

The  near  field  directional  properties  of  the  ax i symmetric  jet  ct  1000°F  have  been 
discussed  in  detail  in  the  previous  section,  and  the  effect  of  temperature  and  velocity  on 
these  directional  properties  was  discussed  briefly  in  Section  III,  as  well  as  in  Appendix  C. 

As  was  noted.  It  is  still  impossible  to  give  a  quantitative  generalization  of  the  effects  of 
velocity  and  temperature,  singly  or  combined,  on  the  noise  at  an  arbitrary  point  in  the  near 
field  of  a  jet  issuing  from  a  convergent  circular  nozzle.  Therefore,  it  is  clearly  not  possible 
to  be  quantitative  regarding  the  directional  properties  of  the  noise  from  a  jet  of  completely 
arbitrary  geometry. 

However,  certain  basic  relationships  were  observed  which  can  be  applied  to  arbitrary 
jet  flows.  The  directivity  of  noise  generated  within  any  region  of  a  jet  will  be  related  to  b/fo 
where  b  is  a  width  parameter  characterizing  the  width  of  the  local  jet  flow  gradient.  For 
constant  value  of  b/x^.1  the  directivity  function  will  sharpen  and  rotote  toward  the  direction 
of  jet  flow  as  the  local  Mach  number  in  the  flow  increases.  For  constant  b/^^  5  and  constant 
local  flow  Mach  number,  the  directivity  function  will  broaden  and  rotote  away  from  the  direc¬ 
tion  of  jet  flow  as  the  temperature  is  increased.  Furthermore,  as  jet  velocity  is  decreased, 
the  noise  radiated  to  the  far  field  decreases  much  foster  than  the  turbulent  fluctuating  pressures 
themselves.  Hence,  the  sound  pressures  along  the  jet  boundary  do  not  decrease  with  velo¬ 
city  at  a  rate  commeasurate  with  the  decrease  in  total  acoustic  power.  A  more  exact  state¬ 
ment  of  theso  conclusions  must  await  the  solution  of  the  jet  noise  equation  developed  in 
Appendix  8. 

It  is  recognized  that  the  statements  in  the  previous  two  paragraphs  represent  the 
weakest  link  in  the  prediction  of  near  field  noise  reduction,  and  clearly  represent  an  area 
where  concentrated  experimental  and  theoretical  work  is  required.  However,  except  for 


regions  within  approximately  one  diameter  of  the  jet,  the  reductions  predicted  for  various 
nozzle  and/or  flow  configurations  based  on  total  radiated  acoustic  power,  are  expected  to 
dominate  for  the  range  of  practical  engine-nozzle  configurations  and  jet  flows.  Therefore, 
the  concepts,  methods  and  conclusions  developed  relative  to  various  noise  suppression 
devices  in  the  remainder  of  this  section  are  considered  to  give  usable  results  for  preliminary 
design. 

The  various  devices  which  have  been  proposed  and/or  developed  for  the  reduction 
of  jet  n^.  rse  act  to  reduce  the  gross  jet  velocity,  shift  the  jet  frequency  spectrum  or  both. 

The  following  subsections  discuss  these  devices  in  terms  of  their  noise  generating  character¬ 
istics  in  an  effort  to  determine  methods  for  predicting  these  characteristics. 

Coaxial  Turbofan  Flows 

Although  velocity  reduction,  per  se,  involves  loss  of  efficiency  and  increased 
engine  weight  flows,  it  is  possible  to  attain  the  objective  of  lower  exit  velocity  with  gains  in 
overall  performance  with  bypass  or  turbofan  engines.  This  principle  was  first  pioneered  by 
Greatrex  of  Rolls  Royce  with  the  Conway  engine,  and  has  been  extended  in  this  country  by 
the  General  Electric  CJ  805-23  "aft  turbofan"  and  the  Pratt  and  Whitney  JT3D  and  JT4D 
"forward  turbofons. "  The  bypass  ratio,  which  is  the  ratio  of  fan  weight  flow  (w$)  to  primary 
weight  flow  (Wp),  varies  for  these  engines  between  1  for  the  Conway  to  approximately  1.5 
for  the  latter  three,  all  at  takeoff  thrust. 

The  modification  of  a  pure  jet  to  a  turbofan  configuration  increases  the  total  thrust 
by  utilizing  part  of  the  energy  ordinarily  exhausted  through  the  nozzle  to  drive  an  additional 
turbine  stage  which,  in  turn,  drives  a  compressor  for  the  bypass  air.  Although  the  thrust  of 
the  primary  jet  is  reduced,  the  additional  bypass  f'ow  more  than  compensates  for  this  reduction, 
resulting  in  a  significantly  greater  net  thrust.  Furthermore,  the  engine  components  can  be 
matched  for  cruise  power  to  achieve  a  lower  fuel  consumption  per  pound  of  thrust  (SFC).  This 
lowered  SFC  reduces  the  weight  of  fuel  required  to  achieve  a  fixed  range  and  more  than 
compensates  for  the  increased  engine  weight  associated  with  the  bypass  stages  and  their  ducting. 
These  performance  factors  are  examined  as  a  function  of  bypass  ratio  for  a  hypothetical  engine 
in  Section  V. 

The  jet  noise  generated  by  a  turbofan  engine  can  be  estimated  directly  from  equation 
(IV-1),  together  with  figures  63  and  64  when  the  two  nozzles  are  coplanor  and  at  the  same 
velocity.  However,  in  most  cases  the  exhaust  velocity  of  the  bypass  flow  will  be  lower  than 
the  primary  velocity,  necessitating  a  more  refined  method  of  noise  estimation. 

Figure  66  gives  an  example  of  an  annular  constant  density  jet  flow  for  a  secondary 
to  primary  velocity  ratio  (Us/Up)  of  .5  and  a  radius  ratio  r$/rp  of  1.5.  This  solution  was 
constructed  from  the  solutions  for  a  jet  exhausting  into  a  uniform  stream  (the  primary  jet)  and 
exhausting  into  ambient  (the  outer  secondary  flow)  from  the  Section  II  and  Appendix  A  solutions. 
The  velocity  profiles  for  this  example  are  further  illustrated  in  figure  67. 

I?  is  clear  from  the  figures  that  foe  primary  jet  is  wholly  surrounded  by  a  uniform  flow 
of  .5  Up  from  foe  nozzle  to  x/rp  equal  to  3.5.  From  x/r0  of  3.5  to  6.3,  foe  surrounding  flow 
velocity  var.es  between  .5  and  .4  times  U,,,  and  so  on.  Thus,  foe  flow  of  the  primary  jet  is 


constructed  for  each  increment  considering  the  average  value  of  the  velocity  in  the 
surrounding  jet.  Therefore,  the  primary  core  is  not  as  long  as  shown  in  the  Section  II 
example  (in  figure  17)  for  a  constant  value  of  velocity  ratio,  but  is  considerably  longer  than 
if  it  were  exhausting  into  the  ambient.  It  is  clear  from  this  example  that  the  resulting 
velocity  profiles  will  be  a  function  of  both  rs/rp  and  Us/Up.  Further,  it  is  also  clear  that 
the  profile  eventually  approaches,  at  some  axial  station,  the  normal  similar  profile  for  a 
single  jet  of  constant  velocity.  The  flow  downstream  of  this  axial  station  will  be  identical 
with  that  of  an  equivalent  single  circular  jet  of  the  same  exit  momentum.  When  the  asymp¬ 
totic  profile  is  reached,  the  momentum  equation  from  Appendix  A  gives  for  the  combined  flow 


x  d-'n. 


where  the  primes  refer  to  the  equivalent  flow. 


As  this  combined  flow  momentum  must  equal  the  exit  momentum  of  the  annular  jets: 

1 rp'tft  b‘  =  Up  K*  +-  ft^  U»  (n\  -  rtp) 

assuming  that  the  nozzle  pressure  ratio  does  not  exceed  1.89. 
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Figure  66.  Estimated  flow  for  constant  denst.y  annular 

coplaner  jet  with  =  .5  and  rj/rp  =  1  -  5 
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Rearranging  gives: 


Since  equation  (I V— 2)  contains  two  unknowns,  b1  and  Um',  it  cannot  be  solved  exactly 
without  an  additional  relationship  which  demands  an  exact  solution  of  the  form  developed 
for  single  flows  in  Appendix  A.  However,  a  useful  approximate  solution  can  be  directly 
obtained  when  rs/rp  is  large  and  (Us/Up)2(rs/rp)2  is  also  large  compared  to  unity.  In  this 
case,  equation  (I  V— 2)  reduces  to 


^m'  _  rp  Us  rs  Um  _  rs 

TV  '  Up  ^0r  W  “  ^ 


(iV-3) 


But  this  is  exactly  the  relationship  for  a  single  jet  with  exit  velocity  of  Us  and  radius  of  rs, 
as  would  be  expected.  Hence,  at  the  tip  of  the  core  of  the  equivalent  jet,  b  =  rs  =  r'  and 
Um'  =  Us  =  IT.  Thus,  under  these  conditions  the  shrouding  flow  becomes  dominant  in  determin¬ 
ing  the  final  jet. 


On  the  othei  hand,  when  rs/rp  approaches  one,  the  primary  flow  must  be  dominant, 
as  is  easily  seen  from  equation  (I V— 2) .  Here,  at  the  core  tip  Um'  approaches  Up  and  b1 
approaches  Tp,  as  expected.  In  order  to  determine  the  equivalent  single  jet  flow  for  intermed¬ 
iate  values  of  rs/ip  and  Us/Up  several  graphical  solutions  were  constructed,  similar  to  figure 
66.  The  resulting  equivalent  jet  parameters  are  summarized  in  figures  68  and  69.  As  shown  by 
figure  68,  r‘/rp  ~  rs/rp  f°r  ’5  4US/Up  the  maximum  deviation  occurring  in  the  range  of 
rs/rp  =  where  r'/rp  =  1-35.  However,  U'/Up  does  net  approach  Us/Up  until  much  larger 
values  of  rs/rp,  except  when  U$/Up  approaches  one. 

It  is  clear  from  the  flow  profile  that  the  total  acoustic  power  generated  downstream  of 
the  siotion  where  the  flow  profile  becomes  approximately  identical  to  that  of  the  equivalenl 
jet,  is  simply  the  power  expected  to  be  generated  by  the  equivalent  jet.  However,  to  estimate 
the  toia-  power  of  the  entire  jet,  it  is  necessary  to  consider  the  noise  generated  by  the  flow 
between  .he  nozzle  and  the  station  where  equivalence  is  definitely  attained.  Referring  to 
figures  66  and  67,  the  initial  flow  is  characterized  by  two  individual  velocity  profiles,  each  of 
total  velocity  drop  of  .5  Up.  Assuming  that  the  noise  generated  in  the  jet  varies  with  the 
eighth  power  of  velocity,  rhe  noise  generated  in  this  region  is  very  low  compared  with  the 
remainder  of  the  jet.  Now  if  the  flow  at  x/rp  =  24  has  reacned  full  equality  with  an  equivalent 
flow  of  U'/Up  =  .87,  and  r'Ao  =  1*35,  it  is  possible  to  compute  the  axial  position  of  the 
equivalent  jet  nozzle  which  gives  the  x/rp  =  24  flow  conditions.  For  x'f/r1  =  13,  and  for  a 
Mach  1  jet  as  in  Section  !l,  tha  axial  position  of  the  equivalent  nozzle  would  be  at 
x/rp  =  24  -  13  x  1.35  =  6.5.  Coincidentally,  this  axiai  position  is  just  beyond  the  region  of 


the  jet  where  the  velocity  gradient  appears  as  two  individual  gradients  and  thus  is  at  the 
beginning  of  the  portion  of  the  flow  which  might  be  expected  to  contribute  the  major  portion 
of  the  noise.  Hence,  as  o  first  approximation  it  would  appear  that  the  noise  estimated  for  a 
total  flow  of  the  equivalent  jet  would  closely  approximate  the  actual  noise. 

This  hypothesis  may  be  directly  tested  against  model  experiments  by  Lee  et  al  (Ref. 2). 
These  tests  consisted  of  varying  the  secondary  velocity  of  an  annular  coplanar  jet  with 
rs/rp-=:1 .5  and  rszz2.65  inches.  The  total  acoustic  power  levels  relative  to  the  power  level 
of  the  primary  flow,  measured  where  U$  =  0,  are  given  in  figure  70  as  a  function  of  Us/(Jp. 

It  is  clear  from  the  figure  that  the  minimum  power  is  approximately  2  db  less  than  that  of  the 
primary  alone  and  occurs  at  a  ratio  of  Us/Up  =  .5.  Considering  that  the  acoustic  power  is 
proportional  to  a  high  power  of  U,  It  is  not  surprising  that  the  maximum  reduction  occurs  ai 
U$/Up  -•  .5.  In  the  early  portion  of  the  flow  where  the  gradients  of  the  two  flows  are  separate, 
minimum  noise  should  result  when  the  difference  in  velocity  between  the  primary  and  secondary 
flow  equals  the  difference  between  the  secondary  flow  and  the  externa!  flow  or  ambient. 
Otherwise,  the  gradient  associated  with  the  greater  velocity  difference  would  dominate  the 
noise  generation  and  would  result  in  a  higher  value  than  attained  with  equal  differences.  The 
relative  values  in  figure  70  were  also  calculated  for  the  equivalent  jet  utilizing  figures  68 
and  69,  together  with  equation  (IV-1).  The  results  of  this  calculation  appear  in  reasonable 
agreement  with  the  data. 


Figure  70.  Variation  of  the  ratio  of  total  acoustic  power  to  power  measured 
with  zero  secondary  flow  velocity  as  a  function  o'  secondary  flow 
velocity  (U$)  for  rs/rp  *  1.5,T  670°  R;  data  from  Lee  et  ah 
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Figure  71  gives  the  power  spectra  for  one  of  the  three  test  series  from  Ref.  2. 

These  spectra  have  been  scaled  from  the  original  data  to  10,000  lbs.  thrust  to  afford  direct 
comparison  with  other  power  spectra  in  this  section.  For  this  purpose,  the  overall  power  level 
scales  directly  with  area  or  thrust,  the  diameter  scales  directly  with  (thrust), and  frequency 
scales  inversely  with  diameter. 

Examination  of  figure  71  shows  that  the  minimum  power  for  constant  thrust  is  attained 
by  the  velocity  ratio  of  .59  which,  of  those  given,  most  nearly  approaches  a  ratio  of  .5. 

Figure  72  compares  the  power  spectrum  calculated  for  the  jet  equivalent  to  the  Us  =  560  ft/s, 

Up  =  951  ft/s  given  in  figure  71 .  Except  for  the  lowest  frequency  band,  the  error  in 
calculated  power  spectrum  does  not  exceed  1  db,  and  the  general  agreement  appears  good. 
Therefore,  it  is  concluded  that  the  equivalent  }et,  defined  for  coplanar  annular  flows  in 
figures  66  and  67,  affords  a  good  prediction  of  the  resulting  jet  noise. 

Figure  73  illustrates  these  results  for  a  fixed  thrust  jet  of  velocity  ratio  Up/Us  =  .5 
in  terms  of  bypass  ratio.  The  ratio  of  bypass  jet  exit  radius  (rs)  to  zero  bypass  jet  exit  radius 
rp]  increases  rapidly  for  bypass  ratios  up  to  1 .5  and  then  increases  slowly  to  its  asymptotic 
value  of  2  at  very  high  bypass  ratio.  Similarly,  the  total  acoustic  power  (W‘)  of  the  equivalent 
jet  decreases  rapidly  for  initial  increases  of  the  bypass  ratio  up  to  ratios  of  3  to  4  and  then 
decreases  more  slowly  to  its  asymptotic  value  of  -18  db. 
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The  class  of  mixing  nozzles  applies  to  nozzles  which  subdivide  the  jet  exit  flow  into 
many  elemental  jets  which  have  a  total  jet  exit  area  approximately  equal  to  the  original 
circular  nozzle  area,  and  distribute  the  elemental  jets  over  a  base  area  larger  than  the 
original  jet.  Several  examples  (Refs.  6,  7,  8,  9,  10,  11,  12)  of  mixing  nozzles  are 
illustrated  in  figure  74.  As  will  be  shown,  the  noise  reduction  of  this  class  of  mixing  nozzles 
involves  both  a  frequency  shift  and  a  lowering  of  the  velocity  of  the  resulting  "mixed  jet. " 


Since  a  propedy  designed  mixing  nozzle  can  be  substituted  for  a  standard  circular 
nozzle,  the  mixing  nozzle  offers  a  possibility  of  reducing  the  noise  of  existing  engines.  Thus, 
the  development  of  mixing  nozzles  has  been  accomplished  by  aircraft  designers  as  well  as 
engine  designers.  Furthermore,  the  development  and  qualification  of  mixing  nozzles  involve 
much  less  hardware  and  time  than  the  development  of  an  optimum  low  velocity  jet  engine. 
Therefore,  great  emphasis  has  been  placed  in  developing  this  class  of  nozzle,  and  intensive 
experimental  research  toward  this  development  has  been  conducted  in  England  and  the  United 
States  by  both  industry  and  government  laboratories  during  the  past  eight  years. 


The  early  interest  in  mixing  nozzles  stemmed  from  Westley's  and  U I  ley's  (Ref.  13) 
observation  that  insertion  of  toothlike  fingers  Into  a  1-inch  cold  jet  resulted  in  a  considerable 
reduction  of  noise.  The  first  full  scale  application  of  this  concept  was  made  by  Greatrex  of 
Rolls  Royce  (Refs.  6,  14,  15),  followed  by  Callaghan,  Howes  and  North  (Refs.  16,  17)  of 
NACA  Lewis  Laboratory.  The  results  of  the  full  scale  experiments  confirmed  that  the  overall 
noise  was  reduced  on  the  o;der  of  3  to  8  db  at  the  angle  of  maximum  radiation  but,  unfortun¬ 
ately,  the  noise  was  increased  toward  the  side  and  forward  of  the  jet.  Consequently,  the 
reduction  in  overall  acoustic  power  was  only  of  the  order  of  1  or  2  db  (Ref.  17),  with  the 
primary  reduction  occurring  at  low  frequencies,  and  a  slight  increase  occurring  at  high 
frequencies  (Ref-  15). 


The  desire  to  minimize  the  thrust  losses  inherent  with  the  teeth  designs,  and 
simultaneously  to  retain  and  maximize  the  noise  reduction  resulting  from  the  teeth  concept 
led  Greatrex  (Ref.  7)  to  propose  and  test  the  corrugated  nozzle.  His  impressive  results  led  to 
the  extensive  investigations  of  corrugated,  segmented  and  tubular  nozzles  by  Boeing,  Douglas, 
Pratt  &  Whitney,  General  Electric,  NACA  and  others  in  this  country,  as  well  as  Rolls  Royce 
in  England.  These  efforts  culminated  in  flight  hardware,  some  of  which  Is  illustrated  in 
figure  74.  Further,  as  a  result  of  these  research  efforts,  a  considerable  body  of  far  field 
data  is  available,  particularly  from  the  NACA  Lewis  Laboratory,  for  a  wide  variety  of 
nozzle  configurations.  However,  before  examining  these  results,  it  is  instructive  to  develop 
a  method  for  understanding  and  predicting  the  fundamental  noise  generating  parameters  for 
these  complex  nozzles. 


The  most  important  and  difficult  problem  in  estimating  the  noise  generated  by  a 
nixing  nozzle  of  arbitrary  geometry  is  the  definition  of  the  jet  flow.  Once  this  definition 
is  obtained,  the  noise  can  be  predicted  to  varying  degrees  of  occuracy,  depending  on  the 
complexity  of  the  flow,  ana  upon  the  noise  prediction  model  which  is  chosen.  Therefore, 
the  following  paragraphs  begin  with  an  examination  of  simplified  flow  models. 
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21 -tube  nozzle  developed 
by  Boeing  for  use  with 
Pratt  &  Whitney  JT-3 
engine  on  707  aircraft 


Eight -lobe  corrugated  nozzle 
developed  bv  General 
Electric  Co.  for  their 
CJ-805-3  engine  installed 
on  Convair  880  aircraft 


Figure  75  is  a  sketch  of  a  section  of  an  axisymmetric  nozzle  which  has  several 
peripheral  tube  nozzles.  The  upper  half  of  the  sketch  illustrates  the  jet  core,  the 
outer  boundary  of  the  jet  and  the  .6  Um  velocity  contour.  The  distance  (b)  between  the 
.6  Um  velocity  contour  and  the  core  boundary  or,  downstream  of  the  core,  the  center- 
line  of  an  individual  axisymmetric  jet,  was  defined  in  Section  II  and  Appendix  A  as  the 
normalizing  factor  in  the  non-dimensional  radial  parameter  (7^  )  where: 


with  r  the  radial  distance  from  the  jet 
centerline  and  a  the  radius  of  the  core. 


U  =  0 


-  U=  .6  U_  . 


at  nozzle  exit 
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Figure  75.  Sketch  of  two  nozzles  in  the  plane  containing  the  axis  of  a 
simplified  axisymmetric  peripheral  tube  nozzle  illustrating 
development  of  final  mixed  jet  and  the  mixing  of  the 
individual  jets. 


The  lower  half  of  figure  75  shows  a  square  velocity  profile  at  the  tube  nozzle 
exit  where  U  =  Ue3  and  the  asymptotic  tip  profile  from  figure  12.  The  details  of  the 
flow  in  the  center  of  this  nozzle  are  deliberately  omitted,  because  it  is  only  possible  to 
define  them  in  a  rather  arbitrary  manner  for  a  few  very  special  cases.  In  this  conneciion 
it  is  noted  that  a  theoretical  solution  for  the  internal  mixing  must  consider  the  radial 
pressure  gradients  in  the  Navier  Stokes  equation,  which,  as  in  Appendix  A,  are  almost 
always  assumed  to  be  zero  in  turbulent  mixing  theories.  Furthermore,  the  geometry  of 
even  the  simplest  mixing  nozzle  is  rather  forbidding  for  rigorous  theoretical  analysis. 

Fortunately,  the  details  of  this  internal  mixing  near  the  nozzle  appear  to  have 
little  bearing  on  the  noise  generation,  except  when  there  is  insufficient  spacing 
between  the  outer  nozzle  elements  to  permit  the  necessary  Inflow  of  ambient  air  to  the 
center  of  the  jet,  for  example,  nozzles  E  and  F  of  reference  18.  When  the  peripheral 
spacing  is  insufficient,  the  static  pressure  in  the  central  region  of  the  jet  is  probably 
lowered  far  enough  below  ambient  to  cause  a  deflection  of  the  outer  jets  toward  the 
center.  When  this  deflection  occurs,  the  internal  flow  mixing  Is  less  complete  than 
would  otherwise  be  anticipated  for  the  geometry  and  the  low  frequency  noise  reduction 
is  reduced. 

For  the  example  of  figure  75,  the  total  exit  momentum  0|Je 3)  is 

^3  =  nj>e3  Ue32  ^  re32 

where  n  is  the  number  of  peripheral  nozzles. 

Assuming  the  peripheral  jets  are  close  enough  to  combine  into  a  larger,  slower 
moving  jet,  denoted  by  subscript  4,  at  some  downstream  station,')T13  ='1914  ,  by 
momentum  continuity.  Assuming  also  that  the  combined  axisymmetric  jet  can  be  repre¬ 
sented  by  the  asymptotic  profile  of  Section  II,  then 

%  ^  (|v_4) 

where  04  is  the  width  of  the  combined  jet  core 
(which  may  be  zero) 
b4  is  the  normalizing  parameter  as 
before,  except  that  it  applies  to  the 
combined  jet. 

To  the  first  approximation,  the  outer  circumference  of  the  peripheral  tube  jets 
will  mix  at  the  normal  rate  with  the  atmosphere.  In  this  case,  the  maximum  velocity 
(Um3)  on  the  centerline  of  individual  tube  jets  will  decrease  at  axial  distances  beyond 
the  tuba  flow  core  tips  until  at  some  axial  station  if  equals  the  velocity  for  the  internal 
or  central  flow.  At  this  axial  station  the  outer  width  parameter  (b3)  for  the  individual 
tube  flows  should  approximately  equal  the  width  parameter  (bij)  for  the  assumed  combined 
flow.  Then  downstream  of  the  cores  formed  by  the  normal  mixing  of  each  peripheral  jet. 


(IV-5a) 


l 
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continuity  of  momentum  (from  Appendix  A)  gives: 

f-rrbiU.V  TT  re32  ue32 . 

For  a  constant  density  jet,  and  at  the  downstream  staiion  where  b3  is  assumed 
equal  to  b4  and  Um3  =  Um4,  a  simple  substitution  in  equation  (IV-5a)  yields: 

^4  Um4  =  re3  Ue3  (IV-5b) 

Further,  by  noting  that  the  centerline  distance  of  the  peripheral  jets  from  the  final 
combined  flow  is  r2  -  rey  then 

a4  #  r2  "  re3  (IV-5c) 


Solving  equations  (IV-4)  and  (IV-5)  gives  for  this  peripheral  tube  nozzle: 

II  .  =  (nl//2  “  1)  Ue3  l.  —  TJ  -  re3 
m4 


1/2 

and  a4  +  b4  =  (r2  -  re3)  (JL- —  ) 

nl/2  -  1 


(IV-6) 


The  noise  radiated  from  this  simplified  model  should  consist  of  the  noise  radiated 
by  the  outer  portion  of  the  individual  tube  flows  before  they  coalesce  at  some  downstream 
station,  plus  the  noise  radiated  by  the  combined  flow  further  downstream.  The  combined 
flow  can  be  considered  identical,  beyond  the  station  where  the  flows  have  combined, 
to  the  downstream  flow  of  a  jet  which  issues  from  a  nozzle  of  radius  re4  =  04  +  b4  at  a 
velocity  Ue4  equal  to  Um4,  as  shown  in  figure  76. 

Some  of  the  many  assumptions  which  were  made  In  the  preceding  simplified  analysis 
can  be  tested  with  the  data  taken  by  Laurence  and  Benninghoff  (Ref.  19)  on  the  flow  from 
the  three  sector  nozzle  illustrated  in  figure  77.  For  the  purpose  of  these  computations, 
the  actual  three  pie-shaped  segment  nozzles  are  replaced  by  three  1.16  inch  radius 
circular  nozzles  which  have  the  same  area  as  the  actual  nozzles.  The  center  of  each 
circular  nozzle  is  assumed  to  act  at  the  centroid  of  each  segment,  as  shown  in  the  top  of 
figure  77.  From  equation  (IV-6)  Um4  is  calculated  to  equal  .47  Ue3,  b4  equals  2.46 
inches,  and  (a4  +  b4)  equals  4.26.  Following  the  assumptions  of  the  derivation  of 
equation  (IV-6),  the  station  where  the  flows  are  expected  to  be  combined  is  located  at 
the  position  where  the  maximum  or  centeriine  velocity  of  each  of  the  three  peripheral 
jets  would  be  expected  to  equal  Um4.  From  figure  16  in  Section  II  the  axial  velocity  of 
a  circular  Mach  .3  jet  of  radius  re  falls  to  .47  of  its  original  velocity  at  x/re  =  29. 

Since  rg3  =  1.16  Inches,  the  jet  should  be  fully  combined  at  on  axial  distance  of  33.6 
inches  from  the  nozzle,  while  upstream  of  this  station  the  outer  portion  of  the  three  jets 
would  be  expected  to  mix  in  the  normal  manner  for  single  jets. 
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The  two  lower  graphs  on  figure  77  show  typical  velocity  contours  measured  by 
Laurence  at  two  downstream  locations  and  compare  the  profiles  based  on  these  contours 
with  the  profiles  calculated  for  the  individual  flows  from  the  relationships  given  in 
figures  12  and  16.  At  the  x  =  8  inches  station  where  x/re3  =  6.9,  there  is  excellent 
agreement  between  the  calculated  and  measured  profiles,  indicating  that  the  outer 
portion  of  the  individual  flows  does  mix  as  expected  for  a  single  jet  throughout  the  core 
region.  At  x  =  24  inches,  where  x/re3  =  20.6,  the  agreement  looks  good  between  the 
velocity  profile  measured  to  the  side,  but  the  effective  center  of  the  individual  flow 
appears  shifted  radially  outward  by  approximately  .5  inch  from  the  assumed  center  given 
by  the  centroid  of  1.8  inches  from  the  center  of  the  entire  nozzle.  Unfortunately,  no 
profiles  ore  given  downstream  of  24  inches  to  determine  axial  station  at  which  the  three 
flows  become  indistinguishable  in  the  combined  flow. 


A  further  confirmation  that  the  turbulence  in  the  outer  peripheral  mixing  zone  is 
similar  to  that  expected  for  the  individual  flow  is  given  in  figure  78.  Here,  the  normal¬ 
ized  mixing  length  {  £ ),  (see  Section  II)  derived  from  the  measured  turbulence  intensity 
and  velocity  gradient,  is  given  as  a  function  of  x/re3,  and  compared  with  Laurence's 
data  (Ref.  20)  for  a  single  jet  shown  in  figure  9.  The  agreement  gives  further  evidence 
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that  the  mixing  of  the  outer  portions  of  the  flows  from  peripheral  nozzles  is  substantially 
in  accordance  with  the  mixing  expected  for  each  individual  nozzle,  at  least  prior  to  the 
station  where  all  the  jets  coalesce.  Therefore,  the  noise  generated  in  this  mixing  region 
should  have  the  frequency  characteristic  expected  for  the  individual  nozzies,  as  given  in 
figures  64  and  65.  Further,  the  total  acoustic  power  generated  can  be  calculated  from 
equation  (IV-1)  when  the  equation  is  corrected  for  the  total  number  of  nozzles  (n)  and 
the  ratio  of  the  exposed  outer  periphery  (Pe3)  of  these  nozzles  to  the  total  periphery  of 
n  single  nozzles.  Thus,  for  computing  the  total  acoustic  power  in  the  initial  mixing 
region,  equation  (IV-1)  becomes 


PWL,  =  146+10  Ioq  n  de32  Pe3  +  80  log 

(nTT  d^)  1000 

=  146  +  10  log  de3  Pe3  +80  log 

^  Tooo 


(IV-7) 

in  db  re  10"^  watts 


It  shculd  be  noted  that  the  noise  generated  by  the  internal  mixing  can  be  neglected 
because  the  internal  turbulence  levels  are  much  lower  than  th^se  on  the  periphery,  as 
shown  by  Laurence  (Ref.  19).  These  lowered  values  of  turbu  ence  would  be  expected 
since  the  total  velocity  difference  between  the  core  of  the  peripheral  jets  and  that  of  the 
internal  flow  is  significantly  iess  than  the  difference  across  the  outer  mixing  zone,  as 
soon  as  the  internal  flow  is  initially  accelerated. 

The  noise  from  combined  downstream  flow  should  be  approximately  equal  to  the 
noise  generated  in  a  similar  downstream  region  by  a  jet  of  exit  radius  equaling  04  +  b4 
and  of  exit  velocity  equaling  Um4,  as  illustrated  in  figure  76.  For  accuracy,  this  com¬ 
putation  should  recognize  that  the  mixing  region  of  the  combined  jet  between  the  equiva¬ 
lent  nozzle  of  radius  04  +  b4  and  the  downstream  station  is  non-existent.  Hence,  the 
contribution  of  noise  from  this  region  to  the  total  predicted  for  the  combined  jet  should  be 
subtracted.  This  correction  involves  the  high  frequency  portion  of  the  acoustic  power 
spectrum  generated  by  the  combined  jet.  Since  for  practical  mixing  nozzles  the  high 
frequency  power  spectrum  of  total  flow  is  dominated  by  the  noise  from  the  individual  tube 
flows,  this  correction  appears  to  be  an  unnecessary  complication  and  can  be  neglected, 
although  the  methodology  for  the  correction  exists  in  this  report. 

These  methods  have  been  applied  to  the  estimation  of  acoustic  power  spectra  for  the 
three  segment  nozzle,  and  the  results  are  given  in  figure  79.  For  this  particular  nozzle, 
the  acoustic  power  from  the  mixed  flow,  which  is  maximum  at  a  low  frequency,  is  con¬ 
siderably  less  than  the  power  generated  by  the  individual  lobes  because  the  velocity  of 
the  combined  jet  is  iess  than  one-half  of  that  of  the  lobes.  Therefore,  the  primary 
reduction  ir.  acoustic  power  is  evident  in  the  lower  frequency  region.  This  computed 
reduction  of  power  spectra  for  two  nozzle  exit  velocities  is  compared  in  figure  80  to  a 
range  of  power  reductions  reported  by  Lee  et  al. (Ref.  2)  for  a  similar  nozzle  over  the 
same  range  ot  velocities.  It  is  interesting  to  note  that  the  computed  frequency  shift  in 
the  acoustic  power  reduction,  which  results  from  the  velocity  factor  in  the  modified 
Strouhal  number  of  figure  64,  generally  agrees  with  the  data. 
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Figure  79.  Calculated  acoustic  power  spectrum  level  for 
3  lobe  2.83-inch  radius  nozzle  with  exit 
velocity  of  1020  ft/sec.  showing  contributions 
of  lobe  flows  and  final  mixed  flow  in  comparison 
with  equal  area  circular  jet  of  2-inch  radius. 
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Figure  80.  Comparison  of  acoustic  power  reduction  calculated 
for  three  lobe  2.83  inch  radius  nozzle  relative  to 
2  inch  radius  circulor  nozzle  with  data  from  Lee  et 
al.  for  five  pressure  ratios 


The  preceding  discussion  and  equations  regarding  an  idealized  peripheral  tube 
nozzle  and  Laurence's  three  segment  nozzle/  were  limited  to  constant  density  jet  flows. 

In  order  to  extend  the  discussion  to  jets  of  practical  interest*,  it  is  desirable  to  consider 
the  gross  effect  of  temperature  and  density.  Although,  as  discussed  in  Section  II,  no 
usable  theoretical  solution  was  found  for  hot  high  Mach  number  jets,  the  data  indicated 
that  the  core  length  of  the  hot  jet  correlated  well  with  the  core  length  of  a  cold  jet  when 
both  were  compared  at  the  same  flow  Mach  number.  Further,  it  was  inferred  that  the 
mixing  rates  were  also  similar,  and  therefore  that  the  constant  density  velocity  profiles, 
mixing  lengths,  etc.,  had  useful  application  in  the  hot  jet  case. 

Recalling  from  the  preceding  discussion  that  the  momentum  radius  of  the  combined 
jet  was  approximately  equal  to  04  +  b4,  then  the  area  (Ae4)  required  for  a  jet  of  constant 
velocity  Ue4  =  Um4  which  has  the  same  momentum  as  the  original  jet  is  approximately 
TF  (a4  +  b4)2.  The  following  brief  derivation  will  give  Um4  as  a  function  of  exit  temper¬ 
ature  ratio  Te]  /Tc  and  area  ratio  Ae]/Ae4,  where  the  subscript  0  refers  to  the  ambient 
surrounding  air  and  1  denotes  conditions  at  the  nozzle  exit  of  a  standard  engine,  without 
mixing  nozzles. 

Assuming  that  the  mixing  of  the  jet  takes  place  at  constant  static  pressure, 

(Pel  =  Pe4  =  Po)'  anc*  perfect  gas  law  applies: 

j|4.Te4  =/=lTc,  =/„T.  <lv-8> 

From  continuity  of  mass  flow,  the  total  mass  flow  rate  (m4)  across  the  plane  perpen¬ 
dicular  to  the  jet  axis  at  the  station  where  the  flows  are  combined  must  equal  the  nozzle 
exit  mass  flow  rate  (me^ )  plus  the  mass  rate  (m0)  of  the  inflow  from  the  atmosphere.  Thus, 
for  the  hypothetical  exit  of  the  combined  jet: 

JIa.  =-vti£4=  vn0+rn\el  (iv-9> 

From  the  conservation  of  heat,  disregarding  conversion  of  turbulent  energy  to  heat, 
and  assuming  the  ratio  of  specific  heats  is  constant  throughout  the  flow, 

me4Te4  =  me|Tel  +  moT0  (IV-10) 

Eliminating  oJq  from  equations  (1V-9)  and  (IV-10)  and  substituting PAU  for  m  gives: 


P  e4Ae4Ue4(Te4  "  ^o)~  _J^el^el  ^el^el  ”  To) 
Substituting  forjP^  from  equation  (IV-8)  gives: 

T  Ue4TQ 

e4  Ue4-ctf  UeO  ->) 

where  CL=  Ae/A^ 
and  p=  T0/Tei 


(I  V— 11) 
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With  this  result,  which  is  summarized  in  figure  81,  it  is  possible  to  estimate  Ue4/Uel, 
once  06  is  determined. 


As  previously  noted,  Ae4  =  IT  (04  +  b4)2  although  neither  04  nor  b4  can  be 
defined  for  the  hot  jet  with  the  precision  developed  for  the  constant  density  jet. 
Fortunately,  the  majority  of  practical  nozzles  have  the  equivalent  of  many  small 
elements,  which  distribute  the  nozzle  exit  area  relatively  uniformly  over  the  area 
enclosed  by  the  envelope  of  the  mixing  nozzle  defined  by  radius  r2*  From  the  last  of 
equation  (IV-6),  which  is  repeated  below. 


re4  =  a4  +  b4  = 


when  re3/r2  1  and  n  '  ^1,  re4  approaches  r2,  the  envelope  radius,  as  would  be 

expected.  Furthermore,  although  b4  and  04  have  not  been  exactly  defined  for  the 
combtoed  hot  jet  flow,  it  is  reasonable  to  expect  that  here,  too,  re4  would  similarly 
converge  to  r2-  Thus,  to  a  first  approximation,  Ae4  can  be  considered  equal  to  A2 
when  the  effective  radii  (re3)  of  the  individual  peripheral  nozzles  are  small  and  the  actual 
exit  area  is  relatively  uniformly  distributed  over  the  envelope  area  (A2).  This  assumption 
enables  a  direct  calculation  of  the  acoustic  power  estimated  for  the  combined  downstream 
jet  flow. 


In  order  to  estimate  the  noise  produced  by  the  initial  mixing  of  the  elemental  jet 
flows,  it  should  generally  be  sufficient  to  consider  only  those  flows  on  the  outer  peri¬ 
phery,  since  all  inner  flows  mix  at  a  lower  turbulence  level.  As  noted  In  the  discussion 
of  figure  65,  the  total  acoustic  power  developed  by  n  separated  elemental  jets,  each  of 
area  Ae3,  is  simply  n  times  the  acoustic  power  generated  by  one  of  the  elemental  jets. 
Thus,  if  the  total  area  of  the  elemental  nozzles,  nAe3,  equals  the  area  Ae]  of  the  original 
single  circular  engine  nozzle,  the  total  acoustic  power  remains  unchanged,  but  the  fre¬ 
quency  spectrum  is  shifted,  as  illustrated  by  figure  65.  When  the  n  elemental  jets  are  in 
close  proximity,  as  in  a  mixing  nozzle  configuration,  their  total  acoustic  power  genera¬ 
tion  in  the  initial  mixing  region  can  be  considerably  less,  because  only  the  peripheral 
mixing  zone  is  of  major  importance  to  the  noise  generation. 


Therefore,  the  total  power  of  the  initial  mixing  process  is  equal  to  the  total  power 
of  the  original  nozzle  of  diameter  de]  times  the  ratio  of  the  effective  outer  perimeter  of 
the  mixing  nozzle  configuration  to  the  total  perimeter  of  the  n  elemental  nozzles.  One 
reasonable  approximation  of  this  perimeter  ratio  R  for  n  nozzles  each  of  radius  re3  is  given 
by: 

r  =  r2/re3  +  •  18  np  "  1  (IV-14) 

n 

where  np  is  the  number  of  elemental  nozzles  on  the  outer  periphery. 
(n>np  »!) 

With  the  above  approximation,  the  total  acoustic  power  radiated  in  the  initial 
mixing  region  for  n  uniformly  distributed  elemental  nozzles  is: 
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PWL3  =  146  +  20  log  de3  +  10  log  nR  +  80  log  Ue3  ,  in  db  re 


10“13  watt 

and  the  total  power  radiated  by  the  combined  jet  flow  is: 


1000 


(IV-15) 


PWL4  =  146  +  20  log  de4  +  80  log  Ue4  ,in  db  re  10_U  watt 


1000 


(IV-16) 


The  spectra  associated  with  these  two  acoustic  powers  are  derived  from  figure  64. 


The  general  results  of  these  considerations  are  clearly  illustrated  in  figures  82  and 
83.  Figure  82  gives  an  example  of  varying  the  envelope  radius  and  hence  the  ratio 
r2 /re  1  for  an  idealized  tube  mixing  nozzle. for  a  10,000  lb.  thrust  engihe  with  1850 
ft/sec.  exit  velocity.  Here  the  diameters  of  the  elemental  tubes  are  held  constant  at 
4  inches.  When  r2Ae]  =  1,  the  elemental  nozzles  are  distorted  and  not  separated,  and 
the  nozzle  is  simply  the  original  circular  axisymmetrio nozzle.  As  the  ratio  ^Ael 
increases,  increasing  the  separation  between  the  elemental  nozzles  and  increasing  the 
area  A2,  the  noise  spectrum  shows  two  maxima.  The  spectral  peak  at  the  lower  frequency 
results  from  the  combined  jet  flow.  It  will  be  noted  that  as  the  ratio  r2Ael  increases 
indefinitely  the  low  frequency  maximum  disappears,  as  expected  for  elemental  jets  whose 
flows  do  not  combine  into  a  single  larger  low  speed  jet.  In  addition,  as  ^Ael  increases, 
the  level  of  the  higher  frequency  maximum,  associated  with  the  peripheral  elemental  flows, 
increases  until  it  reaches  the  spectrum  generated  by  the  summation  of  alii  of  the  elemental 
flows,  each  acting  individually.  Note  that  the  level  of  noise  from  these  mixing  nozzles 
at  frequencies  above  the  upper  frequency  maximum  is  actually  greater  than  the  levels 
generated  by  flow  from  the  original  single  nozzle.  Thus,  the  noise  reduction  of  the  mixing 
nozzles  is  maximum  throughout  the  middle  frequency  region,  a  fact  which  has  plagued 
nozzle  designers  who  desired  to  reduce  the  level  of  high  frequency  noise  reaching  the 
community. 

Figure  83  illustrates  the  effect  of  varying  the  radius  (re3)  of  the  elemental  nozzles 
retaining  the  ratio  r2Ae]  constant  at  1.5.  It  is  clear  that  a  reduction  of  element  radius, 
increasing  the  number  of  nozzles  to  retain  the  necessary  constant  exit  area  in  accordance 
with  figure  65,  reduces  the  level  and  increases  the  frequency  of  the  high  frequency  peak. 
Thus,  for  a  fixed  ratio  of  r2Ael/  the  nozzle  proposed  by  Tyler  (Ref.  4)  which  had  a  very 
large  number  of  very  small  elements  will  give  the  maximum  noise  reduction  at  the  higher 
frequencies. 

These  methods  were  applied  to  the  calculation  of  the  acoustic  power  radiation  from 
the  two  NACA  twelve-lobe  nozzles  sketched  in  figure  74.  The  results  of  the  estimated 
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power  spectra  for  these  two  nozzles  are  compared  with  the  data*  (Ref.  8)  in  figure  84. 
The  agreement  appears  good  in  the  low  frequency  region,  tending  to  justify  the  assump¬ 
tion  that  re4  =  r2  and  the  calculation  of  Ue4  from  this  assumption  and  figure  81 .  The 
high  frequency  maximum  appears  correct  as  estimated  for  the  nozzle  without  centerbody 
but  2  to  3  db  low  for  the  nozzle  with  centerbody.  These  latter  high  frequency  data  also 
suggest  that  the  re3  chosen  for  the  lobes  as  scaled  from  the  figure  in  reference  8  may 
have  been  slightly  too  large,  resulting  in  a  lower  estimate  of  the  center  frequency  for 
the  elemental  flows.  Further,  the  data  suggest  that  the  periphery  ratio  R,  defined  by 
equation  (IV-14  is  too  optimistic  when  the  spacing  between  the  lobes  becomes  large, 
violating  the  original  assumption  of  a  relatively  uniform  distribution  of  the  elemental 
flows  throughout  the  entire  envelope  area.  Nevertheless,  this  idealized  prediction 
procedure  appears  to  be  a  very  useful  too!  In  determining  the  general  noise-generating 
characteristics  of  an  axisymmetric  mixing  nozzle. 


The  results  of  the  preceding  three  figures  indicate  that  the  reduction  in  power 
spectra  from  practical  mixing  nozzles  will  have  fairly  well-defined  limits,  since  the 
envelope  radius  and  minimum  element  radius  re3  are  both  limited  from  practical  con¬ 
siderations.  The  limitations  are  evident  in  the  range  of  power  spectra  for  ten  examples 
of  tube  and  corrugated  axisymmetric  mixing  nozzles  (Refs.  2,  9,  18,  21,  22),  shown  in 
figure  85.  The  nozzles  include  2,  3,  6,  8  and  12-lobe  configurations  and  11,  19,  21 
and  31 -tube  configurations.  The  maximum  power  reduction  which  is  of  the  order  of  15  db 
at  300  cps  was  obtained  with  the  12-lobe  nozzle  with  centerbody,  discussed  with  figure 
84.  Simultaneously,  the  maximum  increase  in  high  frequency  noise  in  figure  85  resulted 
from  this  same  nozzle. 


Figure  86  gives  a  similar  presentation  for  8  nozzles  with  rectangle  envelopes  of 
varying  aspect  ratios  obtained  from  NACA  experiments  (Refs.  9,  18).  For  the  majority 
of  these  nozzles  the  ratio  A2/Ae]  was  held  approximately  constant,  accounting  for  the 
narrow  range  in  the  low  frequency  power  spectra.  However,  two  of  these  nozzles,  E 
and  F,  had  much  poorer  low  frequency  performance  than  the  other  six,  as  shown  in  the 
figure.  For  both  of  these  poorer  nozzles  the  rectangular  elements  had  much  higher  aspect 
ratios  in  comparison  with  the  other  six,  making  It  difficult  for  ambient  air  to  mix  with  the 
interior  flow.  Hence,  it  is  probable  that  the  internal  flow  could  not  gain  sufficient  mass 
flow  to  slow  down  to  the  velocity  expected  from  the  area  ratio  and,  therefore 

the  actual  equivalent  area  (Ae4)  of  the  combined  flow  is  less  than  A2.  This  result 
emphasizes  the  importance  of  providing  sufficient  access  to  the  central  region  of  the  flow 
for  the  necessary  inflow  of  air  from  the  surrounding  ambient  air. 


*The  measured  data  for  the  nozzles  given  In  figure  84  and  other  similar  figures  have  been 
derived  by  applying  the  measured  reductions  of  acoustic  power  as  reported  in  the  litera¬ 
ture,  to  the  generalized  J57  spectrum  which  has  been  utilized  throughout  this  report. 
This  procedure  is  necessary  to  minimize  the  variations  of  ground  attenuation,  for  example 
figure  53,  which  have  affected  almost  all  measurements  of  jet  noise  in  varying  degrees. 
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Figure  84.  Comparison  of  calculated  power  spectra  with  power  spectra  derived 
from  measured  power  reduction  for  two  12-lobe  nozzles  in  NACA 
TN  4261  with  J57  engine. 
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Range  of  sound  power  spectra  for  ten  axisymmetric  mixing  nozzles  of 
lobe  and  tube  type.  Original  data  noise  reductions  have  been 
applied  to  standard  J57  power  spectrum  .Data  from  NASA  Report 
1387  and  TN  4261,  Gonerol  Electric,  Rohr  and  WEAL  files. 


spectrum 


ozzle  F  V  , 
zzle  E^/' 


NAG  A  Report  1387 


Standard  J57 


80  100 


200  300 


500  800  1000 

Frequency,  cps 


5000  8000  10000 


Figure  86.  Range  of  sound  power  spectra  for  eight  segmented 
retangular  nozzle  arrays.  Original  data  noise 
reductions  from  NACA  Report  1387  and  TN  4261 
have  been  applied  to  standard  J57  power  spectrum. 


Acoustic  power  ratio  W/Wref 


The  effect  of  the  spacing  between  the  nozzle  elements  has  interested  several 
investigators,  including  Rollin  of  NASA  (Ref.  23).  His  investigation  was  entirely  with 
a  rectangular  array  which  does  not  exactly  fulfill  the  assumptions  for  the  derivation  in 
this  section.  Nevertheless,  it  is  interesting  to  compare  the  reduction  in  overall  acoustic 
power  predicted  by  these  methods  with  his  results.  This  comparison  is  given  in  figure  87. 
The  results  of  both  prediction  and  experiment  agree  that  a  minimum  in  total  acoustic  power 
should  be  found  for  this  nozzle  array  at  spacing  ratios  between  one  and  two,  and  that  the 
effect  of  the  width  of  the  individual  rectangular  elements  is  slight.  However,  the  predic¬ 
tions  do  not  account  for  the  effect  of  nozzle  pressure  ratio  on  power  reduction  which 
appears  in  the  data  in  the  neighborhood  of  a  spacing  ratio  of  one.  Nevertheless,  the 
general  agreement  is  sufficient  to  indicate  that  these  prediction  methods  are  approximately 
valid  over  a  considerable  range  of  arbitrary  nozzle  geometry. 


Spacing  ratio  s/t 


Figure  87.  Comparison  of  calculated  variation  of  relative  acoustic  power 

with  spacing  ratio  with  measured  values  for  5-segment  rectangular 
nozzle  for  several  values  of  segment  width.  Data  from  Roiiin, 
NASA  TN  D-770. 

The  total  acoustic  power  from  an  idealized  mixing  nozzle  is  the  sum  of  PWI3  and 
PWL4  given  in  equations  (IV— 15)  and  (I  V— 16)  respectively.  The  results  of  these  equations, 
as  demonstrated  in  figure  84,  appear  fairly  valid  when  the  nozzle  flow  is  relatively 
uniformly  distributed  over  the  entire  envelope  area.  Therefore,  it  is  useful  to  apply  these 
equations  to  determine  a  measure  of  the  minimum  overall  acoustic  power  obtainable  by  an 
idealized  mixing  nozzle.  The  results  are  summarized  in  figure  88,  as  a  function  of  the 
ratio  of  envelope  radius  to  the  radius  of  the  original  single  nozzle  for  several  sizes  of 
typical  elements.  It  is  clear  from  the  figure  that  a  relatively  large  reduction  In  sound 
power  can  be  obtained  by  increasing  the  envelope  radius  from  1 .25  to  1*5  times  the 
radius  of  the  original  jet.  In  this  region  of  radius  ratio,  the  overall  sound  power  is  con- 


trolled  by  the  low  frequency  spectrum  of  the  combined  flow,  as  illustrated  In  figure  82. 
However,  further  increases  of  the  radius  ratio  beyond  approximately  1 .5  do  not  result  in 
further  decreases  in  total  acoustic  power  because  the  total  power  becomes  controlled  by 
the  high  frequency  noise  generated  by  the  initial  mixing  of  the  individual  elements,  and 
this  power  increases  as  element  spacing  increases. 

In  addition,  figure  88  shows  that  the  total  acoustic  power  is  minimized  by  minimiz¬ 
ing  the  size  of  the  elemental  nozzles,  which  was  seen  in  figure  83.  Thus,  for  a  range  of 
elemental  nozzle  diameters  between  4  and  2  inches  for  an  engine  of  the  J  57  class,  or 
more  generally  for  a  range  of  t^/t]  between  .18  and  .09  the  maximum  acoustic  power 
reduction  for  an  idealized  mixing  nozzle  varies  between  5  and  8  db,  and  is  obtained  at 
ratios  of  r2/re'j  between  1.25  and  i.5.  Additional  (eduction  in  power  requires  utilization 
of  smaller  elemental  nozzle  sizes. 


1  1.5  2  3  4  5  7  10 


Ratio  of  envelope  radius  to  radius  of  original  single  nozzle  (^Ael) 

Figure  88.  Computed  variation  in  total  acoustic  power  for  ideolized 
mixing  nozzle  relative  to  standard  nozzle  for  Te]  /J0-3, 
os  a  function  of  the  rotio  of  the  mixed  jet  radius  to  standard 
nozzle  exit  radius  with  element  size  os  a  parameter. 

Ejector  Nozzles 

Another  method  of  reducing  the  velocity  of  the  jet  stream,  ond  hence  its  noise 
generation,  is  by  the  use  of  an  ejector  which  shrouds  the  jet.  If  the  ejector  is  sufficiently 
long  and  has  sufficient  entrance  area,  it  will  pump  enough  secondary  air  to  reduce  the 
averoge  velocity  of  the  jet  at  the  ejector  exit  and  ochieve  a  relatively  square  profile. 

If  these  conditions  ore  fulfilled,  the  approximate  exit  velocity  from  the  ejector  con  be 
obtained  directly  from  equation  (IV-13)  or  figure  81,  which  neglect  the  effect  of  the  axial 
momentum  of  the  secondary  air  flowing  into  the  ejector.  Although  for  most  practical 


purposes  this  secondary  air  momentum  is  negligible  since  the  entrance  velocity  is  low, 
equation  (IV-13)  can  be  modified  to  include  the  effect  of  secondary  air  velocity  (Us) 
with  the  following  result,  as  shown  by  Dyer  et  a!  (Ref.  24): 

Unfortunately,  an  ejector  must  be  many  diameters  long  to  obtain  sufficient  mixing 
to  achieve  a  relatively  square  velocity  profile  and  the  resulting  minimum  U4/Ue]  predicted 
by  Q £  ,  the  ratio  of  nozzle  exit  area  to  ejector  exit  area.  Therefore,  the  apolication  of 
relatively  short  practical  ejectors  whose  length  is  of  the  order  of  one  to  two  times  the 
diameter  of  the  primary  jet  to  standard  nozzles  results  in  almost  no  reduction  in  centerline 
velocity.  Hence,  as  shown  by  North  and  Coles  (Ref.  25),  the  corresponding  reduction  of 
total  sound  power  is  insignificant. 

However,  considerably  better  resuits  are  achieved  (Refs.  26,  27)  when  an  ejector 
is  attached  to  a  mixing  nozzle.  In  this  case,  an  ejector,  approximately  two  standard 
nozzle  diameters  long,  shrouds  the  major  portion  of  the  initial  mixing  region.  Since  in 
this  region  the  outer  periphery  of  the  multiple  nozzle  is  greater  than  that  of  the  standard 
nozzle,  there  is  an  increase  of  momentum  exchange  relative  to  the  standard  nozzle  and  a 
corresponding  increase  in  the  efficiency  of  the  ejector  per  unit  axial  length.  Therefore, 
the  velocity  profile  across  the  exit  of  an  ejector  which  is  combined  wirh  a  mixing  nozzle 
approaches  the  idealized  square  profile  much  more  closely  than  does  the  profile  of  the 
same  ejector  when  utilized  together  with  a  'tandard  nozzle.  Consequently,  the  low 
frequency  portion  of  the  acoustic  power  from  the  ejector-mixing  nozzle  combination 
approaches  the  levels  shown  in  figures  82  and  88  where  T2  is  considered  to  be  the  ejector 
radius. 

In  addition  to  the  reduction  of  low  frequency  noise  generation,  the  addition  of  an 
ejector  can  reduce  the  high  frequency  power  generated  by  the  elemental  jets  on  the  peri¬ 
phery  of  the  initial  mixinn  region.  This  reduction  results  from  the  reduced  velocity 
difference  between  the  cores  of  the  elemental  peripheral  jets  and  the  secondary  velocity 
induced  in  the  ejector  intake,  in  accordance  with  the  footnote  to  equation  ( I V—  1 ) .  For 
example,  if  the  induced  secondary  velocity  Us  Is  one-tenth  of  the  primary  exit  velocity, 
the  sound  power  generated  in  the  shrouded  portion  of  the  initial  mixing  zone  would  be 
reduced  by  approximately  f  80  log  .9^  ,  a  4  db  reduction  compared  to  the  power  gener¬ 
ated  in  this  region  without  me  ejector.  To  maximize  the  net  reduction,  it  appears  that 
the  ejector  should  shroud  the  entire  length  of  the  initial  mixing  region,  where  the  high 
frequency  noise  is  generated.  Therefore,  for  a  J57  engine  v.-ith  a  mixing  nozzle  whose 
typical  elements  have  a  two-inch  radius.  It  would  be  desirable  to  shroud  for  an  axial 
distance  approximately  equal  to  the  first  20  -  25  elemental  nozzle  radii.  This  ejector 
would  be  between  40  and  50  inches  long,  or  of  the  order  of  1 .8  to  2.2  times  the  standard 
nozzle  diameter  of  22.5  inches.  Note  that  the  optimum  length  also  depends  on  the  ratio 
between  ejector  radius  and  the  standard  nozzle  radius,  to  assure  that  its  length  covers  the 
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entire  initial  mixing  zone  as  previously  defined  for  mixing  nozzles. 

These  comments  regarding  the  effect  of  the  ejector  on  the  high  frequency  spectra 
are  confirmed  by  the  data  (Refs.  26,  27)  given  in  figure  89  for  four  ejectors  with  mixing 
nozzles.  Note  that  the  high  frequency  acoustic  power  with  ejector  is  less  than  or  equal 
to  that  of  the  standard  engine,  whereas  the  high  frequency  power  of  the  mixing  nozzles 
alone  shown  in  figures  85  and  86  was  generally  higher  than  that  of  the  standard  engine. 
This  reduction  of  high  frequency  power  resulting  from  the  addition  of  an  ejector  to  a 
mixing  nozzle  is  of  importance  to  community  noise  and  was  probably  a  major  considera¬ 
tion  in  the  design  of  the  Douglas  8-lobe  Daisy-ejector  combination  suppression  nozzle. 

The  ejector  has  additional  potential  for  the  reduction  of  the  high  frequency  noise, 
if  its  inner  surface  is  made  to  be  acoustically  absorbent.  For  example,  if  the  ejector 
covered  the  entire  initial  mixing  region,  and  were  capable  of  absorbing  all  of  the  high 
frequency  noise  generated  in  this  region,  shown  in  figures  82,  83  and  88,  the  radiated 
power  spectrum  should  approach  the  idealized  spectrum  for  the  combined  jet  flow  alone. 
Thus,  the  absorbing  ejector  has  the  potential  capability  of  significantly  improving  the 
performance  of  the  mixing  nozzle  beyond  the  limits  shown  in  figure  88.  Furthermore, 
the  shrouding  effect  of  the  ejector  will  always  be  of  major  benefit  to  structure  in  the 
near  vicinity  of  the  exhaust  nozzle. 

Slot  Nozzles 

It  has  been  shown  in  the  mixing  nozzle  analysis,  figure  82,  that  achievement  of 
maximum  low  frequency  reduction  by  replacing  a  single  nozzle  with  a  large  number  of 
small  nozzles  in  an  axisymmetric  array  with  their  higher  characteristic  frequencies, 
involves  a  relatively  large  separation  of  the  nozzles.  Another  basic  configuration  for 
attaining  the  potential  frequency  shift  shown  in  figure  65,  is  to  distribute  the  elemental 
nozzles  in  a  long  line,  perhaps  along  a  wing  trailing  edge.  A  major  step  in  this  direction 
is  embodied  in  the  high  aspect  ratio  slot  nozzle. 

Several  model  scale  experiments  with  slot  nozzles  have  been  conducted  by  NASA, 
Langley  (Refs.  28,  29)  to  determine  both  the  noise  radiation  characteristics  and  the 
possibilities  of  utilizing  the  slot  nozzle  in  conjunction  with  a  wing  flap  to  augment  low 
speed  lift.  Although  no  acoustic  power  data  were  reported,  the  various  spectra  showed 
the  presence  cf  considerably  more  low  frequency  energy  than  would  be  anticipated  from 
figure  65.  Further,  the  results  indicated  that  the  pressure  fluctuations  on  a  jet  augmented 
flap  would  be  sufficient  to  warrant  special  design  of  the  flap  structure  to  avoid  acoustic 
fatigue. 

A  full  scale  experiment  on  slot  nozzles  of  14:1  and  100:1  aspect  ratio  reported  by- 
Coles  (30)  also  confirmed  that  the  low  frequency  noise  generated  by  the  slot  nozzle,  far 
exceeded  that  expected  for  independent  circular  nozzles  of  diameter  equal  to  the  slot 
width.  The  corrected  power  spectra  for  Cole's  100:1  slot  nozzie,  standardized  J57,  and 
idealized  126  two  inch  diameter  independent  nozzles,  ore  summarized  in  figure  90.  It 
is  clear  that  the  high  frequency  power  generated  by  the  slot  nozzle  in  the  vicinity  of  the 
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Figure  89.  Range  of  sound  power  spectra  for  four  mixing  nozzle  -  ejector 

combinations.  Original  noise  reduction  date  from  NASA  TN  4317 
4261  and  D-871  have  been  applied  to  J57  power  spectrum. 


peak  of  the  spectrurp  for  two  inch  nozzles  exceeds  the  power  generated  by  the  standard 
J57.  Considering  only  this  high  frequency  portion  of  the  spectrum,  and  utilizing  the 
periphery  ratio  concept  which  preceded  equations  ( I V—  1 4  and  IV-15),  it  is  noted  that  the 
periphery  of  the  slot  nozzle  (33.7  ft.)  is  approximately  one-half  that  of  the  126  two-inch 
diameter  jets  (66  ft.).  Therefore,  it  would  be  expected  that  the  total  acoustic  power  of 
the  high  frequency  region  would  be  of  the  order  of  3  db  less  than  that  of  the  standard  engine. 
This  expectation  is  approximately  confirmed  by  the  data. 

However,  the  flat  low  frequency  portion  of  the  spectra  shown  by  the  data  is  not 
expected  from  this  initial  rough  analogy.  In  fact,  it  was  this  apparent  inconsistency 
between  data  and  expectation  which  partially  motivated  the  analysis  of  the  flow  from 
slot  jets  in  Section  II  and  Appendix  A.  The  results  of  this  analysis  showed  that  the  mixing 
region  adjacent  to  the  core  of  a  slot  jet,  and  the  length  of  the  core  were  approximately 
equivalent  to  an  axisymmetric  jet  which  had  a  diameter  equal  to  the  slot  width.  However, 
as  shown  in  figure  20,  the  velocity  downstream  of  the  slot  jet  core  decreased  proportional 
to  x“V2,  a  much  slower  rate  of  decay  than  that  of  the  axisymmetric  jet  where  the  center- 
line  decay  Is  proportional  to  x“1 .  Therefore,  downstream  of  the  core,  the  acoustic  power 
per  unit  axial  length  should  decrease  as  x“4  for  the  slot  jet,  as  compared  to  x"7  predicted 
for  the  axisymmetric  jet  from  the  same  dimensional  analysis.  It  is  believed  that  this  lower 
rate  of  axial  velocity  decay  is  responsible  for  the  jet's  generation  of  the  additional  low 
frequency  acoustic  energy  which  results  in  the  rather  characteristic  slot  jet  power  spectrum 
shown  by  the  data  in  figure  90. 

Figure  90  also  contains  a  first  order  estimate  of  the  power  spectrum  of  the  slot  jet. 

This  estimate  utilized  the  power  spectra  from  figure  46,  and  overall  power  per  unit  axial 
length  from  figure  45,  and  the  x"^  relationship,  all  adjusted  for  the  slot  jet  flow  solution 
of  Appendix  A.  The  resulting  prediction  shows  a  large  enhancement  of  low  frequency 
power  relative  to  the  axisymmetric  case,  which  is  generated  for  a  very  great  distance 
downstream.  However,  although  reasonable  agreement  of  prediction  with  the  data  is 
found  above  1000  cps,  the  estimated  low  frequency  power  spectra  are  considerably  lower 
than  the  measurements.  Although  this  discrepancy  could  result  from  turbulence  generated 
internally  in  the  high  flare  slot  nozzle,  it  is  more  probable  that  the  data  represent  a 
typical  slot  noise  spectrum.  Hence,  the  increase  in  low  frequency  power  of  the  slot  jet 
over  the  circular  jet  is  greater  than  the  estimate.  This  result  is  not  surprising  in  view  of 
the  discussion  of  figures  41  and  48,  where  it  was  concluded  that  the  turbulence  level  in 
the  high  Mach  number  flow  does  not  fall  off  as  rapidly  downstream  of  the  core  as  does  the 
centerline  velocity.  Hence,  as  shown  in  figure  45,  the  x”^  law  does  not  apply  until  very 
far  downstream.  Similar  reasoning  applied  to  the  slot  jet  leads  to  the  conclusion  that  the 
reduction  of  power  per  unit  axial  length  downstream  of  the  slot  jet  core  also  is  probably 
slower  than  the  reduction  of  centerline  velocity.  Consequently,  the  x”^  law  utilized  for 
the  estimate  gives  too  much  reduction  of  acoustic  power  generation  with  axiol  distance. 
Therefore,  the  estimated  low  frequency  slot  jet  power  spectrum,  although  exhibiting  the 
proper  trend,  is  too  low  in  level. 

These  results,  although  for  from  providing  a  good  prediction  method  for  the  slot  jet, 
do  serve  to  confirm  the  flat  low  frequency  characteristic  of  the  slot  jet  power  spectrum. 
Consequently,  the  reduction  in  acoustic  power  spectra  aftai noble  with  a  slot  je*  is  nowhere 
near  that  obtained  by  a  large  number  of  nozzles  whose  diameter  is  equal  to  the  slot  width. 
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Figure  90.  Example  of  noise  from  o  100:1  slot  nozzle  comparing  data  from 
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LIST  OF  SYMBOLS  USED  IN  SECTION  IV 


Roman 

Greek 

A 

a 

Area 

Sonic  velocity 

Q L 

^el  /^e4 

a 

(varies  with  T) 

Core  radius 

JS> 

A 

ToAel 

b 

(varies  with  x) 

Gradient  width  parameter 

8 

Mixing  length 

d 

Diameter 

\ 

Radial  parameter  = 

f 

Frequency 

r  -  a 

(0 

on 

Per  cps 

Momentum 

T” 

m 

n 

Mass 

Number  of  nozzles 

\ 

P 

Wavelength 

P 

P 

Periphery 

Pressure 

Density 

PR 

Pressure  ratio 

Subscripts 

PWL 

Power  level  re  10"1 3  watts 

r 

Radius 

* 

At  nozzle  throat 

R 

Perimeter  ratio 

o 

Ambient 

s 

Rectangular  nozzle  spacing 

-1 

Reference  or  standard 

T 

Temperature 

jet  (P.R.  -1.89) 

t 

Rectangular  nozzle  width 

2 

Envelope  radius  of 

U 

Flow  velocity 

mixing  nozzles 

w 

Power 

3 

Characteristic  of 

w 

Weight 

divided  flow 

X 

Axial  distance  from  nozzle 

4 

Characteristic  of 

xt 

Core  length 

e 

combined  flow 

Exit 

Superscripts 

Equivalent 

Time  derivative 

m 

Maximum  at  a  given 
axial  station 

1 

P 

P 

s 

t 

Perimeter 

Primary 

Secondary 

Total 
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SECTION  V 


PERFORMANCE  ASPECTS  OF  NOISE  SUPPRESSION 
NOZZLES  AND  TURBOFAN  ENGINES* 

The  methods  for  reducing  jet  noise  discussed  in  the  preceding  section  included 
mixing  nozzles,  ejectors,  and  slot  nozzles.  These  devices  have  generally  been  designed 
to  attach  to,  and  "quiet,"  existing  engines.  The  other  way  to  reduce  jet  noise,  discussed 
in  Section  IV,  involves  design  of  entirely  new  turbofan  engines.  This  section  discusses  the 
performance  gains  and  penalties  inherent  in  suppression  nozzles  and  a  modern  turbofan 
engine. 

Suppression  Nozzles 

The  losses  in  static  sea  level  thrust  for  mixing  nozzles  which  have  been  quoted 
range  between  1%  and  12%  (Refs.  1,  2,  3,  4,  5,  6)  with  the  better  thrust  nozzles  averaging 
between  1%  and  2%.  However,  the  12-lobe  nozzles  with  and  without  centerbody,  which 
were  among  the  best  noise  reduction  nozzles  and  were  utilized  for  the  example  of  figure 
84,  had  static  thrust  losses  (Ref.  2)  of  6%  and  3%,  respectively.  The  main  reason  for  this 
thrust  loss  Is  probably  the  induced  low  static  pressures  on  the  rearward  canted  areas  of  the 
nozzle.  In  addition  to  thrust  losses  on  the  ground,  the  mixing  nozzles  are  subject  to  drag 
penalties  in  flight  of  the  order  of  1%  to  3%  of  total  aircraft  drag,  resulting  from  base  drag 
losses,  increased  skin  friction  and  interference  between  peripheral  nozzle  elements  at  high 
Mach  number. 

Unfortunately,  the  optimization  of  a  mixing  nozzle  for  acoustic  performance  generally 
degrades  its  aerodynamic  performance.  Maximum  noise  reduction  requires  increased  nozzle 
envelope  area  and  increased  numbers  of  uniformly  distributed  smaller  elemental  nozzles. 

On  the  other  hand,  optimum  thrust  and  drag  performance,  for  other  than  a  circular  nozzle, 
leads  toward  minimum  segmentation,  such  as  a  high  aspect  ratio  slot  nozzle  or  a  4-lobe 
nozzle  arranged  In  a  cross.  Both  of  these  latter  designs  can  be  designed  with  nozzle  thrust 
coefficients  which  are  practically  equal  to  the  standard  circular  nozzle.  Further,  both  can 
be  designed  to  avoid  the  drag  penalties  of  the  lobed  nozzles  caused  essentially  by  "crowd¬ 
ing"  the  boundary  layer  flow  into  narrow,  rapidly  diverging  passages.  The  slot  nozzle 
offers  excellent  possibilities  in  tail-mounted  twin  engine  aircraft  configurations,  with  the 
nozzle  itself  forming  the  trailing  edge  of  the  engine  pylon. 

The  arrangement  of  a  tail  mounted  engine,  Caravelie  style.  Is  aerodynamically  not 
a  simple  task.  The  interference  between  the  body  and  the  converging  parts  of  pylon 
and  nacelle,  forming  diffusers,  under  the  aggravagated  conditions  of  the  very  thick  body 

*The  major  portion  of  this  section  was  contributed  by  R.  Kosin  of  Norair  Division  of 
Northrop  Corporation,  who  was  the  consultant  for  this  program  in  the  area  of  engine  and 
aircraft  performance  and  design. 


133 


boundary  layer  with  local  Mach  numbers  close  to  one,  cause  high  drag,  many  rimes  us  my,, 
as  the  normal  friction  drag.  This  condition  can  be  remedied  by  enlarging  and  flattening 
of  the  "channel"  between  body  and  nacelle.  This  remedy  results  automatically  from  a 
transformation  of  the  standard  nozzle  or  usual  multilobe,  etc.,  nozzle,  to  a  slot  nozzle 
extending  to  the  fuselage.  Additional  improvement  results  from  a  reduction  of  the  pylon 
trailing  edge  which  in  turn  results  automatically  with  the  configuration.  A  drag  estimate 
without  careful  wind  tunnel  tests  is  of  no  real  value  in  a  discussion  of  the  problem,  but  to 
the  experienced  aircraft  engineer,  the  logic  of  the  thought  will  be  convincing  and  if  he 
endeavors  to  improve  an  aircraft,  he  might  use  wind  tunnel  and  analytical  work  to  clarify 
the  problem  to  rhe  desired  degree. 

The  use  of  the  slot  nozzle  in  connection  with  the  jet  flap  principle  has  been  sug¬ 
gested  by  some  authors,  hoping  to  draw  some  additional  benefit  from  the  jet  flap  application. 
Extensive  wind  tunnel  work  by  NASA  has  shown  that  the  jet  flap  shows  worthwhile  gains 
only  with  high  thrust-weight  ratios  at  which  a  reasonable  glide  angle  cannot  be  realized. 

So,  at  least  at  present,  the  slot  nozzle  in  connection  with  the  jet  flap  has  to  be  discarded 
from  practical  considerations. 

The  addition  of  an  ejector  to  a  mixing  nozzle  can  result  in  an  increase  of  static 
thrust  of  the  order  of  1%  to  3%  (Refs.  1 ,  2)  over  the  thrust  of  the  mixing  nozzle  alone.  Thus, 
If  the  mixing  nozzle  has  a  thrust  loss  of  1%  to  2%,  relative  to  the  standard  circular  nozzle, 
the  ejector  offers  the  possibility  of  recovering  this  entire  loss  and  achieving  static  thrust 
values  equal  to  or  greater  than  achieved  with  the  standard  nozzle.  This  thrust  augmentation 
results  from  the  negative  axial  pressure  acting  over  the  nose  of  the  ejector,  resulting  from 
the  acceleration  of  the  secondary  air  flowing  into  the  ejector.  This  advantage  in  thrust  is 
generally  lost  as  soon  as  the  aircraft  has  attained  a  relatively  low  flight  velocity.  Further¬ 
more,  at  higher  aircraft  cruise  velocities,  the  drag  penalties  of  a  fixed  ejector  are  usually 
high,  up  to  the  order  of  a  7%  increase  in  total  aircraft  drag  (Refs.  1,  5).  Therefore, 
ejectors  for  subsonic  and  transonic  aircraft  should  be  retracted  during  cruise.  Note  that  for 
supersonic  cruise  a  properly  designed  ejector  yielding  a  convergent-divergent  nozzle  can 
be  used  in  flight  to  improve  overall  performance.  From  an  overall  performance  viewpoint 
the  combination  of  a  4-lobe  nozzle  with  retracting  ejector  offers  an  attractive  possibility. 
Further,  from  the  acoustic  requirements,  although  the  4-lobe  design  without  ejector  is  far 
from  optimum,  the  4-lobe  nozzle  with  absorbing  ejector  could  give  good  acoustic  perform¬ 
ance,  particularly  In  the  near  field  in  the  vicinity  of  the  nozzle. 

It  is  clea’  from  this  brief  discussion  that  noise  suppression  nozzles  involve  performance 
penalties,  at  best  on  the  order  of  1%  in  all  categories  (Ref.  7),  for  minimal  suppressors. 
Further,  these  penalties  do  not  include  consideration  of  the  reliability  of  complex  nozzle 
geometry,  relative  to  fatigue.  Therefore,  it  Is  desirable  to  examine  the  turbofan  as  an 
alternative  method  of  reducing  jet  noise. 


Turbofan  Performance 


The  results  of  Section  IV  demonstrated  that  large  reductions  of  jet  noise  can  be 
achieved  by  reducing  the  effective  jet  velocity,  utilizing  the  bypass  or  turbofan  engine. 

The  limit  of  these  reductions  is  given  by  machinery  noise,  combustion  noise,  and  turbulence 
developed  in  the  engine  flow  upstream  of  the  nozzle. 

The  noise  of  the  rotating  machinery  can  be  prevented  from  radiating  from  the  front 
of  the  engine  by  a  sonic  throat  in  the  air  intake,  which  allows  inlet  noise  reduction  to  any 
desirable  practical  level.  A  variable  (sonic)  throat  will  be  required  due  to  the  variation  of 
engine  mass  flow  with  flight  speed  and  power  setting.  However,  this  variable  throat  is  of  no 
greater  complexity  than  ordinarily  encountered  with  a  supersonic  inlet,  and  with  some  aero¬ 
dynamic  development  work  the  efficiency  can  be  kept  high  enough  for  all  practical  purposes. 
For  the  part  of  the  machinery  noise  emanating  from  the  rear  end  of  the  engine,  no  such  simple 
and  effective  solution  Is  available.  If  it  is  necessary  to  reduce  this  part  of  the  engine  noise, 
sound  attenuation  by  duct  lining  with  sound  absorbing  material  or  resonators  seems  to  be  the 
only  practical  answer  at  present.  For  sound  attenuation  with  a  lined  duct,  maximum  duct 
length  is  advantageous,  therefore  front  fans  will  be  favorable  to  reduce  noise  radiated  by 
the  fan  to  the  rear.  However,  the  mechanical  solution  of  front  fans  with  bypass  ratios 
greater  than  1.5  or  2  is  not  simple,  because  with  high  bypass  ratios  the  fan  should  run  with 
a  lower  speed  than  the  gas  generator  compressor  to  avoid  excessive  fan  tip  Mach  numbers  or 
too  low  tip  Mach  numbers  for  the  gas  generator  compressor  will  result,  i.e.  at  least  two  if 
not  three  coaxial  shafts  or  a  gear  assembly  will  be  necessary. 

Generally,  it  can  be  stated  that  the  machinery  noise  will  go  up  with  increasing 
bypass  ratio  (secondary  to  primary  mass  flow  rate)  while  the  jet  noise  in  turn  will  be  reduced. 
Therefore,  the  design  of  a  turbofan  engine  with  a  required  noise  level  is  basically  determined 
by  the  selection  of  the  appropriate  bypass  ratio  if  all  other  means  of  sound  attenuation  are 
equally  applied.  Increasing  the  bypass  ratio  of  a  turbofan  engine  requires  an  increase  in 
energy  drop  in  the  turbine  and  therefore  decreases  the  turbine  exhaust  velocity.  At  the  same 
time,  a  lower  turbine  gas  flow  is  necessary  for  an  equal  total  net  thrust,  both  influences 
tending  to  reduce  the  turbine  exhaust  jet  noise.  On  the  other  hand,  an  increased  amount  of 
bypass  air  is  exhausted  at  lesser  speeds.  The  foregoing  statements  hold  generally  true  only  in 
the  subsonic  regime  with  today's  available  turbine  temperatures  of  about  1500°  to  1600°F. 

With  respect  to  noise  in  the  technically  interesting  regime  the  velocity  effect  out¬ 
ranks  by  far  the  effect  of  the  increased  mass,  such  that  with  increased  bypass  ratio  the  jet 
noise  decreases  rapidly.  The  total  engine  noise  would  drop,  too,  if,  as  mentioned  above,  the 
noise  of  the  rotational  machinery  could  be  held  low  or  attenuated  sufficiently.  A  reduction 
of  the  relative  exhaust  speed  below  a  value  of,  say,  600  to  700  ft/sec.  at  static  or  takeoff 
conditions,  Is  considered  useless  in  view  of  the  machinery  noise  (Ref.  8).  This  would  cor¬ 
respond  to  a  bypass  ratio  of  about  6,  thus,  with  respect  to  noise  reduction,  giving  one 
measure  of  a  practical  limit  to  the  bypass  ratio. 


In  order  to  Investigate  the  effect  of  bypass  ratio  on  the  performance  parameters 
of  an  engine  designed  with  today's  technology,  a  hypothetical  engine  of  15,000  lbs. 
static  thrust,  and  3000  lbs.  cruise  thrust  for  Mach  .8  at  35,000  ft.  altitude,  has  been 
developed.  The  performance  of  this  hypothetical  engine  has  been  approximately  computed, 
assuming  the  present  state  of  the  art  for  materials  and  uncooled  turbines,  and  almost 
optimum  engine  pressure  ratios  for  each  case.  The  curves  should  be  considered  as  repre¬ 
sentative  values  only,  because  the  matching  of  a  multi-spool  bypass  engine  is  a  very 
complex  process,  depending  on  a  great  number  of  variables  and  chosen  conditions. 

Figure  91  shows  the  effect  ot  bypass  ratio  on  both  turbine  and  fan  exhaust  velocities 
for  the  two  flight  conditions  at  the  stated  thrusts.  Note  that  although  both  velocities  drop 
rapidly  with  increased  bypass  ratio,  the  turbine  velocity  remains  considerably  above  that 
of  the  fan.  Therefore,  for  bypass  ratios  less  than  approximately  2,  for  which  rs/rp  is  less 
than  approximately  1.5  the  effective  velocity  for  noise  generation  is  close  to  the  primary 
velocities  as  shown  for  constant  density  coaxial  jets  in  figure  69.  Consequently,  consider¬ 
able  acoustic  improvement  would  result  for  the  low  bypass  ratio  engines  if  the  primary 
nozzle  were  designed  as  an  aerodynamical !y  clean  mixing  nozzle,  for  example  a  4-lobe 
nozzle. 

Figure  92  gives  the  net  thrust  per  unit  nozzle  area,  and  figure  93  gives  the  thrust 
per  unit  mass  flow  and  the  nozzle  areas  as  a  function  of  bypass  ratio.  The  influence  of 
bypass  ratio  on  the  power  plant  v/eight  can  be  derived  from  figures  94  through  96.  Figure 
94  shows  the  specific  static  sea  level  thrust,  i.e.  static  sea  level  thrust  divided  by  weight, 
and  the  total  net  thrust  of  bypass  engines  ot  static  sea  level  and  at  Mach  0.8  at  35,000  ft., 
as  compared  to  the  same  values  of  the  basic  gas  generator.  In  figure  95  the  increment  in 
engine  size  and  weight  due  to  bypass  ratio  for  the  conditions  of  equal  thrust  at  Mach  0.8 
at  35,000  ft.  can  be  seen.  Note  that  in  the  comparisons,  equal  cruise  thrust  has  been 
selected  as  the  best  parameter  to  be  held  constant  for  this  study.  Figure  95  also  shows 
the  accompanying  increment  in  static  sea  level  thrust  and  the  decrease  in  primary  (gas 
generator)  air  flow. 

The  total  air  flow  of  the  bypass  engines,  in  comparison  to  the  air  flow  of  the  jet 
engine  with  the  same  cruise  thrust,  is  shown  in  figure  96,  together  with  the  weight  of  the 
installed  bypass  engines  including  cowling  and  jet  nozzles  in  comparison  to  the  jet  engine 
power  plant  of  equal  cruise  thrust. 

To  summarize  the  bypass  engine  analysis,  figure  97  shows  the  weight  expenditure 
for  installed  engine,  including  cowling  and  cruise  fuel,  per  pound  of  cruise  thrust  for 
three  significant  flight  regimes  (medium  range,  transcontinental  and  transatlantic).  The 
optimum  bypass  ratios  are  about  2,  3  and  4,  in  that  order.  As  can  be  seen  from  figure  95, 
the  takeoff  thrust  shows  a  healthy  increase  with  bypass  ratio  if  the  cruise  thrust  (and  the 
throttle  ratio)  is  kept  constant.  Therefore,  the  takeoff  performances  of  the  airplanes  with 
higher  bypass  ratios  will  be  more  satisfactory  from  the  point  of  view  of  safety  and  increased 
altitude  over  nearby  neighborhoods. 
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Figure  93.  Thrust  per  unit  mass  flow  and  nozzle  orea  for  primary  flow  from  the  turbine 
and  secondary  flow  from  the  fan  os  a  function  of  bypass  ratio  for  a 
hypotHeticd  turbofan  engine. 


A  A  . 


V‘*'S'V'vVvv'>; 


i  2 

3  4  5 

6 

y 

i  ■  1  r - 

Bypass  ratio 

- r  "  1  — i — 

- - c — ; 

rr-rrr - m 

A. 

■  .*.*-*  - 


t. 


s'*  .”** 


A\*V- 

.‘•"a'a 

r  - 


t  Am  '  U  *_m 


V.  uVj’ 


139 


Thrust  ratio  (Ib/lb) 


i  r 


Figure  94.  Specific  static  sea  level  thrust  and  ratio  of  thrust  to  thrust  of  basic 
jet  gas  generator  for  constant  primary  fuel  for  hypothetical 
turbofan  engine  as  a  function  of  bypass  ratio. 
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Bypass  ratio 

Engines  with  equal  thrust  at  35,000  ft  altitude,  Mach  No.  0.8 

Relative  size  and  weight,  static  thrust  and  primary  mass  flow  for  hypothetical 
turbofan  engine  which  has  constant  cruise  thrust  as  a  function  of  bypass  ratio. 
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Rgure  96.  Ratio  of  total  oir  flow,  weight  of  complete  engine  including  cowling  and 
nozzles,  and  fuel  consumption  for  hypothetical  fan  engine  to  basic 
turbine  engine  as  a  funct?  n  of  bypass  ratio  with  cruise  thrust  held  constant 
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In  the  study  for  figure  97,  the  increase  in  drag  with  increasing  bypass  ratio,  as 
well  as  the  decrease  in  fuel  consumption  for  takeoff,  climb,  let  down,  and  loiter,  have 
been  neglected.  It  is  fair  to  assume  that  these  influences  will  balance  or  tip  the  scale 
in  favor  of  the  higher  bypass  ratios.  As  can  be  seen  from  figure  97,  a  bypass  ratio  of  3 
is  no  real  penalty  for  the  medium  range  engine,  a  bypass  ratio  of  4  is  truly  acceptable 
for  the  transcontinental  engine  and  5  for  the  transatlantic  engine,  if  so  desired  with 
respect  to  sound  attenuation. 

Note  that  these  engine  data  are  not  statistical  values,  but  data  which  could  be 
achieved  with  engines  to  be  developed  with  the  state  of  the  art  of  1962  -  1964,  and 
which  do  not  have  the  penalty  of  utilizing  existing  hardware. 

The  acoustic  performance  for  a  hypothetical  bypass  engine  derived  from  figures 
91  to  94  has  been  computed  by  the  methods  of  Section  IV,  accounting  for  the  variation 
of  density  with  temperature.  The  resulting  power  spectra  for  an  engine  scaled  at  10,000 
lbs.  rated  thrust  are  given  in  figure  98  for  various  bypass  ratios,  and  the  variation  of 
overall  power  level  with  bypass  ratio  is  given  in  figure  99  for  several  rated  static  thrusts. 
Note  that  these  predictions  do  not  include  the  effect  of  compressor  and  fan  noise,  or 
combustion  noise  and  internal  turbulence.  However,  for  near  field  application  to  the 
sonic  fatigue  problem  the  discrete  compressor  and  fan  frequencies  are  above  the  frequency 
ration  of  concern,  and  any  noise  resulting  from  internal  turbulence  or  combustion  can 
probably  be  adequately  controlled  to  these  power  levels,  if  the  problem  arises. 


Bypass  ratio  5 

Figure  99.  Estimated  variation  of  overall  sound 
power  level  with  bypass  ratio  for 
hypothetical  turbofan  engines  scaled 
to  various  rated  tnrusts. 
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LIST  OF  SYMBOLS  USED  IN  SECTION  V 


r  Nozzle  radius 

w  Weight 

a  Airflow 

p  Primary  flow 

s  Secondary  flow 
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APPENDIX  A 


JET  FLOW 


The  noise  generating  characteristics  of  a  jet  are  directly  related  to  the  jet  flow  and 
mixing  parameters.  Numerous  investigators  (Refs.  1,  2,  3)  have  investigated  the  mixing 
of  a  jet  issuing  from  a  tiny  hole  or  slot  of  very  small  width.  These  solutions  are,  in  general, 
applicable  to  a  finite  jet  in  the  region  well  downstream  of  the  core,  an  axial  distance  of 
approximately  1.5  to  2  core  lengths.  Since  the  noise  from  the  jet  is  primarily  generated  in 
the  mixing  region  along  the  core,  these  downstream  solutions  are  not  directly  applicable. 

A  fundamental  solution  for  the  core  region  of  an  axisymmetric  constant  velocity  jet 
exhausting  into  ambient  air  was  given  by  Kuethe  (Ref.  4)  and  a  solution  was  given  by 
Squire  and  Trouncer  (Ref.  5)  for  an  axisymmetric  constant  density  jet  exhausting  into  a 
moving  airstream.  A  more  recent  report  by  Warren  (Ref.  6)  considers  a  supersonic  axi¬ 
symmetric  jet  of  constant  enthalpy  exhausting  into  ambient  air.  Both  of  the  latter 
investigators  assumed  a  cosine  distribution  for  the  velocity  profiles  which  has  mathematical 
advantages,  whereas  Kuethe  developed  a  more  exact  asymptotic  expression  from  the  basic 
equations  through  an  iterative  procedure.  Kuethe's  profile  fits  the  data  slightly  better 
than  the  cosine  distribution,  but  is  difficult  to  extend  to  other  cases  of  interest. 

These  previous  solutions  all  contain  at  least  one  empirical  constant,  generally  derived 
from  constant  density  flows  of  very  low  Mach  number.  However,  examination  of  various 
flow  data  demonstrates  a  clear  Mach  effect  which  decreases  the  rate  of  jet  spread  for 
increased  Mach  numbers.  Since  the  jets  of  primary  Interest  are  either  sonic  or  supersonic 
at  the  nozzle,  it  is  desirable  to  incorporate  Mach  effects  into  the  flow  solution.  Further¬ 
more,  since  the  effect  of  a  moving  airstream,  either  from  forward  flight  speed  or  shrouding 
secondary  flow.  Is  of  interest  to  the  jet  noise  problem,  it  is  also  desirable  to  consider  this 
effect  together  with  the  Mach  effect.  In  addition,  it  is  also  desirable  to  compwe  the 
difference  between  the  flow  from  a  slot  nozzle  of  finite  width  and  the  flow  of  an  axisym¬ 
metric  nozzle. 

Therefore,  o  systematic  series  of  solutions  was  undertaken  to  include  the  flow  in  the 
core  region  from  the  following: 

a.  Slot  nozzle,  zero  external  velocity,  constant  density  ,  variable  Mach  number 

b.  Axisymmetric  nozzle,  zero  external  velocity,  constant  density,  variable  Mach  number 

c.  Axisymmetric  nozzle,  variable  external  velocity,  constant  density,  variable  Mach 
number 

The  solutions  for  the  core  region  were  extended  downstream  to  approximate  the  entire 
jet  flow  in  order  to  understand  the  differing  low  frequency  noise  generating  mechanisms 
between  slot  and  axisymmetric  jet.  It  should  be  noted  that  complete  similarity  does  not 
begin  until  approximately  two  core  lengths  downstream.  Hence,  Ibis  extension  of  the 
results  from  the  core  is  only  approximately  correct  In  the  transition  region. 
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PREVIOUS  PACE 


The  results  of  the  following  derivations  are  compared  to  available  data  and 
discussed  in  Section  II  of  the  text. 


Slot  Jet  of  Constant  Density 

The  steady  turbulent  flow  for  the  two-dimensional  slot  jet  must  satisfy  the  equation  of 
continuity: 

l£y  +.  3pv  =  o 

and  the  Navier  Stokes  equation: 
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where  constant  pressure  is  assumed  throughout  the  mixing  region  and  viscous  stresses  are 
neglected  in  comparison  to  the  turbulent  stresses.  These  equations  can  be  combined  and 
integrated  with  respect  toy,  (Refs.  5,  7,  8)  to  give: 

*rj  ~uiwf  ^uiti  =  ^  (a— 3) 


Here,  the  first  term  represents  the  change  of  momentum  with  respect  to  x  in  a  one-half 
section  of  the  jet  bounded  by  the  center  plane  and  the  plane  y=  Y  ,  as  shown  in  figure 
100.  The  second  term  gives  the  inward  momentum  across  the  plane  Y  associated  with  the 
change  in  mass  flow  within  the  region  with  respect  to  x,  times  the  velocity  at  plane  Y. 

The  third  term  gives  the  stress  acting  on  the  region  at  plane  Y. 

This  equation  can  be  solved  directly  by  assuming  an  appropriate  velocity  profile  and 
determining  a  suitable  expression  for't'.  The  experimental  data  show  that  the  velocity 
profiles  In  the  mixing  region  adjacent  to  the  core  are  approximately  similar  when  compared 
to  the  parameter 

for  y  ^  a,  and  0  for  y  <  a 


where  a  is  the  outer  edge  of  the  core 
and  b  is  a  width  parameter , 

The  data  also  indicate,  as  shown  in  figure  12  in  the  body  of  this  report,  that  the  velocity 
profile  is  well  approximate^  by  an  exponential  of  the  form 

U  =  Ue.e^ 

This  profile  was  also  advanced  by  Reichardf  (Refs.  1,  9)  for  the  downstream  region  of  on 
axisymmetric  jet  where  complete  similarity  exists.  Because  the  exponential  profile  agrees 
with  the  data  and  is  convenient  to  evaluate  at  the  outer  boundary,  it  has  been  utilized 
for  all  of  the  flow  equations  in  this  series. 


Jet  bour  dary 


Change  in  rate  in  momentum  / unit  width 
due  to  inflow  of  mass  across  1—  2 


Change  in  rate  of  momentum  due  to  shear  stress  = 


Figure  100.  Balance  of  momentum  rates  in  jet  element. 


The  shear  stress  equals- V  where  u  and  are  the  fluctuating  turbulent 
velocities.  PrandiT's  original  mixing  length  hypothesis  (Ref.  1)  gives: 

IL  v  IT  a/  i 

t  |£p|  -f^V^ 

where  H  ,  the  mixture  length,  equals  cb 
and  k  is  an  empirical  constant  of  the  nature  of  a  correlation  coefficient. 

Equation  (A-3)  can  be  directly  evaluated  at  the  outer  boundary,  Y0,  where 


=  Tm  =  0  giving 
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(A-4) 


At  the  nozzle  Y0  =  h/2  where  h  is  the  nozzle  height,  and  U  =  Ufifor  a  rectangular  profile. 
Therefore,  the  constant  is 

Je  Ue  4"  ’ 

which  is  one-half  of  the  momentum  of  the  jet  per  unit  width. 


Substituting  the  velocity  profile  in  equation  (A-4),  and  noting  that  the  outer  boundary 
can  be  chosen  as  •*>  because  of  the  asymptotic  nature  of  the  profile,  gives 
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which  becomes 
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Integrating: 
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Thus,  the  edge  of  the  core  and  all  lines  of  constant  ore  straight  throughout  the 

region  external  to  the  core. 


In  order  to  obtain  the  relationship  between  o,  b  and  x,  equation  (A-3)  is  evaluated  at 
Y  =  a  +  b,  or ^  =  1,  where  the  velocityU^  is  approximately  .61  Ug  .  Substituting  the 
velocity  profiles  and  the  expression  for-f'  in  equotion  (A-3)  gives: 
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Since  kc2  is  assumed  to  be  independent  of  x,  equation  (A-7)  can  be  directly  integrated 
with  respect  to  x,  and  the  constant  can  be  determined  at  the  nozzle  where  b  =  x  =  0  and 
a  =  h/2.  Therefore,  equation  (A-7)  becomes: 


a  +  ,59b  -  -  .95  k  c2x  +  .5  h 


(A-8) 


Simultaneous  solution  of  equations  (A-8)  and  (A-6)  yields 


b  =  3.2  k  c2x 

and  the  tip  of  the  core,  where  a  =  0 
b  =  .  56  h 
xt  .175 

TT  = 

The  constant  (k  c2)  Is  the  single  empirical  constant  required  to  give  the  entire  flow  field 
adjacent  to  the  core.  It  can  be  directly  obtained  from  a  measurement  of  core  length,  or 
evaluated  indirectly  by  an  empirical  determination  of  JL/h  and  k  separately.  These 
approaches  are  both  discussed  in  connection  with  the  axisymmetric  jet  solution,  and  the 
value  of  (k  c2)  is  obtained  from  axisymmetric  jet  data.  The  resulting  flow  field  is 
illustrated  in  the  body  of  the  report  and  compared  with  Laurence's  measurements  (Ref.  8). 

This  solution  for  the  region  adjacent  to  the  core  for  the  constant  density  slot  jet  can 
be  extended  downstream  through  the  transition  region  with  the  assumption  that  the  expon¬ 
ential  profile. 


(A-9) 


(A- 10) 


remains  valid,  and  where  now 7^  =  y/b,  since  <X.=  0,  andU*  is  the  maximum  velocity  on 
the  centerline.  The  solution  requires  determination  of  U^and  b  as  a  function  of  x. 


In  this  downstream  region,  equation  (A-5'  becomes 


Ctf> 


t  Pe  iCe* 


cb{  =  /i  U, 
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Equation  (A-3),  evaluated  at  Y  =  b,  or  '7  =  1 ,  becomes 

J  ^  <^2  —  ~$Y J  ^  ~  •*-5^  ^  ** 

Since  bU^  Is  a  constant  with  respect  to  x,  the  first  term  is  zero  and  the  remaining 
terms  give 

u7  -Tx  CUn  b)  =  •  T  1 5  le  c4 

Substituting  for  b  from  equation  (A-l  1)  and  differentiating; 

hue  =  —  UV1  tec’dLx 

which  integrates  tofjj^A^  2.54  k  c2  ^  +  constant. 


The  constant  can  be  evaluated  at  the  tip  of  the  core  whereL^ssU^by  equation  (A-10) 
giving 


■ 

u_  _ 


|  4- 


»  - 

b  =  -SU  I  ^Z-SAkc.2- 


It  should  be  noted  that  U^/ije.  varies  as  for  the  slot  jet  rather  than  x-^  ft 

the  axisymmetric  jet.  This  extends  the  axial  length  of  the  mixing  zone  which  is  of 
importance  to  low  frequency  noise  generation  further  downstream  for  the  slot  jet  in 
comparison  with  the  axisymmetric  jet. 


Axisymmetric  Jet  of  Constant  Velocity 

The  steady  flow  for  a  turbulent  axisymmetric  jet  must  saiisfy  the  equation  of  continu 
ity  in  cylindrical  coordinates: 


and  the  Novier  Stokes  equation 


rip  U  2  U 
J  3  X 


4-  h-fV 


3  U  _  3  XK 
3  r.  3  n. 


(A- 15) 


where  constant  pressure  is  assumed  throughout  the  mixing  region  and  the 
viscous  stresses  are  neglected  in  comparison  to  the  turbulent  stresses.  Again,  as  in  the 
case  of  the  slot  jet,  these  equations  can  be  combined  and  integrated  with  respect  to  f~L 
to  give: 


9 

9  Y 


R 


(A- 16) 


This  result  gives  the  summation  of  the  axial  changes  of  momentum  within  the  region 
bounded  by  a  conical  surface  R  in  the  same  manner  previously  described  for  equation  (A-3) 
for  the  slot  jet. 

Following  the  development  of  hie  slot  jet,  the  similarity  parameter  utilized  in 
expressing  the  velocity  profile  is: 

\  ~  ^ --  a-  for/L^a,  and  -  0  for/t  <^a 
and  the  velocity  profile  is  given  by  an  exponential 

U  =  u0  -v-  (.ut-u,1)  e  ^  —  Uo+-  u*eT ^ 

where(J©~  the  freestream  velocity  external  to  the  jet. 


The  shear  stress,^,  is  also  of  the  same  form  developed  previously  ,  where 
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du 


din. 
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=  pkcV|^:|  -&■ 


Equation  (A-16)  may  be  directly  evaluated  at  the  outer  boundary  (R0)  af  the  mixing 
region  where  the  shear  stress  =  0  andlLaU  •  giving 

3  r... .  * 

A.  OA,  =•  0 


d  * 


Hence, 
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(A- 17) 
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At  the  nozzle jt^andOKUefor  a  rectangular  velocity  profile.  Therefore,  the 

con!,an,ls  J=’ue^uc-u„)  ^ 

which  is  the  excess  jet  momentum  M,  divided  by  2  IT.  Substituting  the  exponential 
velocity  profile,  equation  (A-17)  becomes 

f(y.  +  i/e^  )  =^p  UeU*-^-  = 

D 


f 


This  expands  to  give 
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Integration  and  rearrangement  yields 
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or 
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(A- 18) 


where 


u. 


u, 


and 


yy)  —  U  C. 


The  relationship  between  a,  b  and  x  can  be  obtained  by  evaluating  equation  (A-16) 
at  R  =  a  +  b,  where 7)  ~  1 .  Substitution  of  the  exponential  profile  and  the  expression  for 
V-16)  gi\ 


T  in  equation  (A-16)  gives 
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which  expands  to 
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Integration  and  substitution  of  p  and  m  as  in  the  previous  equation  results  in 

■Jy^  .l95'*w«?  +  (*22-+-fc4-"p)ab  +(.oti  4- •  OSS -f>) ^7 

=  -  •■51fec.1(a.  +  b)  (A- 19) 

Simultaneous  solution  of  equations  (A-18)  and  (A—  1 9)  is  considerably  more  awkward 
than  in  the  solution  of  the  final  equations  for  the  two-dimensional  jet  since  (b  +  a)  is  not 
necessarily  Independent  of  x.  It  is  therefore  necessary  to  evaluate  the  dependence  of 
(b  +  a)  on  x  for  each  value  of  the  velocity  ratio  (U^J^)  for  which  a  solution  is  obtained. 

The  remainder  of  this  discussion  of  the  constant  density  jet  is  divided  into  two  sub¬ 
sections,  the  first  considering  the  entire  jet  flow  with  zero  external  velocity,  and  the 
second  considering  the  flow  in  the  core  region  with  freestream  velocity. 

Constont  Density  Jet  with  Zero  External  Velocity 

When  the  freestream  velocity  Uo  is  zero‘d  ■-  0,  and  m  =  1,  allowing  a  considerable 
simplification  of  equations  (A-18)  and  (A-19),  giving 


a }  WT  d-b  +  b*  =  /t£ 


(A- 18a) 
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At  the  nozzle,  equation  (A-18a)  reduces  to  a  =  /te,  since  b  =  0,  and  at  the  tip  of  the 
core,  where  a  =  0,  it  reduces  to  b  =  /t-e .  Solution  of  (A-18a)  for  intermediate  values  of 
a  and  b  demonstrates  that  the  value  of  (a  +  b)  remains  within  3%  of  Hoover  the  entire 
range  of  the  core,  as  shown  in  figure  101.  Hence,  equation  (A-19a)  can  be  integrated 
directly  to  give 

a  +  l.i-s  ab  +  .^9b*=-l.9kc(b-Hx)'X+  E 


At  the  nozzle  b  =  x  =  0  and  therefore  the  constant  E  =  a2  • 
a  =  0  and  b  =  /t^and  x  =  xt  so  that 


•  U  tte  -  1.*)  C*-  X-t 


At  the  tip  of  the  core 


Xt  _  .se 

k 


(A-20) 


In  orc’er  to  give  a  comparison  with  the  flow  from  the  slot  jet  downstream  of  the  core, 
the  solution  for  the  constant  density  axisymmetric  jet  with  zero  external  velocity  is 
extended.  In  this  downstream  region  a  =  0  and  hence ^  Furthermore,  the  center- 

line  velocity U**  is  a  function  of  axial  distance  and  the  assumed  velocity  profile  becomes 

u=u„e^ 


kc2x 


Figure  101 .  Variation  of  a  and  a  +  b  with  kc2x  for  constant  density  axisymmetric 
jet  with  a  zero  external  velocity. 
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With  this  substitution/  equation  (A-18)  gives 


4*  X  Xt 


(A— 2 1 ) 


where  b  is  clearly  the  momentum  radius  at  any  axial  station  downstream  of 
the  tip. 

Similar  substitution  in  equation  (A-19)  gives 

^L(/2>i5lCkT) — e^sU*  bx)= — .vifec** b 

(A-22) 

.  t-  i 

The  first  term  is  zero  since  U,*.  b  is  proportional  to  the  total  momentum  which  is 
independent  of  x.  Substitution  of  equation  (A— 21 )  in  (A-22)  gives 

=  I.Sfc  fee1 

which  integrates  to 

b  =  US&  kc’X  4-  M 

The  constant  N  is  evaluated  at  the  tip  of  the  core  where  x  =  Xj.  andU*=  U^,  with  the 
result  that 


m  m 

b  =  I  +  (/***»-)  A* 


(A-23) 


l  +  l-SC,  kcT  (X-zte'j 


(A-24) 


Constant  Density  Jet  with  Variable  External  Velocity 

Equations  (A-18)  and  (A-19)  can  be  solved  for  any  value  of  ^  I  •  When  the 

ratio  '“/'Ue.  =  1/  no  mixing  occurs  and  only  parallel  flow  exists.  These  intermediate 
solutions  are  of  particular  interest  In  the  region  adjacent  to  the  core,  since  they  give  the 
elongation  of  the  jet  core  for  an  aircraft  in  flight,  enabling  an  estimate  of  the  change  in 
noise  source  location  and  character.  Further,  the  solutions  con  also  be  directly  applied  to 
the  mixing  of  concentric  jet  flows. 

The  equations  were  solved  for  both  core  length  and  width  at  the  tip  for  three  inter¬ 
mediate  values  ofU*,  ,  as  given  below.  In  the  solution,  the  variation  of  (b  +  a) 
with  x  was  found  to  be  approximately  linear,  allowing  use  of  a  simple  linear  correction 
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term  In  equation  (A-19)  prior  to  integration. 


Velocity  ratio 

b  at  core  tip 

length  of  core  ; 

0 

i.oo  nt 

.32Ac2 

.25 

.92  /tc 

A8/kc2 

.5 

.82  /t« 

. 83/kc2 

.75 

.76  rte 

1 . 82 Ac2 

1 

0.00  /te 

It  Is  clear  from  the  table  that  the  presence  of  external  velocity  results  in  a  significant 
elongation  of  the  jet  core,  while  at  the  same  time  reducing  its  width.  Thus,  the  noise 
generated  by  a  jet  of  a  moving  aircraft  is  distributed  over  a  considerably  greater  axial 
distance  in  comparison  to  the  static  jet,  resulting  in  gross  distortions  of  the  near  field 
contours. 
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LIST  OF  SYMBOLS  USED  IN  APPENDIX  A 
Roman 

Distance  from  centerline  to  boundary  of  core 
Velocity  gradient  width  parameter 

IA> 

Width  of  rectangle  slot  jet 
Proportionality  coefficient  for  shear  stress 
Prandtl  mixing  length 
Ue/U* 

U0/Ue 

Radius 

rms  turbulent  velocity 

Mean  flow  velocity  along  the  x-axis 

Mean  flow  velocity  along  the  y-  or  r-axis 

Axial  distance  from  nozzle  exit 

Lateral  distance  from  centerline 

Greek 

Nondimensiona!  radial  parameter 

Density 

Shear  stress 

Subscripts 


Maximum  at  a  given  axial  station 
Condition  external  to  jet  flow 
Core  tip 

Superscript 

(u.  -  u<) 


•  \  •*.  **•  A  •*  .*  S  V,  V  VJ 


^  s  ty 


APPENDIX  A 


References 


1.  Schlichting, Hermann,  "Boundary  Layer  Theory,"  Pergamon  Press,  1955. 

2.  Pai,  S.,  "Fluid  Dynamics  of  Jets, "  vanNostrand  Co.,  1954. 

3.  SzablewskI,  W. ,  "The  Diffusion  of  a  Hot  Air  Jet  in  Air  in  Motion,"  NACA 
TM  1288,  1960. 

4.  Kuethe,  A.  M. ,  "Investigations  of  the  Turbulent  Mixing  Regions  Formed  by 
Jets,"  Journal  of  Applied  Mechanics,  pp  A-87  -  95,  1935. 

5.  Squire,  H.  B.,  Trouncer,  J.,  "Round  Jets  in  a  General  Stream,"  RAE  Report 
No.  Aero.  1904,  January  1944. 

6.  Warren,  W.  R.,  "The  Mixing  of  an  Axially  Symmetric  Compressible  Jet  with 
Quiescent  Air,"  Aero.  Eng.  Lab.  Report  No.  252,  Princeton  University, 
September  1953. 

7.  Goldstein,  S.,  "Modem  Developments  in  Fluid  Dynamics,"  Vol.  I,  Oxford 
at  the  Clarendon  Press,  1938. 

8.  Liepmann,  H.  W.,  Laufer,  J.,  "Investigations  of  Free  Turbulent  Mixing," 
NACA  Report  1257,  1947. 

9.  Reichardt,  H.,  "On  a  New  Theory  of  Free  Turbulence,"  Zamm,  Vol.  21, 
No.  5,  October  1941. 


10. 


Laurence,  J.  C.,  "Intensity,  Scale  and  Spectra  of  Turbulence  In  Mixing 
Region  of  Free  Subsonic  Jet,"  NACA  Report  1292,  1956. 


APPENDIX  B 


DERIVATION  OF  MODIFIED  ACOUSTIC  EQUATION  FOR  TURBULENT  JET  EXHAUSTS 

The  purpose  of  this  section  is  to  derive  a  linear  partial  differential  equation  which 
describes,  for  an  axially  uniform  jet  exhaust,  the  production  of  noise  by  the  fluctuating 
turbulent  flow,  and  the  refraction  and  convection  of  the  resulting  acoustic  field  by  the 
mean  flow,  the  radial  mean  flow  gradient  and  the  radial  mean  temperature  gradient.  The 
derivation  proceeds  in  a  straightforward  mathematical  manner  from  the  equation  of  mass 
continuity  and  the  Navier-Stokes  equations  by  considering  the  pressure,  density  and 
velocity  at  each  point  in  the  flow  to  be  a  linear  superposition  of  means,  fluctuating 
turbulent  perturbations  and  minute  acoustic  perturbations.  In  this  analysis  the  acoustic 
"particle  velocities"  and  their  distortions  by  refraction  in  a  high  shear  gradient  flow  are 
not  neglected  as  is  done  in  refs.  1  and  3. 

It  is  assumed  that  sufficiently  accurate  acoustic  field  equations  can  be  obtained 
from  the  linearization  of  the  fundamental  hydrodynamical  equations,  and  that  the  acoustic 
forcing  function  is  associated  with  the  small  error  or  imbalance  obtained  in  describing  the 
mean  and  turbulent  flows,  less  acoustic  perturbations,  by  these  flow  equations.  The  scat¬ 
tering  of  sound  is  not  considered,  and  this  effect  is  suppressed  by  neglecting  terms  in  the 
acoustic  field  equations  which  describe  the  erratic  refraction  and  convection  of  propagated 
acoustic  waves  due  to  time  variations  in  the  flow  velocity,  velocity  gradients  and  tempera¬ 
ture  gradients. 

In  order  to  obtain  a  single  equation  for  acoustic  pressure  or  density  fluctuations, 
these  linearized  field  equations  must  be  combined  in  such  a  manner  as  yo  eliminate  (not 
neglect)  the  acoustic  particle  velocities.  If  the  particle  velocities  are  neglected,  the  low 
frequency  characteristics  of  the  jet  exhaust  cannot  be  accurately  described.  Since  only 
linear  operators  are  involved,  this  step  can  be  done  systematically.  However,  for  the 
general  case.  In  which  axial  mean  flow  gradients  are  retained,  the  resulting  equation  Is 
so  formidable  as  to  be  of  no  practical  value.  When  the  exhaust  is  axially  uniform,  that  is, 
when  the  axial  gradients  of  mean  flow  velocity  and  mean  density  (associated  with  mean 
temperature)  are  neglected,  this  equation  is  considered.  . y  s»  ,ipl?fied  and  appears  to  be 
amenable  to  solution.  Such  an  approximation  is  reaso  :  e  since  the  width  of  the  gradient 
in  the  thin  mixing  region  around  the  core  increases  linearly  along  the  flow  axis  at  an 
included  angle  of  approximately  5  degrees. 

This  equation  is  a  third  order  partial  differential  equation  containing  terms  which 
depend  directly  on  the  mean  velocity  gradient  magnitudes  and  which  influence  both  the 
refraction  of  sound  waves  and  their  production  within  the  turbulent  flow.  An  order  of 
magnitude  analysis  is  presented  which  shows  that  such  terms  are  not  negligible  ot  the 
characteristic  frequency  of  the  jet  and  are  of  greatest  importance  for  the  low  frequency 
sound  waves.  The  terms  do  become  negligible,  however,  for  frequencies  well  above  the 
characteristic  jet  frequency  depending,  of  course,  upon  the  magnitude  of  the  mean  velocity 
gradient.  When  these  terms  are  neglected,  the  acoustic  equation  reduces  precisely  to  that 
given  by  Ribner  in  Ref.  3. 


The  papers  of  Lighthill  and  Ribner  (Refs.  1,  2,  3)  are  important  contributions  to 
the  fundamental  understanding  of  jet  noise  as  viewed  respectively  from  the  points  of  view 
of  both  quadropole  and  simple  source  generation.  The  above  authors  begin  with  the  same 
equation  for  aerodynamical ly  generated  sound,  namely 


*4  -v*P  = 

3-t1 


3*T  ii 
3Xj 


To  =f>AH  arj  +  (  p  -  a0p)£;  ^  -  (tj  +  A^n) 
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but  their  use  of  this  equation  in  defining  the  jet  noise  field  and  estimating  the  sound  levels 
produced  are  fundamentally  different.  LIghthi 1 1  first  notes  the  quadropole  nature  of  the 


forcing  function, 


9X 


Xii 
j  d y  j 


and  takes  an  integral  solution  of  this  equation  in 


the  form  of  Kirchoff's  retarded  potential  solution  of  electromagnetic  theory;  and,  by 
several  mathematical  manipulations,  obtains  a  free  field  volume  integral  of  3'vT\j • 
By  transferring  to  a  moving  "eddy  convection"  coordinate  system,  Lighthill  obtains  the 
motion  correction  factor  (1  -  Mc  cos  0)“3  for  a  frozen  pattern  of  acoustic  quadropoles 
being  axially  convected  by  the  flow.  Then  for  high  shear  gradients,  Lighthill  approximates 


3V(i  /g- 


by  considering  only  the  predominant  mean  shear  amplification  terms. 
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This  approach  leads  to  several  important  explanations  of  the  jet  noise  phenomena;  however, 
this  approximation  of  the  integral  solution  of  the  above  equation  constitutes  only  a  first 
iterative  solution  for  cold  jets  in  which  refractive  gradient  effects  were  not  truly  accounted 
for.  Ribner,  on  the  other  hand,  deals  primarily  with  the  above  differential  equation,  rather 
than  an  integral  equation.  With  the  assumption  that  the  turbulent  flow  is  essentially  imcom- 
pressible,  Ribner  is  able  to  represent  the  sound  field  as  being  produced  by  a  convected 
field  of  simple  sources,  a  conceptually  simpler  approach.  The  equation  obtained  by  Ribner 
shows  more  directly  than  Lighthill's  the  effects  of  this  convection  and  the  refraction  of  the 
sound  waves  by  the  mean  velocity  gradient.  By  considering  special  space  and  time  correla¬ 
tion  functions  for  the  turbulent  flow  and  by  ol lowing  for  "pattern  fluctuations"  of  the 
convected  sources,  Ribner  obtains  an  explanation  of,  and  is  able  to  eliminate,  the  infinite 
amplitude  peak  associated  with  Lighthill's  motion  correction  factor  (  1  -  Me  cos  8)-®. 
Finally,  applying  the  same  restrictive  assumption  as  does  Lighthill,  Ribner  shows  the 
equivalence  of  the  simple  source  and  quadropole  concepts  in  obtaining  the  famous  U®  law 
and  x°  and  x~7  laws. 


in  summary,  Lighthill's  analysis  adequately  predicts  jet  noise  for  low  Mach  number 
cold  flows  where  refractive  effects  are  a  minimum;  in  particular,  as  shown  in  Ref.  2,  the 
theory  gives  a  good  estimation  of  the  noise  generated  by  homogeneous  turbulence  where 
there  is  no  mean  flew.  Ribner's  equation  (Ref.  3}  is  considered  to  be  accurate  for  high 
frequencies  well  above  the  jet  characteristic  frequency  and/or  for  low  shear  gradient  magni- 


164 


tudes.  In  this  regime,  by  allowing  the  speed  of  sound  to  vary  through  the  flow  according 
to  the  variation  of  the  mean  temperatures  of  the  flow,  the  equation  adequately  describes 
both  the  refraction  and  convection  of  the  sound  waves  for  cold  and  hot  flows.  Presently 
there  exists  a  question  as  to  the  effect  of  Ribner's  assumption  of  incompressibility  on  the 
production  of  noise  In  a  hot  flow.  The  modified  acoustic  equation  derived  in  this  appendix 
extends  the  frequency  regime  down  to  and  below  the  characteristic  frequency  of  the  jet, 
and  adequately  describes  the  convection  and  refraction  of  the  sound  waves  propagating 
through  a  heavily  sheared  hot  or  cold  jet  flow.  The  acoustic  forcing  function  is  more 
complicated  than  those  of  Lighthill  and  RIbner,  but  may  simplify  upon  further  analysis  and 
correlation  with  measured  results. 

In  order  to  clearly  show  the  differences  between  the  equations  of  Lighthill  and 
Ribner  and  the  present  equation,  and  in  order  to  fill  in  the  necessary  background,  a  brief 
mathematical  review  of  the  works  of  Lighthill  and  Ribner  Is  presented  preceding  the  devel¬ 
opment  of  the  modified  acoustic  equation. 

Background 


The  fundamental  equations  defining  the  mean  fluid  flow,  turbulence  and  associated 
acoustic  disturbances  in  a  jet  are  the  equation  of  mass  continuity. 
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and’the  Navier- Stokes  equations. 
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(1=1,2,  3),  the  summation  convention  on  j  being  assumed  In  each  case.  The  turbulent 
motion  of  the  fluid  is  initiated  by  flow  instability  at  the  boundary  of  the  laminar  flow 
region  and  the  propagating  acoustic  perturbations  are  in  turn  generated  in  the  turbulence 
by  small  imbalances  between  fluid  inertia  and  internal  fluid  stress.  Thus,  there  are  no 
external  applied  forces  acting  and  neither  are  there  actual  mass  sources  and  sinks,  so  thot 
the  three  force  components,  Fj,  per  unit  volume,  and  the  mass  flow  rate,  Q,  per  unit 
volume  may  be  set  equal  to  zero. 

Equations  (B-l  and  B-2)  may  be  combined  to  give  the  momentum  equations,  by 
multiplying  (B-l)  by  07] and  adding  the  resultant  equation  to  (B-2),  giving 

lrW+  -ft-*** 

(B-3) 

(1=1,2,  3),  the  Fj  and  Q  quantities  being  retained  only  for  generality. 


A  partial  differential  equation  for  density  fluctuations  in  the  fluid  may  be 
obtained  by  first  differentiating  (B-1)  with  respect  to  time, 
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differentiating  (B-2)  with  respect  to  xj,  and  adding  the  three  resultant  equations  for 

i  =  1 ,  2,  3, 
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the  summation  convention  now  being  applied  to  both  i  and  j  ;  then  subtracting  (B-5)  from 


(B-4),  giving 
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To  complete  the  formulation,  the  quantity  CL  V  jP  is  subtracted  from  both  sides  of 
(B-6),  giving  ° 
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(B— 7) 
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Here,  T: x ,  having  the  units  of  stress,  is  coiled  the  stress  tensor,  and  this  quantity  completely 
defines  tne  state  of  internal  compressive  and  shearing  stress  within  the  fluid. 
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Setting  F;  end  Q  equal  to  zero  to  obtain  now  the  proper  mathematical  model  for 
aerodynamical  !y  generated  sound,  (B-7)  becomes 


1L  - 


a„v.P  = 


^  X  i  3  X  j 


(B— 1 0) 


This  is  the  fundamental  equation  derived  and  used  by  Lighthill  in  Ref.  1  to  study  aero¬ 
dynamic  noise  generation-  It  is  to  be  noted  that  (B— 1 0)  is  the  standard  acoustic  wave 
equation,  with  a  forcing  function,  and  thus  describes  the  generation  and  propagation  of 
sound  waves  in  a  uniform  medium  at  rest.  Since  no  assumptions  were  made  in  obtaining 
(B— 10),  the  acoustic  field  defined  by  (B-10)  is  identical  to  the  actual  sound  field  generated 
by  the  real  turbulent  flow. 

As  discussed  in  Ref.  1,  there  are  two  main  advantages  of  approaching  the  jet  noise 
problem  through  (B-10).  First,  the  acoustic  disturbances  of  the  fluid  are  very  weak  relative 
to  the  fluid  disturbances  generating  the  noise  so  that  no  back-reaction  on  the  main  flow  is 
expected  by  the  sound  field,  and  hence  it  is  reasonable  to  consider  the  acoustic  disturbances 
as  produced  by  forced  oscillations  in  the  flu?'  -  Js  shown  in  (B-10).  Secondly,  the  com¬ 
plex  generation,  convection,  refraction  and  _  aring  of  the  sound  waves  by  the  chaotic 
turbulent  flow  can  be  simplified  conceptually  to  the  well  known  propagation  of  sound  waves 
through  a  hypothetical  quiescent  medium  In  which  equivalent  extt*  .'ally  applied  fluid 
stresses  Tjj  act.  These  equivalent  stresses  are  distributed  throughout  the  volume  of  the 
turbulent Tlow  field  and  are  essentially  zero  in  magnitude  outside  this  field. 

With  the  particular  form  of  the  forcing  function  of  (B-10),  Lighthill  (Ref.  1)  was 
able  to  show  directly  the  essentia!  quadropoie  nature  of  the  generation  of  sound  by  the  flow 
field.  Of  the  nine  terms  in  the  forcing  function,^  X  i  ^  j  /  consider  for 

example  the  term  “j“J  ^  X 1 3 X •  This  term  may  be  written  In  the  expanded  form 
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Noting  that  3  TJ^. X  ,  *3  X  2.  is  comparable  to  ^  Q,  /' 3  "t  'n 

(B-7),  it  is  clear  that  each  of  the  four  terms  in  (B-ll)  is  analogous  to  tne  time  rate  of 
change  of  the  mass  flow  rate  per  unit  volume  of  a  simple  acoustic  source.  The  four  simple 
sources  are  separated  by  distances  (A  Xt  )  and  (A  Xt)  lying  In  the  (x],x2J  "  plane,  and 
hence  form  a  lateral  point  quadropole  in  this  plane.  Similarly  for  each  of  the  terms  in 
*3  "T"ij  ^/^Xi^Xj  '  laferc,l  quadropoles  being  obtained  if  j  and  longitudinal 
quadropoles  being  obtained  if  i  =  j. 

Another  principal  advantage  of  (B—  1 0)  is  that  its  solution  may  be  set  down  immedi¬ 
ately  by  use  of  Kirckhoff's  retarded-potential  solution  of  the  wave  equation,  (Ref.  4). 

This  solution  has  the  form 
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where 


t'  =  t--L  (B-13) 
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r=  -  Y\)2  +  (x2  “  yi)2  +  (x3  -  y 3)^  (B_14) 

and  where  represents  the  volume  of  the  turbulent  flow  field,  outside  of  which  Tj j 
equals  zero.  When  there  are  no  reflecting  surfaces  in  the  sound  field  generated  by  the 
flow,  the  term  Ij  if)  equals  zero,  otherwise  its  value  is  non  zero  and  is  defined  by 
the  equation  (Ref.  4), 
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=  direction  cosine 


(B— 16) 


'Yj  =  coordinate  normal  to  the  reflecting  surface 


and  where  the  two  integrals  are  to  be  taken  over  th  entire  surface  of  the  reflector. 


Now  expand  the  following  integral 
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where  the  divergence  theorem  is  used  to  obtain  the  surface  integral,  and  where  by 
definition 
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In  a  similar  fashion,  the  following  Integral  Is  expanded  to  give 
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Substitution  of  (B-18)  into  (B-12)  gives  the  integral  solution  In  the  slightly 
different  but  useful  form  (to  be  used  in  applying  the  motion  correction  discussed  later). 
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Upon  neglecting  reflecting  surfaces  (B— 19)  simplifies  to 
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which  is  the  form  given  by  Lighthill  (Ref.  1). 

Performing  the  differentiation  with  respect  to  x;  and  x;  on  the  Integrals  in  (B-19), 
which  differentiations  obviously  must  include  the  change  off  with  x^,  gives  the  entire 
integral  solution  of  (B-10)  with  reflecting  surfaces  and  in  terms  of  the  time  rates  of  change 
of  the  fluctuating  stresses,  T;:: 
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(B-21) 


This  solution  (B-21),  which  is  more  general  than  that  given  by  Lighthil I  (Ref.  1),  displays 
separately  the  contributions  of  the  turbulent  flow  to  the  far  field,  near  field,  and  near- 
near  acoustic  field.  The  basic  result  of  the  Lighthill  paper  (Ref.  1)  consists  of  the  first 
integral  term  of  (B-21),  namely 


'  fr.-MiVyj  -Mi)  2LJliik^)  Jlv (h) 


(B-22) 


which  defines  the  tota1  radiated  sound  in  the  far  field  of  the  noise  source.  For  distances 
which  are  very  large  relative  to  the  dimensions  of  the  flow  field,  this  equation  can  be 
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further  simplified  to 


PM  =p.+4-^ 


(B-23) 


where 

x2  -  x^  +  x22  +  x32 


(B-24) 


Following  Lighthill's  approach  of  Ref.  1,  a  significant  simplification  of  the 
acoustic  forcing  function  is  made  based  on  the  following  two  physical  arguments.  First, 
the  effect  on  the  propagated  sound  waves  of  the  small  viscous  stresses  is  essentially  one  of 
damping,  both  in  the  region  of  turbulent  flow  and  in  the  uniform  medium  external  to  the 
flow  field,  so  that  the  equation  (B-8)  reduces  to 

Tii  «  pmnri  +(P-  a\  f)  % u 


Secondly,  if  the  flow  is  unheated,  except  for  kinematic  heating  due  to  pressure  fluctuations, 
the  term  (P—cC'op  Is  nearly  equal  to  zero  in  the  turbulent  flow  and  is  exactly  zero  outside 
of  the  turbulent  flow.  Thus,  for  "cold  flows", 

Tii  ^  (9-25) 

these  being  commonly  known  ns  the  fluctuating  Reynold's  stresses. 

The  turbulence  in  a  high  velocity  jet  flow  is  generated  primarily  in  the  thin  high- 
shear  region  between  the  conical  jet  core  and  the  boundary  of  the  jet  flow.  Statistical 
analyses  show  that  the  measured  random  velocity  fluctuations  in  a  fixed  point  in  this  region 
are  partially  time-correlated  for  small  time  intervals,  and  that  the  velocity  fluctuations 
measured  at  two  points,  oriented  in  the  direction  of  the  jet  flow,  are  space-correlated  for 
small  separation  distances.  The  time  correlation  at  the  single  point  can  be  transformed  to  a 
space  correlation  by  application  of  a  mean  flow  velocity;  and,  it  is  found  that  with  the 
proper  choice  of  this  mean  flow  velocity,  the  resulting  space  correlation  c>.<rve  can  be 
made  to  fit  approximately  the  space  correlation  curve  obtained  from  the  two-point  measure¬ 
ments.  It  Is  concluded  then  that  small,  but  finite,  spaciai  regions  exist  within  the  flow 
over  which  the  turbulent  velocity  fluctuations  are  correlated,  and  that  associated  with 
these  regions,  often  called  eddies,  is  a  mean  downstream  convection  velocity. 

It  is  clear  now  that  since  the  fluctuating  turbulent  velocities  are  correlated  over  a 
so-called  eddy  volume,  the  effective  acoustic  forcing  function,  3  ^  * 

must  similarly  be  correlated  over  the  same  volume,  neglecting  large  fluctuations  in  fluid 
density.  Thus  the  volume  distribution  of  quadropole  radiators  within  an  eddy  volume  may 
be  combined  to  form,  approximately,  a  single  eddy  quadropole  radiator  which  is  convected 
downstream  ct  some  convection  velocity,  i.e.  the  velocity  correlation  curve  over  the  eddy 
volume  is  replaced  by  a  square  form  such  that  the  fluctuating  turbulent  velocities  are 
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perfectly  correlated  within  the  eddy  volume  and  perfectly  uncorrelated  with  the  velocity 
fluctuations  outside  the  eddy  volume.  Each  of  the  point-quadropoles  in  an  eddy  then 
fluctuates  in  phase  with  the  quadropole  at  the  center  of  the  volume  and  can  be  replaced 
by  the  center  quadropole  which  is  now  assumed  to  have  a  strength  equal  to  the  average 
quadropole  strength  times  the  volume  of  the  eddy.  This  approach  is  due  to  Lighthf  1 1  (Ref.  1) 
who  further  visualizes  the  jet  exhaust  as  being  composed  of  a  set  of  non-overlapping  eddy 
volumes  of  the  type  just  defined,  and  hence  as  a  discrete  set  of  point-quadropoles  which 
are  convected  downstream  at  seme  average  eddy  convection  velocity.  The  average  eddy 
volume  may  vary  throughout  the  flow,  and  as  Lighthill  points  out  in  Ref.  1,  allowance  can 
be  made  for  a  variable  convection  velocity  in  the  flow. 
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Fixed  observation  point 


-Convected  eddy  volume 
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(Moving  axes) 


(Fixed  axes) 


High  shear  mixing  region 


Figure  102.  A  diagram  showing  the  transformation  from  fixed  (y^)  to  moving  (?£  ) 
coordinates. 

With  this  idealized  acoustical  model  of  the  jet  exhaust  it  is  possible  to  separate 
out  the  effect  on  the  overall  jet  sound  field,  of  the  convection  motion  of  the  quadropole 
generators.  To  do  this,  Lighthill  transforms  the  integrand  of  the  volume  integral  in  (B-20) 
from  the  fixed  yk  -  coordinate  system  to  an^  -  coordinate  system  which  moves  with  the 
quadropole  radiators  at  a  convection  velocity  aQ7^,  where  the  vector"^  is  the  convec¬ 
tion  Mach  number  (see  figure  102).  Consistent  with  the  retarded  potential  type  solution 


in  a  quiescent  medium,  as  given  by  (B-20),  the  two  coordinate  systems  are  easiiy  related 
as  follows:  an  acoustic  disturbance  generated- at  y  at  time  t*  by  fluctuations  in  the  stresses 
Tij(yk't')  W'H  propagate  at  velocity  aQ  from  y  to  the  observation  point  x  in  the  time 

interval—^—  where  r  =  x  -  y  .  If  the  time  of  arrival  at  "x"  is  t,  then  t  =  t1  +  _ . r _  so 

ao  1  ' 

that  the  retarded  time  t1  is,  as  before 


At  this  time  t,  the  position  of  the  element  of  volume  dVfy^)  from  which  the  acoustical 
signal  emanated  is 

~y  -  (a^fXc)S~  or  y  +  "A^  r.  If  7^  denotes  the  fixed 
ao 

position  vector  of  this  same  element  of  volume  referred  to  the#*  -  coordinate  system,  and 
if  it  is  assumed  that  at  time  t  fhe?^-  and  y^  -  coordinate  systems  are  coincident,  then 

=T+M*c  r  (B-26) 

The  stress  tensor  associated  with  the  moving  volume  element  is  then  written.  In  the 
retarded  time  form  of  (B-12)  as 

"£•)“'  V\'a 

The  moving  volume  element  dV  (y[<)  is  transformed  to  the  volume  element  dV(»£  fixed 
in  the?^-  coordinate  system  by  the  Jacobian  of  the  transformation  in  (B-26).  This  Jacobian, 
J,  equals  the  determinant  of  the  transformation  matrix: 
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Evaluating  the  determinant  gives 
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for  a  reflection-free  sound  field,  (B-20)  becomes 
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In  order  to  obtain  a  form  comparable  to  (B— 22),  the  integrand  of  (B-29)  is  differentiated 
with  respect  to  x;  and  xj  assuming  that  the 7^  k  are  fixed.  From  the  definition  of  r 
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and  from  the  coordinate  transformation 
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Combining  (B-30)  and  (B-31)  gives 
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Combining  (B-32)  and  B-33)  gives 


If  (t|»|T|  (B— 37)  can  be  further  reduced  to 
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where  8  is  the  angle  between  the  jet  axis  and  the  observation  position  vector  x  . 

Comparison  of  (B-38)  with  (B-23)  shows  that  the  effect  of  the  quadropole  convection  on 
the  far  field  acoustic  amplitude  is  represented  by  the  factor  (1  -  Mc  cos  0)“3.  Lighthill 
points  out  in  Ref.  1  that  at  low  Mach  numbers  the  difference  between  the  integral  in 
(B-23)  and  that  of  (B-38)  is  simply  one  of  space  variation  due  to  the  motion  of  the 
-  axes  in  the  latter  case.  He  argues  that  this  space  variation  effectively  constitutes  a 
field  of  octopole  generators.  Since  the  radiation  efficiency  of  these  octopoles  is 
significantly  less  than  quadropoles,  the  former  may  be  neglected  and  the  two  integrals  of 
(B-23)  and  (B-38)  are,  to  sufficient  accuracy,  identical  in  magnitude.  Thus  (1-Mj-cos  0)  J 
is  the  only  correction  factor  that  needs  to  be  considered. 

As  a  final  comment  on  Lighthill's  approach  to  jet  noise,  the  importance  of  the 
velocity  gradients  in  amplifying  the  sound  generated  in  the  flow  is  now  considered.  The 
first  time  derivative  of  T  jj  may  be  written  in  the  convenient  form,  noting  the  form  of  T;j 
given  in  (B-25), 

=  TF  ^  ^  3|Q)  -m,  U-,  _|£.  (b_39) 


For  unforced  aerodynamic  flow  in  which  viscous  stresses  may  be  neglected,  the  equation 
of  momentum  (B-3)  may  be  written  In  the  simpler  form 
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3/i 


=  o 


(B-40) 


Also,  (B-l)  can  be  written  in  the  more  convenient  form 
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=  o 


(B-41) 


Substitution  of  (B-40)  and  (B-41)  into  (B-39)  gives 
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which,  upon  combining  and  rearranging  the  derivatives,  gives 
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(B— 42) 


Since  the  last  three  terms  of  (B-42)  are  pure  space  derivatives,  these  represent  octopole 
contributions  to  radiation  field  and  may  therefore  be  neglected.  The  significant  portion  of 
thus  reduces  to 


dT,ydt 
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(B-43) 


where  G  j  j  *s  ca^ed  the  space  rate  of  strain  tensor. 


(B-43)  has  been  derived  using  the  fixed  y^  -  coordinates  or  its  equivalent,  the 
fixed  x|<  -  coordinates.  It  would  seem  more  correct,  however,  to  use  the  moving^  ^ 
coordinates  since  these  appear  in  the  argument  of  Ty  in  (B-38).  In  this  cose 


4-a0Mcic  -SjL£^  n) 
JX>  'd't  ^X* 


(B-44) 


It  is  seen  immediately,  however,  that  the  added  term  is  a  pure  space  derivative  and  also 
represents  an  octopole  contribution  to  the  sound  field  and  may  thus  be  neglected.  This 
agrees  with  lighthill's  conclusion  that  the  volume  integrals  appearing  In  (B-23)  and  (B-38) 
can,  for  all  practical  purposes,  be  considered  equal,  since  their  inteigrands  are  approxi¬ 
mately  equal. 
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Since  high  velocity  jet  flows  into  free  air  are  characterized  by  large  mean  shear 
gradients,  it  is  clear  that  for  a  mean  flow  of  velocity  U  in  the  xj  direction  and  mean 
velocity  gradients  in  the  X2  and  x3  directions,  the  predominant  terms  in  3 Tjj  /  3f  are 
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where  C-  j  j  Is  the  mean  strain  rate  tensor,  or 
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These  two  tensor  elements  can,  for  say  a  cylindrical  jet  exhaust,  be  combined  into  a 
single  radial  term 

-'err  __  P  _2__u 

Bt:  r  9n, 

H  being  the  polar  coordinate  measured  perpendicular  to  the  axis  of  flow. 


Substitution  of  (B— 45)  into  (B-38)  gives 
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It  is  now  clear  that  the  effect  of  a  mean  shear  is  to  amplify  the  pressure  fluctuations 
generated  at  a  point  In  the  flow,  in  proportion  to  the  magnitude  of  the  mean  shear  gradient. 


If  is  nit  desirable  to  carry  the  Lighthill  discussion  further,  as  it  would  be  necessary 
to  introduce  acoustic  intensity  and  total  acoustic  power  and  to  enter  info  a  detailed  consid¬ 
eration  of  statistical  analysis  of  space  and  time  correlations.  For  the  present  purpose,  the 
above  discussion  presents  the  essential  features  of  lighthill’s  approach  to  aerodynamical ly 
generated  sound. 


The  recent  approach  given  by  Ribner  in  Ref.  3  differs  from  Lighthlli  in  several 
respects  and  presents  several  new  and  interesting  features.  Ribner  begins  with  Ligh*hilPs 
equation  (B-10)  above,  with  T;j  defined  as  In  (B— 25),  that  is 
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Expanding  the  right  side  of  (B-49)  gives 
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(B— 50) 


Now,  differentiating  the  equation  of  mass  continuity  (B-l)  with  respect  to  Xj,  multi¬ 
plying  byO/j  and  summing  over  i,  gives 
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Substitution  of  (B-51)  into  (B-50)  gives 
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Now  define  the  operator  D/Dt  as 
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so  that 
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In  terms  of  this  operator  (B-52)  becomes 
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Consider  a  single  volume  element  of  fluid  at  x|<  moving  through  the  flow  field  with 
velocities'll  .  By  constructing  a  coordinate  system^,  which  moves  with  this  fluid 
element,  and  by  associating  (B-55)  with  the  density  and  pressure  fluctuations  of  this 
fluid,  (B-55)  takes  the  form 

J2l£  -Vp-pf-ayi  +9v?  'ami,  tp 

~j[9nii  ini  -271,]  3t  <»-*) 

Ribner  (Ref.  3)  neglects  the  last  term  on  the  right  of  (B-56),  probably  on  the  basis  that 
following  the  motion  of  the  fluid  this  term  should  be  small.  It  is  to  be  noted,  however, 
that  'd'VTj  /gt  represents  the  acceleration  of  the  particle  of  fluid  and  hence  the  accelera¬ 
tion  of  the  moving  axis.  Ribner’s  form  of  (B-56)  is 
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(B— 57) 


The  density, p  ,  pressure  P,  and  velocity /U~f  ,  are  now  written  as  a  linear  superposition  of 
the  ambient  (o),  incompressible  (°),  and  compressible  (’)  components  of  these  quantities. 
(The  following  is  Ribner's  notation.) 

p  P  -4- 
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Substitution  of  (B-58)  into  (B-57)  gives 
*  (*> 
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If  the  fluid  were  incompressible,  then  for  the  single  fluid  element  being  considered, 
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Equation  (B-60)  defines  the  pressure  field  of  a  turbulent  flow  assuming  that  the  fluid  is 
incompressible.  The  condition  of  incompressibility  is,  according  to  Ribner  (Ref.  3),  a 
close  approximation  even  up  to  low  supersonic  Mach  numbers.  Ribner  imposes  this  condi¬ 
tion  on  (B— 59)  and  further  assumes  that  the  bracketed  terms  can  be  neglected  as  a  small 
contributor  to  far  field  acoustic  noise.  Hence,  following  the  motion  of  the  fluid  element, 

1  (j)  ^  p 

-v  p=—  jL J— — 

t z  (BhS1> 

or  referred  once  again  to  fixed  coordinates 
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The  left  side  cf  (B-62)  contains  terms  which  contribute  to  the  scattering  of  sound  waves  as 
they  propagate  through  the  erratic  eddying  flow.  Such  effects,  which  are  not  of  interest 
here,  are  eliminated  by  replacing  D/Dj.  with  D/15^  where 


_D_=-2_ 

Df  9t 
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(B-63) 


Here  U  is  the  flow  velocity  parallel  to  the  jet  axis,  and  it  is  assumed  that  the  axial  gradient 
2iL  can  be  neglected.  With  these  restrictions,  (B-o2)  becomes 


Sl£ 

D  t 


0) 


D1  P'o) 


O  -t2- 


(B-64) 


Dx  PCo> 
D  t* 


(B -65) 


which  is  Ribner's  finol  equation  for  jet  flow  noise.  As  a  first  order  approximation,  Ribner 
further  assumes  that 
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(B -66)  shows  that  the  far  field  acoustic  pressure  perturbations  are  driven  by  the  near 
field  pressure  perturbations,  P(°),  where  the  latter  are  determined  as  if  the  flow  were 
essentially  incompressible  and  result  from  inertial  effects  within  the  unsteady  eddying 
flow.  Further,  (B-65)  shows  that  the  effective  acoustic  sources  are  simple  sources  result¬ 
ing  from  the  minute  pulsations  or  dilatations  of  the  fluid  at  each  point  in  the  flow.  This 
replaces  Lighthill's  volume  distribution  of  acoustic  quadropole  sources. 

T’;e  effective  convection  of  the  acoustic  sources  within  the  jet  can  be  seen  from  the 
expanded  form  of 
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(B— 67) 


When  the  mean  flow  is  zero, 'V/ and Vj  result  solely  from  unsteady  flow  and  represent  the 
dilatation  of  the  fluid  as  it  moves  through  regions  of  rarefaction  and  condensation.  For 
such  flows,  only  the  first  term  may  predominate;  however,  for  high  velocity  jet  flows,  the 
last  term  is  of  the  same  order  as  the  first  term  (Ref.  3). 

It  is  interesting  to  note  that  Lighthill's  convected  quadropoles  were  analyzed  as  a 
"frozen  pattern"  of  acoustic  sources  moving  along  the  flow  axis.  The  resulting  correction 
factor  for  this  convection  was  shown  to  be  (1  -  Mg  cos  8)“3  (Ref.  1).  RIbner,  on  the  other 
hand,  allowed  for  a  pattern  fluctuation  during  convection  through  the  space  -  time  correla¬ 
tion  function,  and  found  that  this  "softens"  the  directional  effect  due  to  convection  and 
eliminates  the  infinite  amplitude  associated  with  Mq  cos  8=  1  in  Lighthill's  convection 
factor. 

Derivation  of  Modified  Acoustic  Equations  for  Jet  Noise 

With  the  above  background  it  is  possible  to  proceed  rather  directly  with  the  deriva¬ 
tion  of  an  equation  which  more  accurately  accounts  for  the  effects  of  high  shear  flow  gradients 
on  the  production,  refraction  and  convection  (the  latter  two  being  inseparable)  of  sound  In 
high  velocity  jet  flows.  The  partied  velocities  associated  with  sound  waves  are  considered 
in  the  fundamental  flow  equations  ana  eliminated  from  these  basic  equations  tc  give  a 
single  equation  for  the  acoustic  disturbances.  This  equation  is  much  like  that  of  Ribner's 
but  Is  of  a  higher  order  and  contains  additional  terms  which  depend  directly  upon  the  shear 
gradient  magnitudes.  Such  terms  are  shown  to  be  Important  at  frequencies  equal  to  and 
below  the  characteristic  frequency  of  the  jet.  In  addition,  a  minimum  of  restrictions  are 
mode  concerning  the  temperature  and  mean  temperature  gradients  in  the  flow. 
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The  forms  of  the  mass  continuity  and  Navier- Stokes  equations  to  be  used  are: 


_4_ p  van 

dt  r  dU  9*i 


rti?  *T^\ 


ZP  =  R 
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Consider  now  an  acoustic  disturbance,  generated  by  Q  and/or  the  F* ,  to  be  propagating 
through  a  steady  flow  field  having  steady  velocity  gradients  and  with  steady  temperature 
gradients.  The  defining  acoustic  equations  are  taken  to  be  the  linearized  forms  of  (B-68) 
and  (B-69).  These  linear  equations  are  obtained  by  neglecting  products  of  acoustic 
quantities  and  neglecting  all  dc  terms.  Thus  denoting  mean  quantities  with  a  bar  (  )/ 

(B-68)  and  (B-69)  become: 


Further,  assuming  a  mean  flow  velocity,/^  =  U,  along  with  the  x]  -  axis  and  setting 
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The  velocities  ''V?  are  the  acoustic  "particle  velocities".  It  is  desirable  to 
eliminate  these  from  (B-72)  -  (B— 75)  in  order  to  obtain  a  single  equation  for  density  and 
pressure  fluctuations.  This  can  be  done  in  a  systematic  manner;  however,  the  resulting 
equation  is  so  formidable  as  to  be  of  no  practical  value.  However,  if  the  mean  flow 
velocity  U  and  the  mean  density  J3  are  axially  uniform,  i.e.  do  not  vary  with  x], 
this  equation  can  be  greatly  simplified. 
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the  above  equations  reduce  to  the  form 
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The  quantities  (  f^n)  can  be  eliminated  from  (B— 77)  -  (B-80)  by  substitution 
into  the  operational  equation: 
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Noting  that 
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the  resulting  equation  is 
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(B— 84) 


(8-85) 


where  P  and  ^  may  be  related  by  the  local  steady  state  speed  of  sound  in  the  flow, 

P  ~~a2^  (B-86) 

For  on  axially  uniform  Steady  flow  flefd  (cold  or  hot),  (B-85)  defines  the  true  production 
of  sound  for  applied  forces  Fj  and  mass  sources  Q,  and  defines  the  true  convection  and 
rcfroction  of  the  generated  sound  through  this  flow.  It  is  seen  that  the  lateral  gradients 
cC  and  have  a  direct  effect  both  upon  the  distortion  of  the  sound  field  and  production 
of  sound,  the  latter  being  evidenced  by  the  enhancement  of  the  axiol  gradients  of  the 
lateral  applied  force  components. 
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When  the  direct  gradient  terms  are  not  significant,  the  left  side  of  (B-85)  reduces 
to  the  left  side  of  (B-64);  the  resulting  equation  is  then 
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In  order  to  relat  this  equation  (B-85)  to  jet  noise,  it  is  necessary  to  determine  F; 
and  Q  in  terms  of  the  turbulent  density,  pressure  and  velocity  fluctuations.  Define 
P 1  Py  /V~]  as  the  linear  superposition  of  mean,  fluctuating  turbulent  and 
acoustic  components: 
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Now,  substitute  (B— 88)  into  (B-68)  and  (B— 69)  and  assume  that  (a)  products  of  acoustic 
quantities  can  be  neglected  due  to  their  relative  smallness;  (b)  products  of  acoustic  and 
fluctuating  turbulent  quantities  can  be  neglected  as  these  contribute  only  to  scattering  of 
the  sound;  (c)  products  of  acoustical  quantities  and  AJ^and  ,  or  derivatives  of  the 
latter,  can  be  neglected  as  are  quite  small  and  do  not  alter  the  refraction  of 

sound;  (d)  the  mean  flow  velocity  U  and  the  mean  density^3  when  associated  with  acoustic 
quantities  do  not  vary  with  xj,  nor  do  their  lateral  gradients  vary  with  xj . 

The  resulting  equations  are  the  same  as  (B -77)  -  (B-80)  but  with  F|  and  Q  defined 
as  follows: 


F  -  —  P  DAr'  —  9P 
r‘  ~  " p*  9Xi 

Q  =  -r* 

-iff 
=  "If 


P=  P  +  P 

at,  —  a77  +  aF\ 


(B— 89) 


In  performing  the  substantive  derivative  operation, - ,  in  (B-81),  it  is  reasonable  to  use 
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the  more  general  operation  ,  for  F:  and  Q,  since  the  use  of -  on  the  left  side  of 
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(B -77)  -  (B-80)  merely  suppresses  the  scattering  of  propagated  sound  which  Is  not  the  case 
for  Fj  and  Q.  Hence,  the  desired  equation  defining  jet  noise  is 
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with  F:  and  Q  defined  by  {B-89).  For  comparison  with  Ribner's  equation  (B~d5),  the 
quanti  J2S  -  mil,  determined.  From  (B-89) 
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Ribner  obtained  the  first  bracketed  terms  on  the  RHS  of  (B-90)(see  B-65,  B-67)  and 
assumed  complete  cancellation  of  the  second  bracketed  term  with  •  Th°s»  "’'v 


difference  lies  in  the  term 
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which  arises  from  eliminating  the  particle  velocities. 


Because  of  the  small  mc^nitude  of  the  dilatation  0,  and  because  D/VDtis  proportional  to 
Co  rather  than  to  Oj\  os  in  3>/9t\  this  added  term  is  relatively  small  and  may  be 
neglected.  Accepting  the  above  cancellation,  which  seems  sufficiently  accurote  for  the 
present,  (B-90)  becomes 
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can  be  estimated  by  comparing  them  with  — -  _ 5  •  Let  a  typical  Fourier  component  of 
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Let  M  be  an  averoge  flow  Mach  number,  say  1/2  Mq,  and  let^^  *  ^  where  R  is  the 

R 

radius  of  the  exhaust  exit  and  cu  =  2  IT  f.  Then 
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For  an  exit  velocity  of  1020  ft/sec.,  a  flow  Mach  number  of  =  1.0  and  o  nozzle  radius 
of  .5  ft.,  this  quantity  equals 

120 

? 

The  above  ratio  of  absolute  magnitudes  equals  unity  for  f  =  120  cps  and  is  of  the 
order  of  1/2  at  f  =  240  cps,  which  for  the  example  given  is  approximately  the  characteristic 
frequency  of  the  jet.  The  two  additional  grodient  terms,  therefore,  should  not  be  ignored. 

The  solution  of  (8-91)  is  presently  being  Investigated  by  the  author  for  various  types 
of  acoustic  sources,  such  as  simple  sources,  dipoles,  etc.,  and  for  various  gradient  shopes 
and  flow  Mach  numbers. 


LIST  OF  SYMBOLS  USED  IN  APPENDIX  B 


«a 


J 

> 

ml 


r 


^  I 


D/bt 

D/Dt 


Speed  of  sound 

Total  differentiation  operator  =  'fy/d't  +  ^  'd/'dXj 
Total  differentiation  operator  =  'i/'it  +  u  , 

Strain  rate  tensor  =  4-  ^ ^0/3 X i 

Frequency,  cps 
Applied  force  component 
Surface  integral;  (Equation  B— 15) 

Jacobian  of  transformation  from  fixed ^  -  to  moving -  coordinates 
Mach  number 

Coordinate  along  surface  normal 
Pressure 

Mass  source  or  sink  flow  rate 

Distance  between  fluid  element  at  ij  and  observation  point  at  X  *  j  X  - 
Cylindrical  polar  coordinate 


*-m|  -mI*  (7-7) 


Surface  area 


Retarded  time  in  uniform  descent  medium  at— 


Fluid  stress  tensor 


Mean  jet  flow  velocity 


Volume 


Fluid  velocity  component 


*V*V«V‘  *'•  .**  '  •  '*,  .'.v  V/ 


%J  sb  Position  of  fixed  observation  point  with  respect  to  fixed  axes 

ijj  =  Position  of  moving  element  of  fluid  with  respect  to  fixed  axes 

c(,  =  Lateral  mean  velocity  gradient  =  3  U/3^^ 

^  =  Lateral  mean  velocity  gradient  =  3  1^/3  Xj 

ss  Kronecker  delta  =  0  if  i  i j,  =  1  if  i  =  j 
"t[  =s  Shear  viscosity  coefficient 

7^'  —  Compressive  friction  coefficient 

ss  Position  of  moving  element  of  fluid  with  respect  to  moving  axes 
0  =  Dilatation  =  Xj  ;  also,  angle  between  jet  axis  and  observation  vector 

f  —  Fluid  density 


U)  — 


Circular  frequency,  rad/sec. 

i/’i  X,1  +  b/’S  <1  +  3/3  x; 


A  =»  Finite  difference  operator 


Superscripts 


Denotes  near  field  component 
Denotes  far  field  component 
Denotes  mean  or  time  average 
Denotes  vector  quantity 
Denotes  fluctuating  turbulent  component 


Subscripts 


i,  j,  k  Denote  components  of  vector  quantities  along  the  1,  2,  3  axes* 
s  Denotes  surface 

•Note:  The  1  -  axis  is  always  parallel  to  the  mean  flow  axis  of  the  jet  exhaust  for 
both  the  fixed  and  moving  coordinate  systems. 
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Subscripts  (continued 


Denotes  ambient  value 


Denotes  acoustic  perturbation  value 


Denotes  convection 


Denotes  absolute  magnitude 
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APPENDIX  C 


SOUND  PROPAGATION  IN  A  CYLINDRiCALlY 
SYMMETRIC  TEMPERATURE  DISTRIBUTION 


In  an  attempt  to  extract  the  qualitative  features  at  least  of  jet  noise  propagation 
at  the  presence  of  a  temperature  distribution  (hot  jet),  the  propagation  of  sound  from  a 
source  imbedded  in  a  finite  cylindrical ly  symmetric  temperature  distribution  will  be  con¬ 
sidered  here.  (The  atmosphere  is  assumed  to  be  otherwise  homogeneous.)  The  particular 
choice  of  the  distribution  is  dictated  by  that  existing  in  a  hot  jet,  which,  in  general, 
follows  a  bell-shaped  curve.  * 

If  the  z-axis  is  the  axis  of  symmetry,  then  the  temperature  will  depend  only  on  the 
coordinate  fc  in  particular,  the  chosen  dependence  is 


T{n)  «  T“o 


_A _ 

CO?h*dC/L 


(C-l) 


where 


To  -  ambient  (atmospheric)  temperature  -  degrees  Rankine, 

^s  **  ^o  =  height  of  temperature  distribution  -  degrees  F, 

T*  «  T  (o)  =  temperature  along  z-axis  -  degrees  Rankine,  and 
cC  is  a  width  parameter. 


Under  these  circumstances,  the  sonic  velocity  will  be  given  by 


where 


C  c/t)  -  C„  |j  + 

Co  “  speed  of  sound  in  the  homogeneous  atmosphere 

p  =  *  = 


(C-2) 


Equation  (C-l)  has  been  plotted  for  Ts  =  900*F,  T0  =  60aF  in  figure  103,  and  is 
compared  with  o  typical  hot  jet  profile  of  the  same  T$  value.  The  value  afC&,  which 
controls  essentially  the  width  of  the  distribution,  is  fixed  by  demanding  that  T$  -  T0 
falls  to  half  its  value  at  the  some  distance  that  the  jet  temperature  fails  to  half  its 
original  value. 

In  the  subsequent  sections,  an  analysis  based  on  ray  acoustics  and  wave  theory, 
respectively,  is  presented  for  the  particular  variation  of  the  speed  of  sound  given  by 
Eq.  (C-2). 


*lt  is  explicitly  assumed  here  that  the  temperature  distributions  existing  in  a  jet  are  wel 
approximated  by  the  measured  velocity  profiles. 
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Assumed  temperature 
distribution  (Eq.  C-l) 


GEOMETRICAL  ACOUSTICS  APPROXIMATION 


Inasmuch  as  an  analysis  based  on  ray  acoustics  (high  frequency  approximation)  can 
be  helpful  in  gaining  some  insight  in  the  way  sound  propagates  in  a  non-homogeneous 
medium/  this  section  will  present  a  calculation  of  the  ray  paths  and  the  intensity  distribution 
for  rhe  case  under  consideration.  Reflecting  boundaries  near  the  source  will  not  be  considered 
here,  for  in  this  case  the  use  of  the  approximation  becomes  impractical. 

T':e  Ray  Paths 


Consider  a  z  -  /x,plane.  A  ray  starts  at  the  source  (assumed  located  at  the  origin) 
making  an  angle  <|>swith  the  z  -  axis;  let  the  velocity  at  the  source  be  cs (figure  104).  Denote 
he  nngfo  which  the  tangent  to  the  ray  at  a  point  ( fl ,  z)  makes  with  the  vertical  by  <|>and  the 
velocity  at  that  point  by  c  (/t). 


The  ray  path  can  be  determined  from  Fermat's  principle,  which  requires  that  the 
integral  j  ^  be  stationary,  <  w = ° 


where  ds  is  an  increment  of  the  path,  c  the  sonic  velocity,  and  ^  the  symbol  for  the  variation 
of  a  functional.  In  t^e  present  case,  c  =  c(tx),  so  that  j u  u 

U  .[(..)*•(.,)*],  [,. if  ^ 

ri  •  .  11*11 


The  variation  problem  then  becomes 


for  which  there  corresponds  the  Euler  equation 


■il  cxn.  —  o 

Cl  n.) 


r  jl_  _  j  t  ]  - 
L  ^  -jrzuyr-j 


Hence, 


Ui4 


c*  in) 


A 

r  _  f  CLin.) 

Lrzr 


j  V  -A  -  c\/0 


dLA 


where  A  is  a  constant.  Substitution  of  «£  3  =  Co1'<|>yields  A  = 

which  is,  in  fact,  Snell's  law;  the  usual  form  of  the  law  results  if  we  let 

A  =  cs/cos 

Thus,  the  equation  to  be  integrated  is 


*  V 


•y*v* 

k\* 

/j 


.  ;  •  V"* \>  .v  /•  V- 


.  »V  «-.  *\  v.  ^  ■r,  ^  »\  sV-V  V'.V, 


-  .  »*K  "  *  ■»>  ' 


V  J 


i/. 


h:: 


a 

S 


J3 

I 

is 


Alternatively,  Eq.  (C-3)  can  be  derived  by  assuming  Snell's  law  and  making  use  of  the 
geometry  of  the  situation.  Substituting  for  c2(/t)  and  c2s  =  c2(o)  from  Eq.  (C-2)  and 
simplifying,  we  obtain 


ib?  [ 


V 


C C  jn^ot zt 


Co'SKoC/L 


-y 


(C-4) 


where 


Let 

to  obtain 


CL-!±£>  ;  y  -  = 

<Zos*cf>z  *  4--1  P  + 

cdsK^ock.  +J3  =  * 


(C-5) 


f3 


fs«*. 


13 

>1 

P* 


H 

L  «’ 
£‘j 

fi 


5 


B 


L*  i 

t’i 


wa 

-~£~f 

JsaJ 


dw 


W,  W  P  wXi-  §  wKl-  6w) 


or 


Jrw4 
^  _ 

W.V 


<w 


*  (C-6) 


where 


S  =e,  VY  =  -M&±iL  :  «.  =  &  + 

0  “  j2>»v«*4,  1  ^ 


,c- 


(C-7) 


X  =  -=L^  =  iP=i- ;  8=-^5««t- 

ocvh^t  ctVP'+  ^  6 

i  .  i 

The  integral  in  Eq.  (C-6)  can  now  be  reduced  to  a  sum  of  two  elliptic  integrals  by 
a  further  transformation  of  the  variable  {Ref.  2).  Let 

if  =  .  (C-8) 

°  V.W-  S") 

to  obtain  immediately 
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*  < 
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The  integrals  appearing  in  Eq.  (C-9)  are  elliptic  integrals  of  the  first  and  third  kind, 
respectively.  For  numerical  evaluation  purposes,  they  are  now  expressed  in  the  Legendre 
normal  forms,  by  making  the  further  transformation  (Ref.  3) 

<ja  =  4-  c. (C~n) 


which  yields 

oC£  =i/sr  r.  co% 


/bf^hJ+WVKTT (*w.k) 


(C-12) 


where 


-  r  =  r  ou> 
j0Y^rap  jo  ** 

TT(i|v<-A)=  /  - — i3L _ 

Jo  (l-MVnNf) 

le*  =  ^^4 a  (*n©<lulu*)  (c_13) 

2j/  -  olrv~‘  _ V~  'SmK  cen, 

L  {(fi+  ao'b^c<:  Co*cj>.J^ 

The  integral  F(.T(Me)  can  readily  be  evaluated  for  a  given  <j>s  with  the  use  of  available 
tables  (Ref.  3);  the  integralTTQf^bJcan  be  expressed  in  a  number  of  ways,  e.g.  in  terms 
of  Jacobian  Theta  and  Zeta  functions  ;  noting,  however,  that  the  parameter  is  always  greater 
than  unity  (except  in  the  limiting  case  of  sin  0).  Use  can  be  made  of  an  addition 

theorem  relating  the  Integral  of  a  parameter  greater  than  unity  to  that  of  a  parameter  less  than 
unity  (Ref.  4): 

TTty./i.le)  =Fty,k)  -TT 


* 

+l>-oc 


“bav\lr\  '  )~1  -fa.* 


9?'i 


where  M  =  ^~/CC —  ^  ^  (Tables  of "TT^Mi^for  fXC.  \  are  available,  (Ref.  5).) 
Substitution  of  the  above  expression  in  Eq.  (C-12)  yields  the  final  result 

FW^V  Sw^JT^A) 


?=v^ 


l£±L 

CoS 

1  J£> 

1 ». 


*ta.nK_' 


VM( l-le  5m  tf>) 


(C-14) 


L  V/^O-lesurtfO 

It  is  clear  from  Eq.  (C— 13)  and  Eq.  (C— 14)  that  all  rays  emanate  from  the  origin,  and  in  the 
limiting  cases  ofrt*0^5«^othe  rays  are  straight  lines  along  the  z  -  and  rc-axis,  respectively. 

Ray  paths  for  five  values  of  have  been  plotted  in  figures  105,  106  for  =1.61  and 
^  =  5  .077,  corresponding  to  a  temperature  at  the  source  of  approximately  900°F  and  2700°F, 
respectively,  and  an  ambient  temperature  of  60°F.  As  one  would  expect,  the  rays  ere  "bent" 
toward  the  horizontal,  while  inside  the  temperature  distribution;  at  some  dlstancectft,  at 
which  ambient  temperature  has  been  virtually  reached,  they  become  straight  lines.  The 
"exit"  angles  at  this  distance  are  indicated  in  the  figure.  It  is  to  be  noted  that  the  bending 
becomes  more  pronounced  as  decreases.  The  plot  is  symmetric  about  the  (c ot)~  and  (<(.*)- 
axis,  and  the  three-dimensional  picture  results  by  rotating  the/t-fc plane  through fcfrh bout 
the  z  -  axis. 

Intensity  Distribution 

With  the  knowledge  of  the  ray  paths,  the  Intensity  distribution  can  easily  be  determined 
by  considering  the  geometrical  divergence  of  the  rays  (Ref.  6).  The  procedure  used  here 
will  be  that  applied  by  Prldmore-Brown  and  Ingard  to  the  case  of  a  temperature-stratified 
atmosphere  (Ref.  7). 

Let  Q  be  the  total  power  output  of  a  spherical  source  and  the  radius  of  a  small 
sphere  centered  on  the  source.  Consider  two  neighboring  rays  leaving  the  source  at  <^>% 
and^&Kflgure  104).  If  these  rays  are  rotated  about  the  z-axis,  a  certain  amount  of  power, 
lows  through  the  space  bounded  by  the  two  surfaces  formed,  and  it  is  equal  to  the 
power  which  falls  on  the  sphere  zone  cut  out  by  these  rays: 

•  air (a ^  (^ *&)  = 

The  Intensity,  ,  at  any  point  in  the  plane  will  then  be  the  ratio  of  this  power,  §Q,  to 
the  area  through  which  it  flows  at  that  distance  (=2e*\J),$f being  the  perpendicular  separa¬ 
tion  between  the  rays.  Thus, 
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Figure  106.  Calculated  ray  paths 
Vertical  broken  line 
temperature  distribut 


Q.  ^»  in  <K 
-4-<r  ft 


iii 

$* 


where  the  value  of  §£  Is  calculated  with  the  aid  of  figure  104. 


S  i  v\  4>  S  Z  =  Sin 


JLn~ 

djt. 


Since  we  will  b*j  interested  in  the  intensity  as  a  function  of  distance  rather  than  a 
function  of  the  angie  i>  ,  we  now  substitute  in  the  integrand 
From  Snell's  law  1 


C.QS  4>  _  c-Lk) 

Co%  C-* 


cosf4»  +  S4>)  =  cOtVSc 
C0Kt,+ $«M  C» 


Neglecting  ^  £  and  expanding  the  cosines,  one  obtains 


As 

fe-  = 

"Tivi 

s< 

K 

T**  «4> 

Also  using  Eq.  (02),  Snell's  law  reads 

r1+  *  l 

Cos  <p 

1  *  T  co^Kocn. 

eos*«^>s 

whence 

i-8*1] 

whence 


Figure  107 


Figure  i C9  .Intensity  variation  aiong  lines  parallel  to  the  axis  of  symmetry  of  the  temperature  distribution  (?  -  axis).  T  =900°F 


-  *f=  1.0 

- <*/>  =  2.0 

- V-  3.0 

- oty»  =  10.0 

- a/>=  20.0 


The  intensify  variation  along  lines  parallel  to  the  z-axis  (any/j-fplane)  or  on  the 
surface  of  cylindrical  shells  (three  dimensions)/  the  z-axis  being  their  axis  of  symmetry,  is 
shown  In  figure  109.  Again,  regions  of  relative  increase  are  found  in  the  intensity. 

The  polar  distribution  of  the  intensity  is  shown  in  figure  110,  where  the  angle  is 
that  which  the  axis  of  an  observer  makes  with  the  z-axis.  A  slight  shift  of  the  maxima 
toward  =  90°  is  noted  as  an  observer  moves  away  from  the  source.  Finally,  a  comparison 
of  the  directivity  pattern  atcCp  =  20  forjES  =  1.61  and  Jb  =  5.077  indicates  the  increase 
in  the  directionality  of  the  source  for  larger  values  of  ft  (figure  111). 


WAVE  THEORY 


The  Formal  Solution 


In  this  section  the  solution  to  the  wave  equation  will  be  sought  in  the  form  of  a 
Green's  function  expressed  as  a  series  of  normal  modes.  To  do  this,  a  plane  boundary  - 
acoustically  described  by  a  normal  impedance  Z  independent  of  the  angle  of  Incidence  - 
is  introduced.  The  boundary  Is  defined  by  y  =  constant  or  by/c  =  Jl0/ $tn  Q  where 
H. )  6  are  the  cylindrical  coordinates  and/t0  is  the  perpendicular  distance  from  the 
origin  (figure  112).  ^©<0<C'tT) 

Consider  a  spherical  source  of  sound  located  at  the  origin  of  coordinates  and  driven 
at  a  constant  frequency  —  Cu/S.'tf  Th®  sound  field  at  a  point  pfy t,0,?)may  be 
derived  from  a  velocity  potential 


-I  cot 

z)e 


satisfying  the  wave  equation 


v\|> 


aj 


0*00) 


^  = 


(C— 23) 


The  use  of  this  equation  entails  two  inherent  approximations,  namely,  the  usual  require¬ 
ment  of  small  amplitudes  and  that  the  fractional  change  In  any  property  of  the  medium 
must  be  small  within  a  distance  of  a  wavelength.  The  latter  condition  restricts  the  solution 
to  high  frequencies  but  has  not  been  found  stringent  even  for  very  high  temperature 
gradients  (Ref.  7). 

Assuming Eq.  (C— 23)  can  be  separated  in  cylindrical  coordinates, 
yielding 
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£e-  =  -le* 


(C-24) 
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where  k  Is  the  separation  constant. 


If 


-  k*j  R  = 


(C-25) 


(i)R  and  its  derivative  are  continuous  (so  that  the  normal  model  theory  holds), 
(il)  at  the  plane  boundary  the  logarithmic  derivative  satisfies  the  condition 


R 

1 


i  Ui 


•f 


(ili)c6^is  a  continuous  monotonic  function  of  hC  such  that  c(o)  =  cs  and  c 

where  /t  =  -H/Sifl  9  defines  a  "pressure  release"  boundary  to  the  left  of 
the  source  lY^Q^aih) 

(iv)  for  large  values  of  /t,  the  solution  represents  only  outgoing  waves. 


a  set  of  eigenfunctions with  corresponding  eigenvalues  may  be  determined. 
The  space  part  of  the  solution  will  then  be 

(C— 26) 


i|)  C^.z)=V  Aw<£ 


where  A„,  are,  as  yet,  undetermined  coefficients,  and  where  the  solution  to  Eq.  (C-24), 

*  _iw& 

C  in  conjunction  with  the  time  factors  represents  outgoing  waves,  the  positive 

or  negative  sign  being  used  in  the  region^  >  o  or  Z  CO  respectively. 

To  evaluate  the  coefficients  A,,,,  substitute  Eq.  (C-25)  into  the  inhomogeneous  wave 
equation 


^  SM  Sc*) 


(C-27) 


where  the  right  hand  side  corresponds  to  a  point  source  located  at  the  origin.  Here  ^ 
denotes  the  Dirac  delta  function  defined  by: 

S  M  *  o  )  t  =£  o 

=  •  X  =.  o 


A 


in  such  a  way  that  where  the  region  of  integration  includes  the  point  x  =  0. 
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Substituting  Eq.  (C-26)  in  Eq.  (C— 27)  and  using  Eq.  (C-25),  we  obtain 


e  4*  k 


Now,  multiply  both  sides  byi^  ,  integrate  over  a  small  cylinder  enclosing  the  source 
and  letg+Q.  Utilizing  the  orthogonality  relation  for  the  eigenfunctions,  namely 

>*/*  — H/smfi 


*•  E«0t.)  (£„£*•)  d*.  = 

*  *■  /t*/ %\kQ 

we  can  obtain  the  following  expression: 
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where  the  functions  depend  on  the  velocity  profi  le  £(*)  and  the  boundary  condition 
at  the  plane  boundary. 


The  Radial  Equation 


With  the  particular  choice  of  C  (rt.) given  by  Eq.  (C-2),  the  radial  equation 
(C-25)  becomes: 


where  the  approximation 


(  j  4.  — ^  | _ — 

\  c°4>hcCK  •  C.O<>b?cLHs 

will  hold  for  p  44  CO^hot/C  which  in  general  will  be  true. 

Introducing  a  new  dependent  variable  LC  related  to  by 

*R  = 

one  obtains  the  following  equation  in  : 
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(C-32) 


for  1  —  C^/lu^  )  >  O 

-+•  f(-  -}-*^k*-r^-T-'^L*T^- s"'h“n)^ 


where  F  stands  for  "hypergeometric  function"  and 

e=  [-^0-^)  f 

K=  [g  0  "  ] 


(C-33) 


£  ^  c  2.  =  arbitrary  constants. 

The  hypergeometric  functions  appearing  In  Eq.  (031)  and  Eq..  (032)  are  defined 
in  the  circleSi»li**4.lin  which  they  converge.  By  analytic  continuation,  the  solutions 
may  be  extended  to  the  region  of  convergences**^*.?!  /  ^us: 


?oO  = 


c  [  r&yrt-e) _ 

‘l  r(-£-f-)r  <&♦*-£) 


4.  r(t)rc^ - ...  st«vVv,i«..p(-| ~ 

"5” 

+  a'  TV  i  'n 

+  mEM _ .s.-kT*'*-.  Ftr&'i'-iS  _rr 

♦DfcF*F+i5  v  ^ 


[  |  _  (k’lC^/O;  <  O  j  1 


1C  34) 


-»  V 

j,  V 


\y 


Using  the  definition  of  £  and'K,  we  finally  obtain  the  asymptotic  solution  to 
the  radial  equation  In  the  form 


T?(a) 


ft.-**? 


r^te4  ^r-  )0-^-<O) 


T?t*terc*e^  a,;  4^-  iO-J&ge->o) 

“  V  A. 


where!* Ci” C^*'R./^»S)  and  where  the  prime  {  '  )  indicates  substitution  of  by 
in  the  definitions  of  Eq.  (C-36).  The  constant  C2  can  be  evaluated  from  the  boundary 
conditions  at  the  source  and  will  depend  on  the  source  strength. 


The  Eigenvalue  Equation 

The  eigenvalues,^^  corresponding  to  the  eigenfunctions  Rm(K  )  are  determined 


from  the  condition  at  the  plane  boundary,  namely  that 


_  icoj3  at  A.  =  M.,/ 


(C-38) 


in  conjunction  with  the  defining  equations  for  d:  or"K  (Eq.  C-33).  The  eigenvalues 
resulting  from  the  solution  of  Eq.  (C-38)  will,  in  general,  differ  depending  on  the  region 
of  convergence  considered,  i.e.  on  the  choice  of  Rm  to  be  used  in  Eq.  (C-38),  and  can, 
in  principle,  be  computed  for  a  particular  mode  by  retaining  the  appropriate  number  of 
terms  in  the  hypergeometric  series  which  appear  in  the  solutions  Rm. 


Consider  the  region$i«kAft^iand  the  case  of(l-j*^*|CO,  for  which  the  appropriate 
eigenfunctions  are  given  by  Eq.  (C-31). 


With  the  use  of  Eq.  (C-29)  and  the  identification 
S\aV\'<C*-  s  £ 
s  a. 

=  b 

si  C 

the  solutions  take  the  form 


”*  4: 

1 

2. 


'pC*.')  =^c»*woC>t|c,FCat*>iCj -*5)  +  C^§" 
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which  when  substituted  into  Eq.  (C-38)  yield  (for  a  perfect  reflector)  the  following 
equation  in  £  ; 


(5aC 


**  J|  C* 


rc*y 


^  Roo  ret) 


rc^^)  Hbvn)  X. 

r*  cc+^) 

r Cc4 1) 

+v?  rRry  r(eV  j 
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rtc*o  a\  \ 

where  the  ratio  C]/C2  is  given  by  Eq.  (C-37)  and  where  the  relations 

-J-F(a6 ,JV>^2)  =r  F(0C  +  l^4.|,CVl,^) 

2  rc-z)  -  ru4i) 


(C— 39) 


and  the  definition 


P(<t,a  vt)=  re?)  V’  rcctf. 
r  ’  ’  ;  r&t)  v(p)  r 


Virus*") 


C'v*^ ) 


w  ! 


hove  been  used.  The  evaluation  of  Eq.  (C-39)  is  to  be  carried  out  ot/t=  A  similar 

yet  more  complicated,  equation  results  If  Eq.  (C-34)  Is  substituted  In  Eq.  (C-38)  In 

T,'  °Jc'0i,dJ‘3m  s°lu,l<>"  e  (  h<>"c.W)  f«  ooy  jlveo  mode  coooo.  be 

easily  obtained  by  ordinary  means,  except  in  limiting  cases. 


Such  a  case  results  if  we  consider  distances  from  the  source  for  which  the  asymptotic 
limit  of  the  hypergeometric  functions  can  be  used  to  a  fair  approximation*.  Consider,  for 
example,  the  regionSnlttCA)).  Substituting  unity  for  the  hypergeometric  functions  in 
Eq.  (C-34)  and  using  Eq.  (C-37),  the  eigenfunctions  take  the  simple  form 

"IhPC/t)  =  /C*  Cosh  C  4~ 

whence 

so  that  Eq.  (C-38)  yields 

~  C“0^'  “^CC  *U.  - 

or 


where  "U.  s  <X.K0/^\t\  Q  (<3^  0  C  "O') 


Eq.  (C-40)  relates  the  eigenvalue  ,  appropriate  for  the  assumed  limiting  case,  to 
the  temperature  parameters  06  and  JS  the  frequency  ■£-  uj^,  and  the  boundary  impedance 
2:  •  A  si  ml!  jr  expression  to  the  above  Is  obtained  in  the  case  of  ,  . .  by  substituting 

iK/A  f“  ^  •  0-^)>° 

Thus  far  in  this  section  the  wave  equation  has  been  formally  solved  (Eq.  C-28),  v,nd 
the  solution  to  the  radial  equation  has  been  obtained  in  a  certain  region.  Its  validity  de¬ 
pends  on  the  particular  choice  of  C^fand  £  values  for  which  the  approximations  evoked 
in  solving  the  radial  equation  are  valid.  An  equation  satisfied  by  the  eigenvalues^ has 
also  been  obtained  (Eq.  C-39).  However,  before  the  above  results  can  be  put  to  practical 
use  the  eigenvalues  appropriate  for  the  region  considered  must  be  determined  in  a  closed 
form,  for  it  is  precisely  the  knowledge  of  %*which  will  Indicate  how  the  sound  is  propagated 
in  the  medium.  This  will  not  be  carried  out  here.  Instead,  any  qu-litatlve  conclusions 
which  are  to  be  drawn  regarding  the  effect  of  the  presence  of  the  temperature  distribution 
considered  in  this  appendix  are  to  be  based  on  the  results  of  the  analysis  presented  in  the 
section  Geometrical  Acoustics  Approximation,  and  will,  of  course,  be  valid  only  In  the 
high  frequency  range. 


aximatlon  one  order  higher  than  that  In  which  only  the  first  made  is  assumed 


to  survive. 
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LIST  OF  SYMBOLS  USED  IN  APPENDIX  C 


A,  oL  /  J& 

E  (f ,  k) 

F  (Y  /  k) 

F  (c,b,c,x) 

I 

k 

Q 

r,  0  ,  z 
T 

W,  q, 

Z 

F(x) 

$ 

S(x) 

Tf(y,/*,k) 


Temperature  distribution  parameters  (Eqs.  C-l,  C-2) 

Speed  of  sound 

Elliptic  Integral  of  the  second  kind 
Elliptic  integral  of  the  first  kind 
Hypergeometric  function 
Acoustic  intensity 
Modulus 
Acoustic  power 

Radial/  circumferential  and  axial  coordinates,  respectively, 

of  polar,  cylindrical  coordinate  system 

Temperature 

Transformation  variables  (Eqs.  C- 5,  C-3) 

Impedance 
Gamma  function 

Variation  of  a  functional;  also,  constant  (Eq.  C-7) 

Dirac's  delta  function  (Eq.  C-27) 

Parameter 

Elliptic  integral  of  third  kind 


Superscripts 


Subscripts 


Reciprocal 

Time  rate  of  change 


Unperturbed  or  ambient 
At  the  source 


Definitions 
the  text. 


of  other  symbols  used  in  algebraic  manipulations  may  be  found  in 


APPENDIX  C 
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